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ABSTRACT 


Man has always been associated with soil. The responsibility 
of designing foundations for structures has forced the 
Geotechnical Engineer to study the suitability of the soil in 
taking the load from the super -structure. Often it becomes 
imperative to improve the soil properties to prevent its failure 
under the load. Reinforcing the ground with inclusions in the 
form of steel, bamboo, nails, fibres etc., helps in increasing the 

load carrying capacity of the soil- The last decade has seen the 

« 

advent of geosynthetics in the form of geotextiles, geomembranes, 
geogrids, geocells, etc., and their varied applications. Geo- 
synthetics are used to design steep slopes and retaining walls, 
and for improving the performance of the soi 1 -foundat ion systems. 
The need for designing the reinforced soil rationally has kindled 
research activity in this field. 

While a large number of model test results have been 
reported, only few analytical models are available to study the 
performance of reinforced foundation beds. Most of the analytical 
studies deal with the increase in the bearing capacity of the 
system whereas very few of them deal with the reduction of 
foundation settlements due to reinforcement embedded in the soil. 
Further, the interaction between the soil and the reinforcement, 
has not been studied fully or rationally in many of the analyses. 
In the present work, an attempt has been made to analyse the 
soil -strip reinforcement interaction in order to compute the 
interfacial stresses and eventually, the settlement reduction of 
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the surface, using the elastic continuum approach. The membrane 
effect of the reinforcement, however has not been considered. 

The mechanisms that contribute to a reduction in the surface 
‘juttlements may be enumerated as follows. The applied surface 
load causes the soil below it to move downward and outward. The 
r e inforcement , placed at some depth below this load, experiences 
vertical deformation as well as axial extension. Depending on its 
vertical and horizontal stiffnesses, the reinforcement resists the 
soil displacements and consequently, normal and shear stresses are 
mobilised at the interface. These stresses, in turn push the soil 
up thus bringing about a surface heave which eventually is the 
reduction in the surface settlements. In the thesis, the rigidity 
of the strip reinforcement in the horizontal (Chapter 3) and 
vertical (Chapter 4) directions is initially considered 
separately. Following the analyses for rigid strips, the same is 
extended (Chapter 5) to extensible and flexible ones. For high 
values of the respective stiffnesses the results compare well with 
those obtained for rigid strips. The combined effect of the 
horizontal and vertical rigidity of the strips is then considered 
and allowance made for possible slip along the strip -soil 
interface. Axi -symmetr ic problems for rigid and extensible sheets 
are studied and reported in Chapter 6.^ 

The investigation is carried out for a strip reinforcement 
placed at a given depth, centrally below an uniformly loaded 
rectangular area. The Boussinesq'^s equation for displacements 
within a soil mass due to a point force on the surface is 



numerically integrated over the loaded area to calculate the 
horizontal and vertical displacements of points along the length 
of the strip. The reinforcing strip resists these displacements. 
Shear and normal stresses are mobilised at the soil -strip 
interface. In chapter 3, the horizontal displacements of these 
paints due to the shear stresses are computed by integrating 
numerically the Mindlin's solution for horizontal displacements 
due to a horizontal force within an elastic continuum over the 
area of the reinforcing strip. Satisfying the compatibility of 
displacements at the points, a set of simultaneous equations are 
obtained, the solution of which yields shear stresses mobilised 
along the strip. For a rigid ( inextensiblel strip the net 
horizontal displacements are zero whereas they are equal to the 
elongations of points along the strip for an extensible one. For 
the vertical displacements of the points along the strip (Chapter 
4) due to normal stresses on it, Mindlin^'s solution for vertical 
displacements due to a vertical force within an elastic continuum 
is integrated over the area of the strip- “The compatibility 
condition for a rigid strip is that the net vertical displacements 
of points in the soil adjacent to the strip are equal to the rigid 
body translation of the strip. For a flexible strip the soil 
displacements arc equal to the flexural deformations of the strip- 
For the case where the rigidity of the strip in the horizontal and 
vertical directions are considered simultaneously (Chapter 5), the 
effects of the vertical stresses on horizontal displacements and 
vice versa, are also included. Slip along the strip is 
incorporated by limiting the shear stresses mobilised to the 
frictional resistances available at the -interface due to the 


normal stresses including the overburden stresses at that depth- 
The reduction of settlement of points on the surface due to the 
shear stresses along and the normal stresses on the strip, are 
computed by integrating Mindlin'^s equation for vertical 
displacements due to horizontal and vertical forces respectively, 
over the area of the strip. A similar analysis is presented 
(Chapter 6) for an axi -symmetr ic case. In the case of strip 
reinforcements, pairs of strips are also considered and the 
appl icabi 1 ity of the principle of superposition examined. 

A parametric study brings out the effects of the aspect ratio 
of the loaded area, length, depth and the distance of the strip 
from the centre, and the Poisson^'s ratio of the soil for strip 
reinforcements, on the stresses mobilised at the interface and the 
reduction of settlements of points along the surface. The effects 
of the radius of the extent of reinforcement and depth of its 
placement for the axi -symmetr ic case, on the quantities mentioned 
are also studied. The predicted reduction in surface settlements 
are compared with those from selected test results available in 
1 i terat ure. 

The elastic continuum approach is thus used to model the 
complex problems of reinforced foundation beds. The analysis 
quantifies the reduction in foundation settlements due to strip 
and circular sheet reinforcements. An effective method to predict 
the reduction in surface settlement due to strip and circular type 
reinforements placed beneath a loaded area is developed. 
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CHAPTER 1 
INTRODUCTION 


From time immemorial, man has been associated with soil. 
Any civil engineering structure, whether a compound wall or a 
multistorey building has to be founded in or on the surface of the 
earth. Foundations have to be designed with a proper understanding 
of the soil characteristics. There are situations wherein the 
soil is extremely weak or soft to take any load. Design of 
foundations on soft soils involves heavy expenditure as this 
requires large foundation dimensions. Embankments on soft soils 
call for flatter slopes, thus occupying large areas which may not 
be available. Many a times it may be required to construct steep 
slopes for highways. Further, unpaved roads on soft soils require 
a thick course of aggregate spread over a large area so that the 
wheel loads are distributed over a wide area. This again would 
result in an uneconomic design. 

The increased cost of these conventional designs and 
numerous environmental constraints greatly encourage the in-situ 
improvements of weak subsoils. Environmental constraints include 
the depletion of supply of sites with good foundation conditions 
due to rapid urbanisation. Under such circumstances it becomes 
imperative to enhance certain soil properties, thus resulting in a 
significant improvement in the bearing capacity of the soil. To 
achieve this, various ground improvement techniques are in vogue 
which may be listed as follows. 
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1 . Dynamic Compaction, wherein the soil is compacted and 

to consolidate by repeated dropping of a heavy weight on the 
surface of the soil 

2- Installation of sand drains to expedite the consolidation 
process 

3. Use of lime piles to increase the shear strength of the soil 
and reduce its compressibility 

4. Use of stone columns/Granular piles 

5. Chemical stabilization 

6. Thermal stabilization 

7. Consolidation by preloading 

8. Heating and freezing 

9. Soil Reinforcement 

Soil Reinforcement is one of the latest and fast growing 
techniques of ground improvement. Although the systematic study 
of reinforced soil did not start until very recently, the practice 
of building houses and roads on fibre-reinforced soil and 
constructing earth walls with different types of reinforcing 
inclusions is an age old art. Use of rope fibres and bamboo 
strips to strengthen rural base roads has not been uncommon. 
There still is a practice of using rectangular grids of bamboo 
strips as the central core around which mud walls are built in 
rural areas. 

It was not very late when the potential of soil reinforce- 
ment was realised. The infrastructure for systematic study of 
reinforced soil was only lacking until Vidal (Jones, 1985) gave 
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the geotechnical coinmunity a basis for a rational approach to 
reinforced soil design. Research in this area started picking up 
momentum and as a consequence, various types of reinforcing 
methods and forms have emerged. The major structural effects of 
reinforcing inclusions can be tension, compression, bending and 
shearing- Broadly, reinforcements can be classified into three 
types, viz., vertical, inclined, and horizontal reinforcements. 

Reinforcement techniques such as stone columns and 
micropiles can be categorised as vertical reinforcements. A stone 
column is a vertical column of highly compacted sand, gravel and 
aggregates. It acts as a ''granular pile^ and its role is to 
increase the resistance and stiffness of the foundation soil. The 
merit of a stone column is its ability to adjust to the applied 
load and thereby redistribute the load where stress concentration 
occurs. Unlike piles, the mechanism of interaction is that of a 
restrained expansion in the surrounding soft soil. The inherent 
advantage is that there is no collapse because an overloaded stone 
column automatically gets relieved of the stresses as it deforms. 

A micropile is made of a bar or a tube of a few centimeters 
diameter surrounded by grout all along its length. Micropiles 
help in increasing the factor of safety of a slope against sliding 
and transmit the shear force from a possible sliding mass to the 
stable soil. Inclined reinforcements are used to stabilise slopes 
and embankments. Soil-nailing, a recent technique for stabilising 
slopes is an example of the same. This method is based on the 
principle of reinforcing soils by means of tension carrying bars 
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(so called 'nails'). Soil-nailing is being used at present to 
stabilise natural slopes, cuts and excavations in granular soils 
with some cohesion, stiff clays and soft rocks. 

Horizontal reinforcements in the form of strips and sheets 
are routinely used in almost all soil -reinforcement projects. 
Sheet reinforcement used in highway construction helps in reducing 
the thickness of the required subgrade. Strips are also used to 
reinforce embankments and retaining walls. Strips, as well as 
sheets are used under footings to increase the bearing capacity of 
the soi 1 — f oundat ion system and to bring about a reduction in 
foundation settlements. The advent of Geosynthetics has resulted 
in their extensive application to soil reinforcement. Geosyn- 
thetics in the form of geotext iles, geomembranes, geogrids, 
geomats, etc. have found varied applications in geotechnical 
engineering. Geosynthetics used in pavement design help in 
redistr ibut ing the wheel loads over a larger area and hence the 
thickness required for the subgrade is reduced. 

The basic mechanisms of reinforced soil can be explained as 
follows. When a sample of cohcsionless soil is loaded vertically, 
it experiences lateral strain. However, if horizontal reinforcing 
elements are placed within the soil mass, the reinforcements will 
prevent lateral strain because of friction between the reinforcing 
elements and the soil, and the behaviour will be as if a lateral 
restraining force or load had been imposed on the element. Thus 
it is evident that the reinforced soil concept is essentially 
based on the mobilisation of the inter facial shearing resistance 
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between the 5oil and the reinforcement which in turn restrains the 
lateral deformation of the soil. Beneath a footing the 
reinforcement checks the lateral displacement of the soil by 
introducing lateral confinement. Further, the footing load also 
causes vertical movement of the soil. If the soil is reinforced 
with strips with some flexural stiffness, these strips will resist 
the vertical deformations because of the mobilisation of the 
intcrfacial normal stresses between the soil and the reinforcing 
strip. These stresses, in effect push the soil mass above upward. 
In the earlier case, the shear stresses at the interface push the 
soil upward. It can hence be concluded that shear and normal 
stresses are mobilised at the soil -reinforcement interface as a 
consequence of its rigidity in the axial and vertical directions 
which in turn result in the reduction of surface settlements. 

The present work aims at computing the stresses mobilised 
along the soi 1 -reinforcement interface and hence to predict the 
settlement reduction of points along the surface. Strip 
reinforcements are considered below rectangular loaded areas and 
circular reinforcements below circular loaded areas. The elastic 
continuum approach is resorted to model the reinforced foundation 
system. The following chapter reviews the literature available 
on reinforced foundations. 

Chapter 3 deals with the rigidity of strip reinforcement in 
the horizontal direction. The strip is considered to be 
inextensible. The horizontal displacements of points along the 
soil -strip interface due to the surface load are calculated using 
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Boussinesq's equation (Poulos and Davis, 1975J for radial 
displacements at any point within an elastic half-space due to a 
point load on the surface. The strip beinq inextensible 
counteracts these displacements and hence shear stresses are 
mobilised at the strip-soil interface. The displacements of the 
points along the strip due to these stresses are evaluated using 
Hindi in's equation (Poulos and Davis, 1975) for horizontal 
displacements due to horizontal force within an elastic continuum. 

The shear stresses are computed by equating the net horizontal 
displacements of the points along the strip to zero. The 
settlement reduction of points along the surface due to these 
stresses are calculated using Hindi in''s equation for vertical 
displacements due a horizontal force within the elastic continuum. 

The rigidity of the strip in the vertical direction 
(inflexible strip) is considered in chapter 4. In this chapter the 
vertical displacements of points at the soil -strip interface are 
evaluated using Boussinesq'^s equation for vertical displace- 
ments. The strip being inflexible experiences a uniform rigid 
body translation and as a result normal stresses are mobilised at 
the interface. These stresses are the difference of the normal 

/ 

stresses acting on the top and bottom faceC. of the strip- The / 
vertical displacements due to these stresses are computed using 
Hindlin^'s equation for vertical displacements due to a vertical 
force at depth. The compatibility of displacements at the 
interface yields equations for the normal stresses and the rigid 
body displacement which are solved for the same. The settlement 
reduction of points on the surface due to these stresses are then 
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Following the analyses for rigid strips, the same is 
extended (Chapter 5) to extensible and flexible ones. For an 
extensible strip the net horizontal soil displacements arc equated 
to the elongation of the points along the strip. For a flexible 
one the net vertical soil displacements equal the flexural 
deformations of the points along the strip. The shear stresses 
(for extensible strips) and normal stresses (for the flexible one) 
are computed from the equations obtained on satisfying the 
displacement compatibility and force equilibrium. The coupled 
for the horizontal and vertical rigidity of rigid strips is 
also considered and allowance made for possible slip along the 
strip-soil interface. In the case of slip, the shear stresses are 
restricted to the strength available at the interface. 

Using the elastic continuum theory, an analysis for 
predicting the reduction in surface settlements due to the 
presence of a circular reinforcement at depth is presented in 
chapter 6. The problem considered is axi -symmetr ic in nature. 
The effects of the radius of the reinforcement, and depth of its 
placement on the shear and normal stresses mobilised along the 
reinforcement -soi 1 interface and the subsequent reduction in 
surface settlements are studied. The horizontal and vertical 
rigidity of the reinforcement sheet are considered separately. 
Subsequently, the extensibility of the sheet is also accounted 


for . 



CHAPTER 2 

LITERATURE REVIEW 


2.1 Introduction 

Improving the soil properties by reinforcing it has been in 
vogue from early times and in particular from the time Vidal 
introduced this concept to Geotechnical Engineering. Research in 
this area to study the mechanisms of reinforced soil has been 
progressing with great momentum over the last few years. 
Innumerable model tests using various forms of reinforcing 
materials have been conducted to have an insight into the 
reinforcing mechanisms of these materials. Some analytical 
approaches have been presented based on limit equilibrium methods 
and earth pressure theories. Finite element analysis has helped 
in quantifying the increase in the bearing capacity of reinforced 
foundation beds. Very few models based on elastic theories are 
available in literature. The present chapter gives a broad 
insight into the state-of-the-art of reinforced foundation beds, 
the various methods of analysis available and the mechanisms 
proposed by various authors. Since the present study uses the 
elastic continuum theory, the next paragraph gives a short 
description of the elastic solutions available. 

Boussinesq has given equations for calculating the stresses 
and displacements of a point within an elastic half -space due to a 
point load acting on the surface. Stresses within the elastic 
continuum for various types of surface loading are available 
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CPoulos and Davis, 1975; Teferra and Schultze, 1988). Mindlin has 
given equations for the stresses and displacements of points 
within an elastic half-space due a point force acting either 
vertically or horizontally within the continuum. The soil in both 
these cases (Boussinesq or Mindlin problems) is assumed to be 
linearly elastic, homogeneous, isotropic, and semi -inf inite . 

This chapter deals with the available literature on 
reinforced foundation beds. The various aspects considered are 
the theoretical approaches, and model and field test results. 

2.2 Theoretical Solutions 

2.2.1 Analytical Solutions 

Binquet and Lee (1975a) were the pioneers in carrying out an 
analytical study on the bearing capacity of reinforced soil beds. 
Using the Boussinesq''s stress distribution they have calculated 
the shear stress distribution on the plane of reinforcement and 
hence the maximum tensile stress in the reinforcement. They have 
also defined the locus of the maximum shear stresses with depth. 

The distance, x , along the length of the reinforcement where the 

o 

shear stress is maximum, is a function of the depth of 

reinforcement and the soil properties. They have suggested values 

of maximum shear stress as 0.3q at a depth of 0.25B and around 

o 

0. Iq at a depth of 2B, where q is the intensity of the applied 
o o 

loading and B, the width of the footing. The normal force on the 
reinforcement due to the loading is obtained by integrating the 
vertical stress over the area of the reinforcement. The authors 
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have considered three failure mechanisms, viz., pull-out and 
tension failure of the reinforcements and shear failure of the 
soil above the reinforcement. They assume that the tension in the 
reinforcement is inversely proportional to the number of 
reinforcing layers. Based on this assumption they have given 
expressions for the tension as a function of the shear stresses 
and normal stresses at that depth, and the BCR, defined as the 
ratio of the average contact pressures for the reinforced and 
unreinforced soil. 

Giroud and Noiray (1981) have made an analysis of the 
unpaved road behaviour with and without geotextile as rcinforce- 
n*6nt. They have considered the effect of geotextile in increasing 
the bearing capacity of the soil. This method is based on the 
consideration of the ''tension membrane'' effect for the case of 
plane strain reinforced unpaved road deforming under the action of 
a single application of dual wheel load. Assuming an allowable rut 
depth, load spread angle, suitable geometric configuration and 
firm anchorage of the reinforcement, the strain in the 
reinforcement is calculated. On knowing the strain, the tension 
in the reinforcement is easily computed. The authors have given 
design charts for the reduction in road base thickness, due to 
reinforcement effect of geosynthetic for various traffic 
intensities. 

Ingold and Miller (1982) have proposed a theory for bearing 
capacity of footings on reinforced clay which is based on an 
equivalent undrained strength of reinforced clay. Number of 
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layers of reinforcement spaced at a distance^ S, are considered- 
The vertical stress at any depth is calculated using a 2(v):l(h) 
stress distribution* Considering the force equilibrium of an 
element of soil spliced between two layers of re inf orcement they 
arrive at an expression for the equivalent undrained shear 
strength as 

C' = + D/2)/(2S(2 + n))j (2.1) 

where is the undrained shear strength, a is the adhesion 
factor , D is the depth to the bottom most layer of reinforcements 
and B is the width of the footing. The proposed theory has been 
verified by conducting model footing tests on clay reinforced with 
geogrids under plane strain conditions. The test results show 
close agreement with the predicted response for S/B, ranging from 
0.4 to 0.8. It is observed that BCR increases significantly with 
the increase in the number of reinforcing layers and decreases 
with increase in the depth ratio, d/B, d being the depth to the 
top layer of reinforcement. 

Based on their test results and the observed failure modes 
for the reinforced ground, Huang and Tatsuoka (1988) have 
presented an analysis to predict the increase in the bearing 
capacity of reinforced ground. The failure modes observed have 
been described in the section dealing with model tests. They have 
obtained expressions for the BCR based on the failure of the block 
mentioned with respect to the failure modes. The predicted values 
of BCR match well with the test results. The results indicate an 
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increase in BCR with increase in length o-f strips and the number 
of layers upto 3. The authors have defined a term called the 
Covering Ratio (CR) which is the ratio of the width of a strip (W3 
to the sum of the width (W) and the horizontal spacing (S) between 
two strips. Their results indicate that for CR>0.9 failure occurs 
according to mode 1. As CR increases the BCR also increases. 

Madhav and Ghosh (1988) have presented two and three 
parameter models to quantify the effect of tension in the membrane 
on settlements. They have also presented solutions for embankment 
and axi -symmetr ic loadings. Their results indicate that in a two 
parameter model, greater values of tension in the membrane reduce 
the settlements below the centre of the embankment. It is further 
seen that, in a three parameter model the maximum settlement 
reduces and settlements tend to become uniform with increasing 
tensile force in the membrane. Further, an increase in the 
stiffness of the fill reduces the maximum settlement. 

Bridharan et. al. (1988) have presented a method of analysis 
for a reinforced sand mattress on soft soil. In their analysis 
they assume that the total load carried by the reinforced system 
is composed of two components, one carried by the unreinforced 
soil and the other by the reinforcement due to frictional 
interaction. The load taken by the reinforcement is a function of 
the normal force acting over it due to the applied loading and the 
overburden pressure. The total load is expressed as 


P = P + P 

8 r 


( 2 . 2 ) 



where is the load taken by the unreinforced soil and P is the 

r 

load taken by the reinforcement. P is considered proportional to 

r 

the BCR, the area of the reinforcement , the normal force and the 
interfacial friction angle. To examine the validity of their 
analysis, the authors carried out model footing tests. Saw dust 
was used to model soft soil. Their test results indicated that 
there was practically no increase in the load carrying capacity of 
saw dust due to the placement of reinforcement since the 
inter facial friction angle was very low. Provision of a sand 
mattress over the geotextile increased the bearing capacity 
considerably. There was good agreement between their predictions 
and test results. One of their results also shows that the 
geotextile performs more as a separator than as a reinforcement if 
provided at the interface. This result seems true because they 
had used a locally available woven geotextile. The performance of 
a geotextile largely depends on its stiffness and it so turned out 
that the fabric they used was well suited for separation than for 
reinforcement. Further, an inclusion in any form at the interface 
improves the performance of the soil -foundat ion system and as a 
result an increase in the load carrying capacity of the reinforced 
system was observed. 

Madhav and Poorooshasb (19891 have proposed a model consist- 
ing of Pasternak shear layers, Winkler springs and newly proposed 
rough membrane to represent the mechanical response of a granular 
fill-geosynthetic-soft soil system. Satisfying the force 
equilibrium, coupled differential equations for displacement of 
soil and tension in the fabric are arrived at. These equations 
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are solved numerically and iteratively. Their results show that 
the settlements of the soil decrease and tend to become uniform 
with increasing rigidity of the shear layer. As anticipated, the 
tension is maximum at the centre of the fabric and reduces with 
increasing distance. The results based on the coefficients of 
friction on the top and bottom of the fabric indicate that it is 
advisable to place the geosynthetic within the granular fill 
rather than laying it directly on soft ground. The influence of 
the geofabric in reducing settlements increases with the load 
intensity. The results obtained from this model compare well with 
reality. 

Houlsby and Jewell tl990) have presented an analysis for 
reinforced unpaved roads in which tension -membrane effect is not 
considered since it is insignificant at small rut depths. In this 
analysis, the loaded area is treated as axi -symmetr ic. It is 
assumed that below the circular footing active pressures are 
developed as the fill tends to spread out laterally. Passive 
conditions are generated outside the footing. Based on these 
assumptions the minimum average shear stress at the fill-clay 
interface is calculated. The presence of a reinforcement at this 
interface prevents the subgrade from taking this inward stresses 
upto the limit load. The authors comment that although there is 
an increase in the bearing capacity for reinforced road as 
compared to an unreinforced one, the increase is not simply by a 
constant factor as suggested by other authors. An important 
feature of this analysis is that the reinforcement, unlike in 
previous design methods is perceived to carry the outward shear 
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stresses which would otherwise be transfnittcd to the subgrade. 
These stresses automatically put the reinforcement into tension. 
This model seems quite appropriate and simpler as emphasis is not 
placed on the role of the reinforcement acting as a curved tension 
membrane. They have presented design charts which demonstrate 
that for certain combination of parameters, enhancement of bearing 
capacity can be achieved by the use of reinforcement. The required 
reinforcement force can also be determined from the charts. 

Sellmeijer (1990) has presented a model for the behaviour of 
a soi 1 -geotext i le-aggregate system, based on combined membrane 
action and lateral restraint. The aggregate behaviour is modelled 
by elasto-plastic shear theory, the geotextile by membrane action 
and lateral restraint, and the subsoil by its bearing capacity. 
The vertical and horizontal equilibrium and the compatibity of 
displacements for both the aggregate and geotextile are satisfied. 
Consequently, differential equations are derived for the 
geotextile deflection and tension. The novelty of this model is 
that the lateral restraint and membrane action of the geotextile, 
the concept of mobilised friction and the bearing capacity of the 
subsoil, are used. 

2.2.2 Finite Element Models 

Brown and Poulos (1981) have demonstrated how a finite 
element model of reinforced earth can be used to investigate the 
increase in bearing capacity and stiffness of a foundation due to 
placement of reinforcement in the soil. The model examines the 
effect of reinforcement on the load -settlement behaviour of a 
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strip foundation. The analysis considers the various components 
of the reinforced soil structure separately and incorporates an 
^ soil model obeying Mohr-Coulomb failure criterion. 
The reinforcement strip is treated as elastic with zero flexural 
stiffness transmitting axial forces only. The model accounts for 
slip at the reinforcement -soi 1 interface governed by a Mohr- 
Coulomb criterion. Two cases, a footing on a homogeneous found- 
ation, and a footing on a foundation incorporating a cavity or a 
pocket of extremely soft material are studied. The results of the 
finite element study show similar patterns of behaviour to the 
laboratory model test results. The reduction in the bearing 
capacity of the unreinforced footing due to the presence of cavity 
from the finite element analysis agrees well with experimental 
observation. The analysis shows that the quantity of reinforce- 
ment required to produce a significant increase in the bearing 
capacity is high, and since the limit state of the reinforcement - 
soil bond is reached at an early stage, the reinforcement -soi 1 
interface slip (rather than the stiffness of the reinforcement) is 
the governing factor. The analysis brings out that the provision 
of the reinforcement helps in spreading the load, thereby causing 
mobilisation of soil resistance over wider area and at shallow 
depth. 

Andrawes et. al. (1982) have presented a finite element 
analysis to predict the load -settlement behaviour of soil- 
geotextile systems. The variational approach is used to obtain 
the stiffness matrix for the soil elements. The geotextile is 
represented by linear elements which have no bending stiffness. 
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The soi 1 —geotext i le interaction is assumed to be only frictional 
and simulated by spring elements of zero length connecting the 
nodes of soil elements to those of geotextile element. Results 
obtained from this analysis are compared with those obtained from 
laboratory model tests which consisted of a geotextile reinforced 
sand layer loaded to failure by a strip footing. The measured and 
predicted data for the load-sett lement curve show that for the 
case where the geotextile was placed at a depth equal to half the 
width of the footing (0.5B) the difference in these values are 
within lO/C upto peak footing load. For the case where the 
geotextile is placed at a depth equal to the width of the footing 
(l.OB) the difference is again within lOV. but only upto 85% of the 
peak load. Beyond this, the predicted behaviour rapidly diverges 

and grossly overestimates the bearing capacity of the footing. 
The model is able to predict the behaviour upto 85% of the peak 

load but fails beyond this load. This is because the model is in- 
appropriate as local failure in the soil occurs which the authors 
say could not be accommodated in the finite element procedure. 

In order to verify the results of their model tests. Love 
et. al. C1987) formulated a finite element program. This program 
was able to handle large diplacements and strains induced in the 
actual model. The subgrade was modelled as an elastic-perfectly 
plastic material with the maximum shear stress equal to C . The 

U 

fill material was modelled as an elastic frictional material 
obeying Matsuoka yield criterion. The reinforcement was treated 


as perfectly rough, so 

that 

any failure 

would 

occur 

in soil 

elements adjacent to 

the 

reinforcement 

rather 

than 

at the 
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interface. Yielding of the reinforcement was not consideredr as 
it was not observed in the model test, and compressive stress was 
not allowed in the reinforcement. The results of this analysis 
compare well with those obtained from model tests for reinforced 
and unreinforced cases. The authors have commented that this 
model may be used with some confidence to perform accurate 
predictions for full-scale structures. They have assumed that no 
compressive stresses are developed in the reinforcement but such a 
condition is possible for reinforcements placed very close to the 
footing. Further, as pointed out by Binquet and Lee (1975a) 
strips at shallow depths fail due to breakage which the authors 
could have included if this model is to be used for predictions 
for full-scale structures. 

Burd and Houlsby (1989) have presented a finite element 
model capable of computing the behaviour of a reinforced unpaved 
road deforming in plane strain. The calculations are based on the 
use of a Von Mises plasticity model for the clay subgrade and a 
model based on the Matsuoka yield function for the granular fill. 
The reinforcement is modelled using elastic membrane elements. 
The formulation is based on a large strain, large displacement 
approach. The finite element model is used to back-analyse 
laboratory scale model tests in which a footing is Jacked 
monotonical ly into a model reinforced unpaved road. There is good 
agreement in the load -displacement curve for reinforced and 
unreinforced systems between the model test results and finite 
element results. The authors say that the finite element model 
may be used to obtain values for parameters that are necessary for 
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a -full understanding of the behaviour of reinforced unpaved road 
but are difficult to measure in physical tests. The authors give 
plots of two such parameters, viz. strain in the reinforcement and 
normal stresses acting on the reinforcement. They conclude that 
the finite element model is useful in estimating values of those 
parameters that are difficult to measure through model tests. 

Burd and Brock lehurst (1990) have reported a small 
displacement finite element parametric study on unpaved roads 
subjected to single monotonic loading. The load is applied by a 
rigid, rough footing of width, 2B, placed on the surface of the 
fill. An elastic-perfectly plastic model is used to model the 
clay subgrade and an elastic-perfectly frictional model is used 
for the fill. The reinforcement is modelled as a perfectly elastic 
material with an infinite yield stress in tension and zero yield 
stress in compression. It is also assumed that no slip occurs at 
the soi 1 -reinforcement interface. It is seen that the variations 
in reinforcement stiffness have a significant effect on the 
magnitude of shear stresses on the upper and lower surfaces of the 
reinforcement and on the tension developed in the reinforcement. 
The authors conclude that for static loading excessively stiff 
reinforcements are not very beneficial since large shear stresses 
are developed, which they say, do not contribute to an improvement 
in static bearing capacity of the road. Further, these large 
stresses induce large values of reinforcement forces. There appear 
to be some reservations with this model since it considers only 
small displacements. Reinforcements may be subjected to compress- 
ion at shallow depths. Further, higher shear stresses help 
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of a line element with no bending stiffness. They carried out a 
series of analysis in which the reinforcement characteristics were 
linear and the elastic modulus, J (defined as the force/unit 
width/unit strain), varied widely (J=0, 150, 1500, 15000KN/m) . 
Their results show that the reinforcement helps in reducing 
deformations if it is sufficiently stiff and strong. Shear 
stresses are mobilised along the interface which eventually induce 
tension in the reinforcement. The effect of the reinforcemnt is 
limited either by its yield or by slip at the interface. 

2.3 Laboratory Model Tests 

A large number of model tests are available in literature 
which aim at studying the behaviour of reinforced soil and 
understanding the mechanisms involved. Since the material avail- 
able is exhaustive only a few test results have been cited here. 

Binquet and Lee (1975b) investigated the mechanisms and 
potential benefits of using reinforced earth slabs to improve the 
bearing capacity of granular soil. They conducted model tests 
with strip footings on reinforced foundations for three 
conditions; (1) Homogeneous deep sand; (2) sand above an extensive 
layer of very soft material; and (3) sand above a finite size 
pocket of very soft material. The reinforcing material consisted 
of 13mm wide strips of household aluminium foil. Tests were 
conducted for a Linear Density Ratio (LDR) defined as the total 
width of the strips per unit width of the footing of 0.425. Three 
series of tests were conducted and results of these tests showed 
that the load-settlement and the ultimate bearing capacities of 
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the footing can be improved by a factor of about 2-4 times above 
the same load-settlement or bearing capacity of an unreinforced 
soil under identical conditions. Results were obtained to study 
the effects of the depth of the top layer of reinforcement and 
number of layers of reinforcement. It was observed that the 
bearing capacity continued to improve with increasing number of 
layers upto six to eight, beyond which there was little additional 
improvement. It was also seen that the optimum location of the 
top strip was at a depth of 1.33B, B being the width of the 
footing and with further increase in this depth the BCR value 
decreased. 

Akinmusuru and Akinbolade (1981) have presented results of 
laboratory model tests with square footings on a deep homogeneous 
sand bed reinforced with flat strips of the rope fibre material. 
The parameters they studied were the horizontal spacing of the 
strips, vertical spacings between horizontal layers of reinforce- 
ment, number of layers, and depth below the footing of the first 
layer, on the bearing capacity of the footing. The fibre strips 
were 10mm wide and thick. The strips of length Im were placed in 
both directions one over the other. The loads were applied in 
small increments and the load -displacement curve to failure for 
each experiment was obtained. The bearing capacity ratio (BCR) was 
defined with the difference that it was restricted to only the 
ultimate failure condition. The study of the effects of the 
various parameters on the BCR revealed the following conclusions. 
1. The BCR decreased with increase in the horizontal and the 
vertical spacings of the strips. 



2. There were three types of failure dependinQ on the vertical 
spacing of the strips which were in agreement with the modes of 
failure identified by Binquet and Lee. 

3. The maximum BCR was observed for the depth of the first layer 
equal to 0.5B, B being the width of the footing. 

4. For number of layers, N, increasing upto three, there was a 
sharp increase in bearing capacity but with further increase in 
N no appreciable improvement was observed. 

The results from the tests conducted by these researchers are in 
agreement with results reported in literature. 

Fragaszy and Lawton (1984) conducted a series of laboratory 
model studies to determine the influence of soil density and 
reinforcing strip length on the load-settlement behaviour of 
reinforced sand. The strips were placed in such a way that the 
linear density ratio (LDR) defined by Binquet and Lee (1975) was 
477. which was slightly higher than the values of the latter. 
Tests were conducted on unreinforced and reinforced sands and the 
effects of strip length and soil density on the bearing capacity 
ratio (BCR) were studied. Their results showed that when BCR was 
calculated at a settlement equal to 1051 of the footing width, the 
BCR was independent of soil density. When calculated at a 
settlement of 4X of the footing width, the percentage increase in 
bearing capacity appeared to be less for loose sands than for 
dense sands. Failure of rectangular footings on dense reinforced 
sand occured at a larger settlement than for an identical footing 
on unreinforced sand at the same density. As the strip length 
increased from three to seven times the footing width, the BCR 
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incrGased rapidly. Additional strip length did not appear to 
affect significantly the bearing capacity. 

The authors compared their results using Binquet and Lee's 
(1975a) theoretical analysis of reinforced earth slabs. Their 
results agreed well with the prediction of Binquet and Lee (1975a) 
that the bottom layer of reinforcement breaks first. The authors 
comment that the equations developed by Binquet and Lee are very 
sensitive to the value of friction angle at failure. 

Guido et.al. (1985) conducted laboratory model tests to 
study the bearing capacity of shallow foundations reinforced with 
geotextiles. Their experimental set-up included square sheets of 
geotextile placed concentrically under a square footing. Their 
results show that at smaller deformations, the full benefits of 
the presence of the geotextile are not exhibited. They comment 
that for the settlement to footing width ratio greater than 
0.017, the fabric deforms sufficiently to mobilise its tensile 
stress, thereby increasing the load carrying capacity. For a given 
value of the depth of the top reinforcement, the BCR increses with 
decreasing values of the vertical spacing of the reinforcements 
and increases steadily with increase in number of layers from 0 to 
3. There is little improvement in the BCR with further increase 
in the number of layers of reinforcement. Guido et. al. (1986) 
have made a comparative study on the behaviour of the reinforced 
earth slabs using geogrids and geotextiles as reinforcements 
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Love et. al. (1987) conducted small-scale model tests to 
study the effectiveness of geogrid reinforcement, placed at the 
base of a layer of granualar fill on the surface of soft clay. To 
make the modelling as realistic as possible, the various 
components were scaled by a factor four. The reinforcement used 
was a version of Tensar SS geogrid. In the tests, monotonic 
loading was applied by a rigid footing, under plane strain 
conditions, to the surface of reinforced and unreinforced systems 
using a range of fill thicknesses and subgrade strengths. 
Deformations of the subgrade and of the geogrid were measured from 
photographs. One set of tests was performed using a double 
footing. The results of the model studies indicated that the 
geogrid reinforcement tends to reduce the shear stress transmitted 
to the surface of the clay subgrade. The authors also comment 
that the failure mechanisms in the clay are mobilised at small 
deformations of the fill, and therefore large deformations are not 
necessary for the benefits of reinforcement to be felt. At large 
deformations, where these are permissible, additional benefit is 
obtained from the membrane action of the reinforcement. 

Setty and Rao (1987) conducted model tests by reinforcing 
locally available lateritic soil with synthetic fibres. The 
reinforced soil contained discontinuous, discrete, and randomly 
oriented fibres with an aspect ratio of about 25. It was observed 
that there was some reduction in the angle of internal friction 
and a substantial increase in the value of cohesion due to the 
presence of fibre reinforcement. There was an overall improvement 
in the CBR values as well as the tensile strength of the soils. 



26 

The expansion ratio of the soil sample decreased. The authors 
infer that fibre reinforced soil can be effectively used in the 
lower layers of pavement to improve its tensile resistance and in 
the subgrades of rural roads where surface cracking becomes a big 
hurdle. 


Miura et.al. (1988) conducted model and field tests of 
reinforced pavement by polymer grids on very soft ground. Field 
tests were performed on a 300m long test road with six sections of 
different types of pavement. Cyclic load tests were carried out 
on reinforced and unreinforced model test pavements. The model 
and field test results showed that polymer grid placed at the 
interface of the subbase and the subgrade suppresses the 
settlement under repeated loading. From the field data obtained 
after one and three months of construction it was observed that 
the modulus of subgrade reaction was smaller than that of 
conventional pavement. It is reasoned that this might be partly 
due to the overestimation of the grid function and partly from 
insufficient compaction of base owing to the action of the polymer 
grid. The results from a finite element analysis indicated that 
the polymer grid is not effective in suppressing the surface 
settlement of the pavement. The authors suggest that interlocking 
effect in polymer grids is an important function which is not 
accounted for in the analysis and that while using polymer grids 
as reinforcing material for pavements, care should be taken while 
compacting so that no gap is left between the grids and the 


underneath layers. 



Dembicki and Alenowicz (198B) conducted model tests on 
two~layer subsoil consisting of a layer of sand fill underlain by 
soft soil. They have investigated the influence of a geotextile 
inclusion on bearing capacity of the system. In order to simulate 
a truck axle under plane strain conditions, the surface of the 
two-layer system was subjected to static loading applied by two 
steel plates. The model tests also aimed at observing the failure 
mechanism of reinforced and unreinforced systems. The geotextile 
used was needlcpunched nonwoven. From the load -sett lement curves 
for the unreinforced and reinforced systems it was observed that 
at small plate penetration values there was little difference 
between both the systems. On the other hand, at large 
displacements of the plates, influence of the geotextile on 
measured pressure was significant. They have compared the 
indicated load-carrying capacity due to inclusion of a geotextile 
with the reduction in pressure on the weak subsoil due to the 
fabric calculated from the analysis given by Giroud and Noiray 
(19811. The results show that the calculated reduction of pressure 
differs from the increase in the load-carrying capacity indicated 
in the tests. The authors attribute this difference to the fact 
that tests consider the changes in failure mechanisms as well as 
tensile force induced in the geotextile, while the calculated 
values show reduction connected only with the latter effect. 

Huang and Tatsuoka (1988) performed a scries of model tests 
to investigate the reinforcing mechanisms for strip reinforced 
foundation beds. Strips made of phosphor bronze were placed 
beneath a strip footing. The tests were aimed at studying the 
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reinforcing effect of short strips having the same length as the 
footing, the effects of length of reinforcement, number of 
reinforcing layers and the density of reinforcing strips. Tests 
were also conducted on deep footings. The researchers have given 
figures for the strain fields for the cases of unreinforced and 
reinforced grounds. It was seen that the ground reinforced wjth 
short strips behaved very similar to the deep footing with depth 
of the footing equal to the total depth of reinforcing layers 

. This trend indicated that the failure was progressive and 
the compressive strength of the active zone beneath the footing 
controls the bearing capacity of the ground. The authors call 
such an effect as ^the deep footing effect". Based on the test 
results they suggest two modes of failure. The first mode is 
observed for densely reinforced conditions wherein the shear bands 
starting from the edges of the footing extend downward to the 
depth, D , then form a wedge beneath the reinforced zone. The 

r 

second mode is observed for lightly reinforced conditions wherein 
the shear bands start from the footing edges and form a wedge 
within the reinforced zone. In this case, the bearing capacity of 
the reinforced ground is controlled by the strength of the block 
immediately beneath the footing. 

Model tests were conducted by various researchers 
(Sreekantiah, 1988; Shankar iah and Narahari, 1988; and Sargunan 
and Hussain, 1988) on square and strip footings on sand beds 
reinforced with different types of reinforcements like aluminium 
foil, 6. I., and bamboo strips. Their experimental investigations 
revealed that the bearing capacity of the foundation increases due 
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to the introduction of reinforcement. There is no appreciable 
increase in the BCR for depths of placement of reinforcements 
beyond 1.4B, B being the width of the footing. With the increase 
in the number of layers of reinforcement, the BCR increases. 

To study the effect of a geotextile placed at the interface 

of a granular fill and a weak clay layer. Das (1989) performed a 

number of model tests on strip and square foundations. He has 

studied the effects of the height of the granular fill to width of 

the footing ratio (H/B), the depth of embedment of footing CD ) 

f 

and the width of geotextile on the ultimate bearing capacity 

of the foundation. The length of the geotextile was equal to the 
length of the foundation. The test results indicate that the 
ultimate bearing capacity increases to a maximum at H/B=0.75 and 
decreases to a constant value at H/B=2 for D^/B=0.5. The author 
has indicated that the ultimate bearing capacity for a given H/B 
ratio can be expressed as 

q = q + Aq + Aq C2-3) 

or 

BCR = 1 + ACBCR) + ACBCR) C2.4) 

8 O 

where q =ultimate bearing capacity with geotextile interface; 

u(HXB> 

Aq =net increase of ultimate bearing capacity due to the sand 

\Jk\aty 

layer of height, H, only; Aq =net increase of the ultimate 

u<a> 

bearing capacity due to the inclusion of geotextile; and BCR is 

the bearing capacity ratio. Using the test results the variations 

of ACBCR) and ACBCR) were computed from the equations which 
s a 

indicated that the beneficial effect of the geotextile reaches a 


maximum value at H/B=0.75. The geotextile has practically no 
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effect for H/B>1.5 for strip foundations and for H/B>1 for square 
foundations. The computed results also indicated that without the 
inclusion of the geotextile, the BCR increases upto a maximum of 
about 2.8 at H/B=l.£>, while the BCR increases to about 3 at 
H/B=0.75r for a strip foundation if a geotextile is placed at the 
sand-clay interface. The corresponding values for a square 
foundation are 1.75 at H/B=l without a geotextile and 1.9 at 
H/B=0.5 with a geotextile at the interface. Also, the most 
economical width of the gcotextile for a strip and a square 
footing are 4B and 3B respectively 

Shimizu and Inui (1990> carried out model tests on soft soil 
reinforced with geotextile cell wall. A six-sided cell of 
geotextile wall frame was buried in the subsurface of the soft 
ground and vertical loads were applied at the centre of the cell. 
The loading tests were carried out to find the necessary 
conditions for realizing the increase in the bearing capacity. 
Furthermore, in order to investigate why and how the bearing 
capacity increases. X-ray photographs were taken of the 
displacements of small lead pieces embedded in the model. The 
load -set t lement behaviour showed that installation of the cell 
increases the bearing capacity of the ground. The critical value 
of the height to width ratio of the cell, is about 1.5 for which 
the bearing capacity is maximum. The stiffness of the geotextile 
affects the bearing capacity too. Furthermore, it was observed 
that if the horizontal tensile strength is not high enough to 
resist the expansion caused by the ground deformation, the cell 
would be ruptured. The X-ray photographs showed that at 
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relatively small settlement, the movements of the particles is 
restricted by the cell wall. 

In order to study the behaviour of reinforced soil beds with 
reinforcements within the size of the footing, Sridharan et. al. 
(1991) conducted model plate load tests- They used two types of 
reinforcements, viz., mild steel grids and aluminium strips and 
conducted tests with varying number of layers. To compare the 

behaviour of this foundation system they conducted tests on 
unreinforced sand bed with embedded footings with the depth of 
embedment equal to the depth to the bottom most layer of 

reinforcement. The test results indicated that for stiff and 

sufficiently closely spaced CS <B/4) reinforcements, the increase 

V 

in load bearing capacity was equal to that of embedding the 
footing to a depth equal to that of the bottom most layer of 

reinforcement with the maximum depth of embedment being equal to 
B. B is the width of the footing and S , the vertical spacing 

V 

between the strips. These results are similar to the results 

obtained by Huang and Tatsuoka (19881. They have given an 

expression for the bearing capacity ratio (BCR) by modifying the 
expression given by Binquet and Lee (1975) to account for the 
embedment effect as 

BCR = (P + kAP + P )/P (2.5) 

S S r S 

where AP is the increase in load carrying capacity due to the 
s 

embeddment of the footing and k=(EA) . /(EA) with a 

maierval 

minimum value of (EA) , equal to 60,000 kg per cm width of the 

•t— L 



footing. and are the load carrying capacities of the soil 
and the reinforcement respectively. 

2.4 Field Tests 

Jones and Dawson (1990) have studied the performance of 
reinforced granular foundations on soft clay in situ to exploit 
this system for foundations for structures. Footings were 
constructed in a 3.5mx2.5mxl.5m trenches. The reinforcement 
consisted of welded mesh with PVC coated galvanised high tensile 
wire. The footings were loaded and their performance monitered 
with a range of instruments. From the load -set tlement curves it 
was evident that the total settlement was within tolerable limits 
even under high loads. The authors comment that reinforced 
granular soil is a promising technique for lightly loaded 
foundat ions. 

Lindh and Ericksson (1990) performed field tests on two 
stretches of pavement by incorporating a 100mm layer of sand in 
the pavement mixed with short plastic fibres. The performance of 
the pavement was studied for over a year and it was observed that 
no rutting took place. It was also seen that there was no 
difference in the stiffness of the layer of stretches with and 
without reinforcement. The authors comment that the performance 
will be observed for further three years. They conclude that sand 
reinforced with short fibres is an interesting method to create 
cohesion in non-cohesive soils. This method also brings about an 
improvement in the strength of the soil. 


The state-of-art of reinforced soil indicates that lot of 
work has been in progress to understand the mechanisms of 
reinforced soil and to predict the increase in the bearing 
capacitiy of reinforced soil. Inumerable model tests have been 
performed using various types of reinforcing materials to improve 
the performance of footings on soft soil. Attempts have been made 
to put forth theories to predict the behaviour of soils reinforced 
with strip and sheet form of reinforcements. Few theories are 
available, most of which either deal with the limit equilibrium 
techniques or are based on the Finite Element approaches. 
Further, most of them aim at predicting the increase in the 
bearing capacity of reinforced soil. Many a times and especially 
for pavements, settlement governs the design criteria. It hence 
becomes very essential that the reduction in surface settlements 
due to reinforcement is known. From the review of the literature 
it is evident that not much effort has been made to devise 
theories to predict the reduction in surface settlements due to 
reinforcements. The present study aims at predicting the 
settlement reduction of the surface due to strip and sheet form 
(axi -symmetric) of reinforcements placed beneath rectangular and 
circular (for the axi -symmetr ic case only) loaded areas. The 
elastic continuum theory is resorted to. The tension-membrane 
action of the reinforcements has however, not been considred. 



CHAPTER 3 

SOIL-RIGID STRIP REINFORCEMENT INTERACTION - SHEAR STRESSES 


3.1 Introduction 

As mentioned in chapter 1, the surface loading causes the 
points to move downward and outward. The points along the 
reinforcing strip placed at depth tend to move vertically and 
horizontally, away from the centre. The strip, in this chapter is 
considered to be rigid Cincxtensible in its axial direction!. As 
a result, the horizontal movements of the points along the strip 
are prevented by the shear stresses mobilised at the interface. 
The displacements of the points along the strip, due to the 
surface loading are calculated using the Boussinesq^'s equation for 
horizontal displacements of points within an elastic continuum due 
to a point force on the surface. The horizontal displacements of 
points along the strip-soil interface due to the mobilised shear 
stresses are computed using Hindi in"'s equation for horizontal 
displacements of points within an elastic continuum due to a 
horizontal force within the continuum. The con^atibility of 
displacements at the nodes is satisfied by equating the net 
displacements to zero. This results in a set of simultaneous 
equations, the solution of v^ich yields the desired shear stresses 
mobilised at the interface. These stresses, by restricting lateral 
displacements tend to push the soil above, resulting in the heave 
of the surface which, in effect is the settlement reduction of the 


points on the surface- 
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3.2 Single Strip 

3.2.1 Problem Definition 

The problem consists of a rectangular area of size, 2Lj.x2B^, 
transmitting a uniform load of intensity, q, on the surface of a 
semi -infinite continuum (Fig. 3.1). A rigid reinforcing strip of 
sizEj 2L x2B , is placed centrally within the soil and below the 

r r 

loaded area at a depth, U - The width of the strip is relatively 

o 

small and its thickness, t , negligible. 

r 

Since the width of the strip is assumed to be relatively 
small, the shear stress developed are assumed to be constant over 
the width and to vary only along its length. In this chapter only 
shear interactions are considered. Hence it is assumed that the 
displacements of the strip in the vertical direction (i.e, with 
respect to the normal stresses) have no effect on the shear 
stresses mobilised along the interface. The elastic continuum 
approach is used to obtain these shear stresses by carrying out an 
interaction analysis. These stresses, acting inward on the soil, 
cause the soil above the reinforcement to move upward. This 
results in the upward displacement of the points along the 
surface, the origin showing the maxmium displacement. Hence, the 
net settlement of the surface is reduced due to reinforcing action 
of the strip. 

3.2.2 Analysis 

Since the analysis is based on the elastic continuum 
approach, the soil is assumed to be homogeneous, isotropic and 



linearly elastic- Firstly, the horizontal displacements along the 
soi 1 ^re inf or cemen t interface due to the applied load are 
calculated- This is done by using Boussinesq*^s solution for 
horizontal displacements at any point within an elastic half -space 
due to a point load acting on the surface of the elastic 
half -space. The Boussinesq^s problem is as shown in Fig- 3-2- 
The displacement of point as shown in Fig- 3-2, along the 
radial direction due to load, P, is given by 


P(l+i> ) rz )r 

s j _ 

2nE R ^ (R + z) 

s 


(3.1) 


where r and z are the radial distance and depth of point A, with 
respect to the load P, f<=-/ r* + z* , E and v are the modulus of 

S S 

deformation and Poisson^'s ratio of the medium. 


Eq. 3. 1 is integrated over the loaded area to obtain the 
horizontal displacement at point A, (Fig. 3.2), within the 
half-space due to the uniform load of intensity, q, with 
appropriate modifications for the direction. That is, the radial 
displacement, , at A due to load intensity, q, on an elemental 
area AA is 

qAA(l+v> ) rz (1-2 p )r 

Ap = 5 — r — 

2TrE R ^ R* (R+z) 

s 

The lateral displacement, Ap^, along the x -direct ion is 

Ap = Ap Cosoi 


(3.2) 


(3.3) 





where ot is the angle between the radial line -from the elemental 
area to point A and the x-axis. The total lateral displacement, 
p , of point A is 

^ J ~ 1 <3-41 

a"^ a 

The horizontal displacement, p^, along the reinforcing strip due 
to the uniform loading, q, is then evaluated as 




pf = r r rjlL - 

J J 2nE R ^ (R+z) ^ 


Ir. 
s 
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<3.51 
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In practice, it is found convenient to carry out the 
integration of Eq- 3-5- numerically. For this purpose the width 
of the loaded area is sub-divided into n sub-elements and length 

B 

into n sub-elements- Due to symmetry only half the length of the 

Li 

strip is considered. This half-length is divided into N elements 
and over each element the shear stress is assumed to be constant. 
The horizontal displacement, P . .t of the centre of the i^*^ 

XVJ 

element (Fig. 3.31, due to loading, q, on an elemental area, AA, 
of element j, on the surface is 
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where r 




and 


R=-/~r‘ 




U being 
o 


the depth to the 


reinforcing strip 



The horizontal displacement, the node i, due to the 
rectangular surface load is obtained by integrating Eq. 3-6 over 
the entire rectangular area numerically, as 


I 


q(l+v ) rU il-2v )r. 
s r o s 


^ 2nE R 

’•f = ‘ = ‘ 


[ I 

— ^ - 


Coso dA 


i3.7) 


(R+U 


Non-dimensionalising all length parameters with half-width of the 
loaded area, Eq. 3.7 is expressed as 



(3.8) 


where is a dimensionless influence coefficient that depends on 

X 

the aspect ratio of the loaded area, the depth of placement 

of the strip, U /B,, the Poisson ''s ratio, w , of the soil and the 
location of the node i. The horizontal displacements of all the 
nodes along the half-length of the reinforcing strip are obtained 
by evaluating Eq. 3.8 for all the nodes, 1 to N. The vector of the 
displacements for all the nodes is expressed as 

{ ^ { I' } , 

where the vectors \<] and are of size N. 


The horizontal displacements due to the surface loading are 
counteracted by the shear stresses mobilised along the 


soil -reinforcement interface. In order to obtain these shear 



stresses, the horizontal displacements due to these stresses are 


to be computed. It is achieved by using Mindlin's solution for 
horizontal displacements due to a horizontal force acting beneath 
the surface of a semi -inf inite medium. 

The horizontal displacement of a point, A(x,y,z>, due to a 
horizontal force, Q (Fig. 3.4), acting at a depth, c, within a 
semi -inf inite elastic continuum is 


Q (1+v )x 
s 

BnE (1-v ) 
s s 



(3. 10) 


where F 


(3-4v> ) 
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4(1 -i> )(l-2v ) 
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X 


R*(R +Z+C) 
2 2 


] 


R = / x*+ (z-c)* 

1 

R = -/x*+ y^+ (z+c)* 

2 

and X, y, z are the co-ordinates of point A (Fig. 3.4). The 
lateral displacement of any point along the strip due to shear 
stress, T, on an elemental area, dA, is 




TdA(l+v )x 

s 

BnE (1-v ) 

SI SI 



(3.11) 


In order to obtain the lateral displacement, p^, of points 
along the reinforcement due to the mobilised shear stresses. 
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reinforcing strip. Only half the length of the strip is considered 
due to symmetry and this length is discretized into N elements. 
The lateral displacement, of node i due to shear stress, t., 

x>-J j 

acting on clement j (Fig. 3.5) is 





Cl+v ) X 
s 

6rrE (1-v ) 
s s 



T. dA 
J 


(3. 12) 


where x and y are the horizontal distances of node i with respect 
to the elemental area, dA, of the element. The other 
parameters are as mentioned in Eq. 3.10 with c = = z- 


Eq. 3.12 is also integrated numerically to obtain the 
influence of shear stress, t., on the horizontal displacement of 

J 

node i. Non -dimensional ising all length parameters with Eq. 

3.12 is written as 


r 

1 


XL j 




C3. 13) 


r 

where I ^ is a dimensionless coefficient that depends on the 
^ j 

length, L /B , and the depth of placement of the strip, U /B , the 

r f u I 

Poisson'^s ratio, v , of the soil and the location of the node i 

s 

and element, j. For every element j, there exists an image j^, 

(Fig. 3.5), shear stress on which is of the same magnitude but 
acting in a direction opposite to the one acting on element, J. 
The influence of the stresses acting on element j' f 
displacement of node i, is written as 


on the lateral 



2 f - 2 

P Z= I -jr ^ 


(3- 14) 


where is the influence coefficient for the influence on the 

displacement at node, i due to the shear stress on element j' . 

jp 

2 

^ij obtained from Eq. 3-12 wherein x and y are the distances of 
node ir with respect to the elemental area, dA, of the 
element- It is to be noted that the horizontal displacement, 

r 

2 

p . T will be in the opposite direction with respect to the 

Kt J 

r 

displacement, p .. 

XIJ 


Combining Eqs. 3.13 and 3.14, i.e., combining the influences 
of the shear stresses on the element and of its mirror image, 
the lateral displacement of node, i is 



(3. 15) 


To obtain the displacement at node, i due to the shear stresses on 
all the elements along the full length of the strip. Eg. 3.15 is 
integrated numerically over the N elements along the half-length 
of the strip as 



t3. 16) 


The horizontal displacements of the centre of all the elements 
along the half length of the strip is 



the vector of horizontal 


Expressing Eq. 
displacements. 


3.17 in 


matrix form, 
is written as 




s 


[■■]{') 


(3. 18) 


where vectors 




are of size 


N, [ I' ■ 


IS a square 


matrix of size N. An interesting property of this matrix is that 


'•j j*- 


At this juncture, the lateral displacements of the centres 
of all the elements along the strip due to the applied surface 
loading CEq. 3.9) and those due to the mobilised shear stresses 
(Eq. 3.18) are known. The net soil displacement of node i, is the 
difference of the displacements due to the surface load and the 
interaction stresses, i.e.. 


net 




(3. 19) 


The net soil displacements should be compatible with the 
displacements of the reinforcing strip. The solution of the shear 
stresses is obtained by satisfying the compatibility of 
displacements at the soi 1 -reinforcement interface. For a rigid 
( inextensiblc) strip the net lateral displacements arc zero, i.e. 





(3.20) 


Combining Egs. 3.9, 3.18 and 3.20 one obtains 






(3.21) 


or 


_ r ] { T/, } = { l' } 


(3.22) 


Eq. 3.22 gives N equations for r/q, the normalised shear stresses. 
It is solved by the Gauss elimination technique for the shear 
stresses on all elemental areas. 

These shear stresses are developed at the interface of the 
soil and strip. The thickness of the strip being negligibly 
small, it can be assumed that the stresses above and below the 
strip are equal to half the calculated values. The stresses 
mobilised in the soil in turn develop tension in the reinforcing 
strip. The tension is evaluated as under (Fig. 3.6) 

AT. = T AA = T. 2B dl (3.23) 

V V t r 

where AT. is the elemental tension on the element i, and dl is the 

X 

length of element i. The total tension at the centre of element i 
is 

N 

T = i T AA + y T. AA (3.24) 

i 2 i A J 

N 

= 4 T. AA + y T. AA 


T 


(3.24) 
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The tension at the centre of the strip, i.c., at x=0 is 

' N 

'^x=o = 5 = 2 ’■j (3.25) 

j=i 

where AA is the area of element j. The tension at x=L i.e., at 

J r 

the edge of the strip is zero. 

The primary aim of this analysis is to estimate the 
reduction in settlement of the surface due to the presence of a 
reinforcing strip inclusion within the soil. As mentioned earlier, 
the shear stresses mobilised help in causing an upward 
displacement of the soil at the surface as depicted in Fig. 3.7. 
The shear stresses cause the soil above it to move upward and 
inward. As a result, the centre or the origin, 'O', experiences a 
maximum upward displacement. The surface displacement profile is 
as shown in Fig. 3.7. Once, this profile is evaluated the problem 
is solved completely. 


In order to arrive at this profile, Mindlin‘'s solution for 
vertical displacement at any point within an elastic half-space 
due to a horizontal force within the half -space is used. 
Referring to Fig. 3.4, the vertical displacement of point A due to 
a horizontal load, Q, is given by 


Q (l+v> )x , z-c 


8rrE 


s''' f 
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(3.26) 
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Fig. 3.6 Stresses on the Strip 





Fig. 3.7 Settlement Reduction of Surface 


Non -dimensional ising 


the shear -nodulus of soil (G =E /2(l+v )) 

s s s’ 

length parameters with Eg. S.29 is rewritten as 


r B, r 

1 f , 4 

P , = cT- I , T . 

ikj G Skj j 


(3.30) 


where is the vertical displacement of point k, due to shear 
stress on element j. Displacement at point k, due to shear stress 
on the image element i' , on the left half of the strip is 


Xkj 




I , . T. 

Skj j 


(3.31) 


r r 

1 2 

where are dimensionless displacement coefficients 

for the influence of shear stresses on the elements j and j' on 

the displacements of node k, and are functions of the ratios, 

L /B^r O /B,, V and the location of the respective elements, 
r f o f e 

Combining Eq. 3.30 and 3.31, the vertical displacement at node k, 
on the surface is 


*kj 


B ^ r r ^ 

" 'skJ J 


(3.32) 


The displacement of node, k due to shear stresses on all the 
elements along the strip is 




(3.33) 


J = i 


Expressing Eq. 3.33 in matrix form as 
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(3.34) 


The width of the lo-^ded ^res is divided into N points where p is 

( z 


evaluated. The vector 
size N and the matri 




is of size N^, vector 




of 


is of size N.xN. 


Eq. 3.34 is rewritten as 


{ 1 g' { } 


(3.35) 


The elements of this vector are termed as Settlement Reduction 
Coefficients (SRC). SRC at any node k, on the surface is 


P , G 

SRC, = I , = — 

k Sk q 


(3.36) 


where p , is the vertical upward displacement of node k, due to 

3K k 

the mobilised shear stresses. The SRC at the origin, 0, or the 
centre of the loaded area is defined as I 

Sc 

The vector of the SRC along the surface is 

where J is obtained as 
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(3.38) 


where the matrix. 



and vector. 



known. 


3.3 Two Strips at Depth, U ' 

3.3.1 Problem Definition 

In this case, two strips of size, 2L x2B , are placed 

r r 

symmetrically about the x-axis at a depth, U^, within the soil 
mass (Fig. 3.8). Each strip is placed at a distance, S^, from the 
x-axis. The mechanism through which shear stresses are developed 
along the soil -strip interface is the same as that mentioned in 
the section for a single strip. The discretizations of the 
surface loaded area and the strip are the same as in the previous 
case. Since the strips are placed equidistant from the x-axis, 
the lateral displacements along both the strips due to the surface 
loading are equal. Consequently, the shear stresses mobilised at 
the interfaces arc also equal. 

The horizontal displacements along the length of one strip 
due to shear stresses arc not only influenced by the stresses 
along itself but also by the stresses mobilised along the second 
strip. In calculating the horizontal displacement of node i, due 
to stress on element j, (Fig. 3.9), in addition to the stresses 
on element j, and its image j", stresses on elements i" and its 
image j"'r along the second strip arc also to be considered. 
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Inc iden t al 1 y r the stresses on all these elements are the same due 
to symmetry » 


3.3.2 Analysis 

Due to symmetry^ the horizontal displacements only alon 9 
half the strip length due to applied loading are calculated. The 
displacements are computed from Eq. 3.7 with the modification, 
X + CS^-y) . Non-dimensional ising all length parameters with 
the half-width of the loaded area, Eg. 3.7 is written as 



(3.39) 


where is a dimensionless influence coefficient that depends on 

the ratios, L,/B,, U /B,, S /B,, v and the location of node i. 
f f of y f s 

f> . is the horizontal displacement of node i due to the surface 

XV 

loading. The vector of the displacements of all nodes is expressed 
as in Eq. 3.9 as 





q 


where vectors 


{ p' } and I l' } 


are of size N. 


(3.40) 


The horizontal displacements due to shear stresses along the 
strip are calculated using Hindi in'^s solution as described in 
Eq. 3.10. The displacements of node i due to shear stresses, 
acting on elements j, J'r J''' a^id as explained in Fig. 3.9 
are obtained from Eq. 3.12 and arc expressed in a dimensionless 


form as 
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where 

I 

'■j 

I ^ 

i.“, 

'■J 

and 

>• J 


coefficients for the influence of 


j» j' T S" and j"' respectively. 


U /B,, S /B,, V and the location 
Of y f s 


D4 

C3.41) 

(3.42) 

(3.43) 

(3.44) 

are displacement influence 
the shear stresses on elements 
and are functions of L /B,, 

r f 

if the respective elements. 


Combining equations 3.41 through 3.44, the lateral displacement of 
node i, due to shear stresses on all N elements is 



j = * 



(3.45) 


where the summation is carried out for stresses acting on all N 
elements. The vector of horizontal displacements is 



[■■]{•} 


(3.46) 


where the vectors^ ^ ^ size ^ j| ^ 

square matrix of size N- The difference in the matrix, I J, in 
this case as compared to a single strip is that the matrix here is 
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a function of the distance, of the strips from the x-axis 

also, while it was not in the former case. 

The net lateral displacement for a ri9id strip is zero. 
Satisfying compatibility of displacements as in Eg. 3.20, shear 
stresses are obtained from Eg. 3.22, written again as 

[ ^ ] { ^ } t3.47) 


Eg. 3.47 is solved for the N values of r/g on all the elements. 

Once the inter facial shear stresses are computed, the 
tension along the strips and settlement reduction coefficients for 
points on the surface are obtained based on the formulation 
discussed for the single strip. The tension in the strips is 
arrived at using eguations 3.22 through 3.25. 

In calculating the settlement reduction coefficients, 

appropriate modification for the influence of the two strips, i.e, 

for the distance, S , is to be made. The vertical displacement at 

y 

any point, k, on the surface due to the shear stress on element j, 

on the strip is given by Eg. 3.29 wherein yj^^ now includes the 

parameter, S , while all other parameters remain unaltered, 
y 

The vertical displacements of point k, due to shear stresses on 

elements j, j , j and j-'" 


(Fig. 3.9) are given by 



(3.48) 
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(3.50) 


(3.51) 


where I_* » I_f r I_f.r and I are dimensionless coefficients 

skj Skj Skj Skj 

for the influence of the shear stresses on the elements j, j' , j'" 

and y" respectively and arc functions of L /B~ U /B,, S /B,r 

r f Of y f 

V and the location of the respective elements. Combining 

equations 3.48 through 3.51, the vertical displacement, . of 

zk 

node k, due to shear stresses on all N elements is 



j = i 


(3.52) 


The vector of the vertical displacements of points along the 

surface is 




The SRC at point k on the surface, due to the stresses along the 



C3.54) 
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The vector of the SRC along the surface is expressed as 


{ '. } = -57 { p. } 


(3.551 


where 




is obtained as 




(3.56) 


where the matri 


ix. [ ] 




and vector, i t J- are known. 


3.4 Two Strips - One Below the Other 


3.4.1 Problem Definition 

Two strips of size, 2L x2B , placed beneath the centre of 

r r 

the loaded area along the x-axis, one at a depth, U^, and the 
other at a depth, are considered (Fig. 3-10). The interference 
of one strip on the other is studied. It is assumed that the 
horizontal displacements at depths, and U^, due to the surface 
loading remain unaffected by the presence of the other strip, i.e 
the Boussinesq'^s equation is applicable for both the depths. 

Since the strips are placed at two different depths, the 
horizontal displacements at the two levels are different and 
consequently, the stresses mobilised are also different. The 
stresses are symmetric about the z-axis. The stresses mobilised 



along each strip influence the horizontal displacements of points 
along the length of the other. In computing the horizontal 
displacement at node i (Fig. 3.11), of strip 1, the influences of 
stresses, t‘, on elements, j and its image j", along strip 1, and 
the stresses, t^, on elements, j and its image j'', along strip 2 
are to be considered. Similarly, for horizontal displacements of 
points along strip 2, the influence of stresses along strip 1 and 
strip 2 are to be accounted for. 


3.4.2 Analysis 

Eq. 3.7 is used to calculate the horizontal displacements of 
points at depths, and U^- The vector of horizontal 
displacements of the N nodes along the half-length of strip 1 is 




(3.57) 


where 




is a dimensionless influence coefficient vector of 


size, N, that depends on the parameters, L^/B^r 
location of the nodes. 


The vector of horizontal displacements at N nodes along the 
half-length of strip 2 is 

where ^ ® dimensionless influence coefficient vector of 

, N, that depends on the parameters, 


size 


and 
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Rg. 3.11 
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.1 


Strip 1 


• f ^ » " l 

i j N 

jf 

■ - i Strip 2 

j 


of the Strips 




the location of the nodes. 


Eq. 3.10 is used to compute the horizontal displacements of 
points along the lengths of the strips due to the mobilised shear 
stresses on them. The horizontal displacement of node i of strip 
1 (Fig. 3.11)t due to shear stresses on all elements along strip 1 

is 


li 


X -V 


C3.59> 


j »= i 


and that due to stresses on all elements along strip 2 is 


r I - r r - 


(3.60) 


r r 

where I.^ and I.^ are dimensionless coefficients for the influence 

'-j '-j 

of shear stress, r*, on elements j and its image j' respectively, 

j 

r r 

of strip 1. I ® and I.* are dimensionless coefficients for the 
'■j '■J 


influence of shear stress, t., on elements 


and its image 


respectively, of strip 2. Apart from being functions of 

r 

and the locations of the elements, I. . and I. . are functions of 

i-j '•J 

r r 

the depth ratios U /B and, I.* and I. . are functions of the depth 

of VJ VJ 

ratios, U /B„ and U /B . From Eg. 3.59, the vector of horizontal 
of ^4 f 

displacements, of all N nodes along strip 1 due to shear 

stresses along strip 1 is written as 




(3.61) 



where 


shear 



stresses. 


the 


coefficient matrix 
, along strip 1 on 


for the influence of 
displacements of nodes 


along strip 1 itself. 


The vector of horizontal displacements, 1 ^ “ }’ of N nodes along 
strip 1 due to shear stresses along strip 2, from Eq. 3.60 is 




(3.62) 


where 


IS a 


stresses 


[ ] 


coef f icient 
along strip 


matrix for the influence of shear 
2, on the displacements of the nodes 


along strip 1. 


Similarly for strip 2, the vectors, | | and | p^“ 

horizontal displacements are written as 


}■ 


of 



{ '’1” } 




(3.63) 

and 
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(3.64) 

where 

[i-i 

and 1 

I ** j arc 

coefficient matrices 

for the 

L J 

influence of 

shear 

stresses, 

along strip 1, and j 

along 

strip 2 

respectively, on the 

displacements of the 

N nodes 

along 

strip 2. 

All 

the matrices. 

viz. [ I “ ]. [ I “ ]■ 

[ ''1 


and [ I square matrices of size N, and the displacement 


and stress vectors are of size N 


Combining Eqs. 3.61 through 3.64 the vectors of horizontal 
displacements of N nodes, each along strips 1 and 2 are 


{ ° el [ ^ } 


C3.651 


where 


{ p' } and { T I 


are vectors of size 2N and 


[^1 


IS a 


square matrix of size 2N. 


The vector; 




compr ises 


of vectors, ^ and ^ and 


is written as 


} = I J 

where { p/ } “ { } + { p“ } 

and 


(3.66) 


The vector comprises of vectors I" 


written as 


{ ”" } = L L '"‘J L J = I ^ \ 


(3.67) 


The matri 


X, [ r ] 


is written as 



(3. 68) 
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Combining equations 3.57 and 3.58 the vector of 
displac ementf. of the nodes along strips 1 and 2, is 




where 


{ p' } a.,d 1 1' j. 


are vectors of size 2N. 


Vector 




is exp'ressed as 


(p' f = I Ip'JIp' J J = Lp' p p ■ 
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comprises of 


vectors 


{ } 


For rigid strips, the net lateral displacement is zero 
compatibility of displacements at all nodes, it can 
that 




hor i sontal 

(3.69) 


p J 

(3.70) 


Sat isf ying 
be written 


(3.71) 


(3.72) 



The solution of these 2N equations gives the shear stresses on N 
elements along strips 1 and 2 respectively. 

The tensile forces along the strips and the settlement 
reduction coefficients due to these stresses are calculated on 
similar lines as explained in the earlier sections. 


3.5 Multiple Strips 

The performance of multiple strips are studied here. Two 

pairs of strips of the same dimensions C2L x2B ) are placed below 

r r 

the loaded area, symmetrically about the x-axis at a depth, U^, 
(Fig. 3.12). One pair consists of two contiguous strips each at a 
distance, S =B , from the x-axis. This pair is equivalent to a 

y± T 

single strip of width, 2B , placed symmetrically below the loaded 

r 

area (the width of each strip in the pair is half that of a single 
strip). Each strip of the second pair of strips is placed at a 
distance, S , from the x-axis. The interference of one strip on 

y2 

the other is studied in this case. 


3.5.1 Analysis 

As explained in the previous section, the horizontal 
displacement vector of N nodes along the length, L^, of the first 
pair of strips is 


{'-i-fci-;) 


where | j and | 

non-dimensional vector of 


(3.73) 


are vectors of size N, 
influence coefficients 



that depends 


a 

on 



the parameters, L/B,, aad the loeation of 

node i. Again, due to symmetry only one strip of each pair is 
considered . 


The horizontal displacement vector of the N nodes alonQ 
half~length, L^, for the second pair of strips at a distancer 
from the x-axis, in non-dimensional terms is 


the 


S 

y2 


where 



dimensionless vect 
parameters, 
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^ sre vectors of size N and ^2 ^ 

or of influence coefficients that depends on the 
U /B,, S /B,, V and the location of node i. 

O f yZ f S 


Combining Eqs. 3.73 and 3.74, the vector of horizontal 
displacements for both the pairs of strips is written as 



1 _ 

f " E“ ^ 
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(3.75) 
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|- and f 


are vectors of size 2N. -I p^ 
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^ comprises 

of vectors. 


■ and 

p^^ while ^ comprises 

of vectors 

{ 1 ' } and j 






The horizontal displacement of node i (Fig. 3. 13> along the first 
pair of strips due to the shear stresses on the N elements of the 
first pair can be written as mentioned in Eq. 3.45 as 






(3.76) 
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r^ r^ r 

where influence coefficients for 

the displacement of node i due to shear stress, t‘, on elements j, 
j'r J and j ^ , the meanings of which are as discussed in 
section 3.3.2. The vector of horizontal displacements of all N 
nodes then is 


{ 






(3.77) 


where 


[ ] 


is a square matrix of size N and 


[ ] 


IS a 


matrix of the influence of shear stresses along the first pair of 
strips on the displacements of N nodes along the first pair of 
str ips. } and pj'*' ^ are vectors of size N, ^ being 

the vector of shear stresses along the first pair of strips. 


Similarly, the vector of horizontal displacements of N nodes along 
the first pair of strips due to shear stresses along the second 
pair is written as 


where ^ j is a square matrix of 
stresses, -f t* T, along the second pair 


(3.78) 


the influence of shear 
of strips on displacements 


of N nodes along the first pair 



The vector of horizontal displacements of N nodes along the second 
pair of strips due to shear stresses along the first pair is 



square matrix of 
along the first 


t3.79) 


the influence of shear 
pair of strips on the 


displacements of N nodes along the second pair of strips. The 
vector of horizontal displacements of N nodes along the second 
pair due to shear stresses along the second pair is 
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where 


[ 


is a 


stresses 


' ] ‘ 


square 

along 


matrix of the influence of shear 
the second pair of strips on the 


displacements of N nodes along the second pair. 


All the matrices mentioned in Eg. 3.77 through 3.80 are 
functions of the parameters, L /B,, U /B,, v and the respective 
relative distances between the pairs of strips under 
consideration. Eqs. 3.77 through 3.80 arc combined to form a 
vector of horizontal displacements along both pairs of strips as 


(3.81) 



where 
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{ p :} = UoLoJ 

where { } = { p/‘ } + { p[‘' } 

and 


(3.82) 


Compatibility of displacements at 2N nodes is satisfied for rigid 
strips by equating the net displacements to zero, i.e. 




(3.83) 


Substituting for the vectors, a set of 2N simultaneous equations 
are obtained as 


[ ] { r/., } = { l' } 


(3.84) 


The stresses along both pairs of strips are obtained by solving 
Eq. 3.84. 

The settlement reduction coefficients of points along the surface 
due to these shear stresses are obtained as discussed in the 


earlier sections. 



3.6 Results and Discussions 

The main objective of this analysis is to evaluate the 
interaction stresses between a reinforcing strip and the soil, to 
obtain the distribution of the tensile forces in the strip, and to 
estimate the reduction in settlement due to these stresses. As 
described in the section on formulation and analysis, Boussinesg-'s 
equation for displacements due to a point load on surface and 
Hindi in^s equation for horizontal displacement due to a point load 
within the soil mass are integrated over the respective areas. 

For the purpose of numerical integration of the Boussinesq''s 
equation, the rectangular loaded area on the surface is divided 
into sub-elements of size, 0.025 B^k 0.025B^. For centrally placed 
strips the strip half -width, B , is kept constant at 0.05B, 

r f 

whereas for strips which are not centrally placed, it is 
maintained at 0-0253^.. The half width of O.OSB^ was selected in 
order to highlight the influence of the interacting shear stresses 
acting in the longitudinal direction of the strip only. 

The influence coefficient, for the horizontal 

'-j 

displacement of node i due to shear stress on clement j is 
evaluated by dividing the element j into 20x20 subareas for 
|i-j|<2. Coarser subdivision was found to be adequate for cases, 
|i-j|>2. Settlement reduction coefficients obtained considering 
elements of length equal to 0. were compared with those with 
element lengths of O.OSB^ and found to agree very well (Table 
3.1). Based on this, the lengths of the elements along the strip 
are maintained at O. IB^. The size of the elements are. 


therefore 



O.IBj«0.1Bj for centrally placed strips and 0. IB^rO. OSB^, B being 
the half-oidth of the loaded area, for strips which are off-centre. 
Thus, an element of length equal to 0. IB, would imply that a 

strip with length, L/B,=2.0, will have 20 elements along half its 
length, i.e N = 20. 

Td check the accuracy of the numerical integration, SRC 
obtained at the origin for uniform stresses at depth are compared 
with those from the exact solution given by Vaziri et. al., 
(Table 3-2) and found to agree closely. 

A parametric study is carried out in order to bring out the 
effects of length, depth, and distance of the reinforcing strip 
from the x-axis, the aspect ratio of the loaded area and the 
Poisson^'s ratio of the soil on the interaction stresses and the 
settlement reduction coefficients of points along the surface. 
The ranges of the parameters considered are 

1. : 1 to 10 (aspect ratio of loaded area) 

2. : 1 to 5 (length of reinforcement) 

3- - 0.25 to 2 (depth of reinforcement) 

4. S /B : 0.025 to 2 (distance of reinforcement from x-axis) 

y f 

5- : 0. 1 to 0-5 (Poisson^'s ratio of soil) 

6. : 0.025 for S /B^. > 0 

^ (width of reinforcement) 

: 0.05 for S /B = 0 

y f 

Except for cases where the effect of Poisson ^s ratio, is 

studied or unless mentioned, all results are for = 


0.3 



Table 3-1 Comparison of Settlement Reduction Coefficients for 


Different Element Sizes 



Table 3.2 Comparison of Computed Values of I with the Exact 

SC 

Solution for Uniform Stresses (Vaziri et. al., 19823 



I for L /B=l, 

sc r f 

B /B =0.1 

r f 

c 

0 

m 

Exact Solution 

(Vaziri et . al . 3 

Computed Values 

0.25 

0.0115016 

0.011520 

0.50 

0.0113095 

0.011315 

1.00 

0.0068968 

0.006898 

1.50 

0.0041129 

0.004113 

2.00 

0.0026271 

0.002627 



3.6.1 Shear Stresses - Single Strip 

Fig. 3.14 shows the variation of the normalised shear 
stress, T/g, with the distance, x/B^, along the strip for various 
depths of placement, U /B , and for a strip of length, L/B=2.0, 
placed below a square loaded area (L^/Bj=1.0}. The stresses in 
this figure and all other figures depicting the stresses, are the 
total stresses mobilised at the soil-strip interface. These 
stresses may be assumed to be acting in equal proportion over the 
top and bottom surfaces of the reinforcing strip. Due to symmetry, 
stresses along half the length of the strip are only considered. 
A positive shear stress is one that prevents outward lateral 
movement and as such is directed inwards. 

In the figure, as intutively felt, the shear stress at the 
centre, i.e. x/B^=0, is zero. For strips at shallower depths, 
i.e. U^/B^=0.25, the shear stresses are negative over a large 
portion of the strip. Positive shear stresses are mobilised only 
for x/Bj. in the range 0.85 to 1.1. With increasing depths of 
placement, the shear stresses are positive over most of the length 
of the strip. The sharp increase in shear stresses observed at 
the extreme end of the strip, i.e at x/B^=2.0, is because the 
strip is assumed to be rigid. This is in consonance with infinite 
contact stresses observed at the edges of rigid footings on soils. 
For depths of placement, U^/B^, upto 1.0, the shear stresses 
increase upto a distance, x/B^=1.0, along the strip and then 
decrease gradually. The increase in stresses is attributed to 
increasing displacements over the distance, x/B^=0 to 1. Larger 
the displacement, higher is the shear stress to be mobilised. 
Beyond x/B^=l, the displacements reduce for depths, U^/Bj.Sl.0, and 



hence the stresses mobilised arc also smaller. For strips placed 
at a depth, 1.0, the stresses mobilised are a maximum and 

the maximum stress, T/q=0.8 is observed near x/B^=1.0. With 
further increase in the depth of reinforcement, i.e. for 

shear stresses are smaller than those for U /B=1.0 
o > of 

over most of the length of the strip. This is because the 

displacements for depths, are smaller than those for 

U^/Bj,= l. Further, there is a gradual increase in displacements 

with distance, x/B^ for strips placed at greater depths. As a 

result, shear stresses increase monotonically with the distance, 

x/Bj., along the strip. A positive peak or a sharp increase in the 

shear stress is observed at the edge for U /B>1.5. 

o f 

The effect of the length of strip, L /B , on the shear 

r f 

stresses for a depth, U^/B^=1.0, beneath a square area (Lj./B^=1.0) 
is depicted in Fig. 3.15. For short strips, i.e. L /B,=1.0, the 

shear stresses increase almost linearly over the length of the 

strip. For strips with L /B,>1.0, it is interesting to note that 
the stresses increase to a value of about O.Sq at a distance, 
x/B^=1.0, i.e beneath the edge of the loaded area and then 

decrease monotonically. For strips of length, L /B=2.0, the 

stresses arc positive over the whole length of the strip. For 
strips longer than 2-5B^ the mobilised shear stresses become 
negative over the distance, x/Bj.>2.5. It is also worthwhile to 
note that for any length of the strip (L^>B^), the stresses are 
positive upto a maximum distance of 2.5B^ and are independent of 
the ratio, L /B . At distances, x/B >2.5 and depths, U^/B^=1.0, 
the soil displacements due to the surface loading are inward i.e. 
towards the centre. As a result the stresses mobilised are also 



Shear Stress r/q Shear Stress r/q 



Fig. 3.14 Variation of Shear Stress with Distance 
for Various Depths of Placement of Strip 
(U/B,=2, U/B,= 1, ^,== 0 . 5 ) 



Fig. 3.15 Variation of Shear Stress with Distance 
for Different Lengths of Strip 
Uo/B,=1, !'.=0.3) 



in the same direction, i.e. negative. A sharp increase in the 
negative stresses at the edge is noted. 

The effect of the aspect ratio of the loaded area, L^/B^, on 

the normalised shear stresses is presented in Fig. 3.16, for 

depth, U^/B^=1.0, and length, L^/B^=2.0, of the strip. For all 

values of L^/B^ the shear stresses increase from 0 to a maximum 

value and then decrease. The maximum stress is observed at a 

distance of about x=Bj. for all curves. The stresses are a maximum 

for a rectangular area with an aspect ratio, L^/B^=2.0, and a 

maximum stress, r/q, of 1.0 is observed at a distance, x/B=1.0. 

f 

Stresses mobilised for areas with aspect ratio, L /B =5.0, are 

f f 

slightly smaller than those for a rectangle with LyB=2.0 but are 

f f 

higher as compared to stresses for a square area (.LVB =1) or a 

f f 

strip with L^/B^=10. Stresses for an area with Lj,/B^=10 are much 
smaller than those for other rectangles over distances, x/Bj,>1.5. 
As the aspect ratio, increases beyond 2, the horizontal 

displacements start decreasing beyond a depth, U /B.=l. This 

O f 

results in reduced stress mobilisation. For areas with L^/B^=10, 
displacements at x/B^>1.5 are negligibly small resulting in very 
low shear stresses. 

Fig. 3.17 shows the effect of Poisson'^s ratio, v , on the 
normalised shear stresses, x/q, for a length, L /B =2, depth, 

r z 

U^/B^=l, and an aspect ratio, 4/Bj.=l of the loaded area. As can 
be intutively felt, the shear stresses increase with increase in 
Poisson's ratio, v . If v is small, the lateral displacement due 
to surface loading is small and hence shear stresses mobilised are 
also small. The lateral displacements increase with and are a 
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Fig. 5.16 Variation of Sheor Stress with Distance 

for Various Aspect Ratios of Loaded Area 
(L./B,=2, Uo/Bf=1. 


Fig. 5.17 Variation of Shear Stress with Distance 
for Various Poisson's Ratios of soil 

(Uo/Bf= 1 , Lr/Bf=2, Lf/Bf=l) 






maximum for p^= 0.5 (undrained case). Consequently the stresses 
mobilised are also a maximum. 

So far, the effect of the various parameters cited earlier 
on the shear stresses mobilised at the interface were studied. As 
stated in section 3.2 these shear stresses in turn produce tension 
in the reinforcing strip. The formulation for the computation of 
tensile forces is dealt with in section 3.2.2. 

The variation of the normalised tension in the reinforcement 
strip, T/qB^r along the length of the strip for lengths, L /B =2, 

I r f 

and aspect ratioy of the loaded area and various values 

of depth ratio, U /B^, are depicted in Fig. 3-18. A negative 

o f 

2 

value of the normalised tension, T/qB^, indicates compression. It 

is clear from this figure that for a depth, U /B,<0.35, the 

reinforcing strip is subjected to compression throughout its 

length. Strips at a depth, U /B=0.5, are subjected to compression 

o f 

beyond a distance, x of 1*4B^. The maximum tension at the centre 

is observed at a depth, U /B=1.0 and is equal to O. 105qB^. This 

maximum value decreases with further increase in depth, i.c. for 

U >B - It is also evident from the figure that strips placed 
o f 

beyond a depth of 0-5B^ are subjected to tension over their full 
lengths. This is corroborated in Fig. 3.14 which shows the effect 
of depth on shear stresses. The stresses are predominantly 
negative for strips at depth, U^/B^^O.35, while they are positive 
for strips at greater depths. Consequently, strips at shallower 
depths (U /B<0.35) are ineffective in providing restraint to the 


soil. 



The effect of the various parameters mentioned at the outset 


of this discussion on the mobilised interaction stresses and 
tension in the strips are brought out through the figures 
illustrated so far. It is due to these stresses that the surface 
settlement is affected. The improvement in the- soi 1 -foundat ion 
response due to strip reinforcement can well be appreciated if the 
effect of the various parameters on the settlement reduction 
coefficients (SRC) are studied. Hence more emphasis is laid in 
presenting the effect of the considered parameters on the SRC. 
From the study of the stresses it is more or less clear that 
strips placed in the depth range (U^/B^) of 0.75 to 1.5 below 
rectangular areas are most effective. 

3.6.2 Settlement Reduction - Single Strip 

To study the effect of shear stresses on the SRC along the 
surface, it is essential to look into the basic solution given by 
Mindlin for vertical displacements due to a horizontal force at 
depth. The vertical displacement of any point, B(x,0,0) (Fig. 3.4) 
is given by Eg. 3.27. Non-dimensional ising all parameters with the 
depth, c of the load, Eq- 3.27 is written as 

Q c - (3.85) 


where G =E /2(l+v ) and all parameters are the same as defined in 
s s s 

Eq. 3.27. I is an influence coefficient dependent on x/c and 

Z 

The plot of I versus the normalised distance, x=x/c, for 
2 

different values of Poisson^s ratiOy is shown in *9 

For all valuaa of f the coefficient, I. increases from zero, 

& 



tension 


0 . 00 - 

0 . 



Fig. 3.18 Variation of Tension with Distance for 
Various Depths of Placement of Strip 
(U/B,=2. L,/Bf=1, ^5=0.3) 


Fig. 3.19 Influence Coefficient for Vertical Surface Displacemen 
due to Horizontal Force at Depth (Mindlin Solution) 





reaches a iriaxirtuiTi valuer and then decreases. For i> SO. 3, I 
becomes negative at larger values of x/c. In all cases I tends 

z 

to zero at very large values of x/c. For undrained soils iv =0.5), 

s 

values arc positive for all distances indicating surface heave 

throughout. The maximum value of I ranges from 0.015 for v =0.1 

* s 

to 0.03 for v^=0.5. The distance x, at which the maximum heave 

occurs increases from 0.5 for v =0.1 to about 0.7 for v =0.5. The 

s s 

distance, x, beyond which settlement occurs for v <0.5 is again a 

s 

function of v . For <0.4, distances beyond which points settle 
Ci.e. negative I > get closer to the point force. 

Z 

Fig. 3.19 gives an insight into the fundamental solution 
given by Hindi in. It is this equation that has been used to 
calculate the SRC on the surface for the various shear stress 
distributions arrived at from the interaction analysis. This 
figure will help in understanding the variation of SRC due to the 
stresses presented in the following figures. 


Fig. 3.20 depicts the variation of SRC along the x-axis on 
the surface for strips of various lengths, and for a depth of 
placement, U /B=1.0, beneath a square loaded area (L /B=1.0). 

Of II 

x,/B,=0 is the origin 0, or the centre of the loaded area, while 

x^/B^=1.0 is the edge. For all lengths of reinforcement, the 

settlement reduction is maximum at the centre of the loaded area, 

0. SRC reduces with increasing distances from the centre i.e. 

with X /B . The SRC increases with increasing lengths of the 
f f 

strips. The improvement in SRC for the length ratio, L^/B^ 
increasing from 1 to 2 is significant. As the length ratio, 


increases from 2 to 5, the increase in SRC is insignificant. The 



SRC at x^/B^=0, for L^/Bj.-2, is 0.. 00335 while it is 0.0035 for 

L^/B^=5. Thus at a depth, U^/B^=1.0, no advantage accrues by 

providing reinforcements of lengths greater than twice the width 
of the loaded area. A strip of length equal to B^ gives a SRC of 

0. 0023 at the origin and -0.0003 at the edge, the latter 

indicating a small additional settlement. For strips of lengths 
equal to 2B^ to 5B^, the SRC at the origin is around 0.0034 and at 
the edge around 0.001. The reasoning for this trend is explained 
as follows. The positive shear stresses along the strip of length 
equal to cause the origin to move up while they try to pull the 
edge down. Hence a small settlement at the edge is observed in 
Fig. 3.20. For strips with longer lengths, stresses beyond a 
distance, x/Bj.=2.5 (Fig. 3.14) are negative. The additional 

positive stresses (as compared to those for L^/Bj,=1.0) beyond a 
distance, x/B^=1.0, help in pushing the points from the origin to 

the edge further upward and hence the 9RC values at all points, 

1. c. x^/B^ varying between 0 and 1.0 are higher as compared to 
those for L /B =1.0. 

r f 

Fig. 3.21 presents the effect of depth, ^he SRC 

for strips of length, L /B =2, and for a rectangular loaded area 

r I 

with aspect ratio, L^/B^=2. It is seen that except for U^/B^.— 0.25, 
for all other depths the SRC is maximum at the centre of the 
loaded area, i.e. at x^/B^=0 and decreases gradually with 

increasing x /B . It is interesting to note that for a depth 

f f 

ratio, U /B=0.25, the maximum SRC is observed at a distance, 
of 

x/B=0.75. This is because, for this depth the stresses are 
f f 

As a result, the points away from the 


predominantly negative. 









oriQir' on the surface are pushed up by a larger anouTrt as compared 

to the centre. Hence SRC increases upto a distance of x/B=0.75 

f f 

The negative stresses beyond a distance, x/B^=l along the strip 

cause the edge of the loaded area (Xj./B^=l) to move downwards and 

hence SRC decreases beyond Xj,/B^=0.75. With increasing depths, 

the SRC are positive and the rate of decrease in SRC with x /B 

f f 

decreascs- 

The effect of aspect ratio of the loaded area, L^/B^., on SRC 
for different depths, U /B,, and for strips of length, L /B =2.0, 
is presented in Figs. 3.22 and 3.23. In Fig. 3.22 the plot of SRC 
versus the distance, along the surface is shown for strips 

at depth, U^/Bj.=1.0- For all values, SRC is maximum at the 

centre and decreases gradually with increase in 

interesting to note that SRC values increase with L^/B^ and attain 
maximum values for rectangles with aspect ratio Lj./B^=2.0. With 
further increase in L^/B^., the SRC values decrease. The maximum 
SRC at the centre of the loaded area or origin is for a rectangle 
with an aspect ratio, Lj./B^=2, and it equals 0.00435 while the 
least improvement is for a square which equals 0.00335. 

For strips of length, at a larger depth of U^/Bj.=2, 

the variation of SRC with distance, >{j./Bj., along the surface is 
presented in Fig. 3.23. In this figure, as compared to the 
previous one, the variation in the SRC with or with the 

aspect ratio of the loaded area, (L^/B^) is smaller. Further, SRC 
is maximum for a longer rectangle with L^/B^=5. The SRC at the 
centre of the loaded area, i.c. at Xj/Bj=0, for rectangular areas 
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Rg. 3.22 Voriotion of SRC with Distance along Surface 
for Different Aspect Ratios of Loaded Area 
(U/Bf=2, Uo/B,=1, i^s=0.3) 
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Rg. 3.23 Variation of SRC with Distance along Surface 
for Different Aspect Ratios of Loaded Area 
(U/B/=2. Uo/B,=2. >'s=0.3) 




with is 0.00385 and that ■for a square (L^/B^=l) is 

0.00185. 


Fr om 

these two 

figures. 

viz. 

Figs. 3.22 

and 3.23, 

the 

settlement 

reduction 

for strips 

of length, L /B =2, 

r 1 

is observed 

to 

be maximum 

if placed 

at intermediate 

depths, i.e. 

for U /B, 
o f 

of 


around 1 beneath rectangular areas with aspect ratio (L /B ) of 2. 

f f 

At shallower depths the negative shear stresses are responsible 
for low SRC values whereas at greater depths the shear stresses 
are too deep to bring about any significant improvement in the SRC 
values. These results emphasize the fact that strips of length, 
L /B =2, placed at intermediate depths (U /B =1.0) are 

r f Of 

advantageous as compared to those placed at very shallow or at 
great depths- 

Thus far, the variation of SRC with the distance, 
along the surface was studied and the effects of the various 
parameters on SRC were looked into. From this study it is clear 
that except for strips placed at very shallow depths, the centre 
of the loaded area, 0, experiences the maximum settlement 
reduction. In the next couple of paragraphs the effect of the 
parameters on the SRC at the centre of the loaded area defined as 
I„ is analysed. 

Sc 

Fig. 3.24 shows the variation of with depth, fo*' 

various aspect ratios of the loaded area and for a strip of 
length, L/B=2.0. It is observed that for square 

r f 

tL/B=1.0), I increases with depth and reaches a maximum 
f f sc 


value 



A siiriilar trend is observed for 


of 0.0034 at a depth of 0.9B . 

{ 

areas with L /B =2.0 for which the maximuin value of I is 0.0044 
* ^ sc 

occuring at a depth of 1.1 In the case of areas with aspect 

ratio of 5. the versus U^/B^ curve is flatter for depths upto 

0-5B beyond which the values of I increase upto depths in the 
range 1.4B to 1-5B and then decrease with U /B . A maximum 

II Of 

value of I equal to 0.00465 is observed for areas with L /B =5 

SC. f f 

at a depth of 1 . 4B . In case of long rectangular areas (LyB=10), 

1 f f 

the reduction in settlement due to a single strip is less than the 

reduction for an area with L,/B=5.0. The depth of the strip to 

f f 

cause maximum reduction in settlement however, is more C=1.5B^). 

These trends can be explained with the help of figure 3.25. 
At a depth of 0.5B^ a maximum positive stress of 0.35q and a 
maximum negative stress of -1.4q are observed for areas with 

L^/B^=5. On the other hand, areas with Lj,=10Bj. show a maximum 

positive stress of 0.2q and a maximum negative stress of -1.9q. 
These stresses lead to a lesser settlement reduction for areas 
with Lj./Bj.= 10 as compared to areas with L^/Bj.=5, at the same depth 
of O.SBj.. 

The change in the rate of variation of with 

areas with L /B =5 and 10 at shallow depths can be explained with 
f f 

the help of Fig. 3-19. It can be seen from this figure that for 

V =0.3, the coefficient I is negative for distance, x>3c 

indicating settlement of soil. This implies that for a strip at a 

depth of 0.25B^, negative stresses beyond a distance of 0.75B^ 

(3x0. 25B ) push the soil at the origin upwards while for a strip 
f 
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Fig 3 24 Effect of Depth of Placement of Strip on 
for Different Aspect Ratios of Loaded Areo 
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Fig, 3,25 Variation of Shear Stress with Distance 

for Different Depths of Placement of Strip 
(U/B,=2, J^s=0.5) 




at a depth of 0.5B^, the negative stresses beyond a distance of 

1.5B^ cause upward movement of soil at the origin. From Fig. 3.25 

it is evident that for an area with Lj,= 10B^r the negative stresses 

for a strip at a depth of 0.25B^ arc much larger than those for a 

strip below an area with an aspect ratio, L^/B^=5, and hence give 

rise to a higher I value. 

sc 

Fig- 3.26, 3.27 and 3.29 present the variation of I with 

sc 

length of the strip, various values of the aspect 

ratio, L,/B,, of the laoded area. The effect of LVB on I for 

f f 1 f sc 

a strip placed at a depth U^/B^=0.25 is depicted in Fig. 3.26. For 
all L,/B, values, the I increases gradually with increase in the 

f f SC 

length, L /B., of the strip. I for a short strip (L/B=l), for 

r f SC r f 

all L,/B, values is negative. It is also evident that the I 

f f sc 

values increase with increase in the aspect ratio, It is 

worth noting that for smaller aspect ratios of the loaded area 
tL,/B,<2) increase in the length of strip does not improve the I 

f f SC 

values whereas for longer rectangular areas (L^/B^>3) longer 

strips show higher I values. For square loaded areas (L /B =l)r 

sc * * 

I r for strips of length, L /B =1 and 5 are -0.0008 and 0.0012 
sc r f 

respectively. For longer rectangles (L /B =10), I for strips of 

I t SC 

length, L /B =1 and 5 are -0.002 and 0.006 respectively. 

r f 

In Fig. 3.27, the plot of I^^ versus the length, of 

the strip for different aspect ratios, L^/B^., of the loaded area 
and for a depth ratio, U^/Bj.=l, is presented. For L^/B^<3, 
increases upto a length, L of 2Bf. With further increase in 

r 

L /B , no significant increase in I values is observed. In case 

IT' I SC 
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Fig. 3.26 Effect of Length of Strip on' Ijc for 

Different Aspect Ratios of Loaded Areo 
(Uo/B,=0.2^ i's=0.3) 
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Fig. 3.27 Effect of Length of Strip on Isc for 

Different Aspect Ratios of Loaded Area 
(Uo/Bf=1, i's^O.3) 


of areas with L^/B^>3, increases for strip lengths upto 2B^, 

is almost constant for L^=2Bj, to L^=3B^, and then shows a rapid 
increase beyond this length. Further, for L/B=3.0, I values 

f f SC 

are a maximum. There is a drop in I values with further 

£» O 

increase in L /B values, while the rate of increase in I values 
’ • sc 

beyond Lj./Bj.=3 is higher. 

Fig. 3.28 helps in explaining these patterns. It is 
observed frowt this figure that beyond a distance of 2.7B^ for 
areas with L /B =3, and 2.25B for areas with L =10B , the shear 

SI I 1 X 

stresses are negative. The negative stresses act over longer 
lengths for longer strips. Again, from Fig. 3.19 it is evident 
that at a depth, U =B , negative stresses beyond a length of 3B 

Of I 

cause the origin on the surface to move up, thus facilitating 
settlement reduction. Hence a continuous increase in I values 

SC 

with increase in the lengths of the strip for is observed. 

Fig. 3.29 presents the variation of I with the length, 

SC 

L /B„, of the strip for strips placed at a depth, U /B =2. In 

this plot, unlike the earlier plots no significant improvement in 

I with increase in L /B beyond 3 and for all L./B. i'S observed, 
sc r f f 1 

Further, I is maximum for L /B =5. The value of I for strips 
sc f f 

of length, L /B =5 and for aspect ratios of the loaded area, 

r f 

L/B^=l, 5 and 10 are 0.0024, 0.0052 and 0.0050 respectively. 

From the results presented in Figs 3.27 and 3.29 it can be 
noted that longer strips placed under longer rectangles prove more 
beneficial. Increasing the length of the strips beyond 2.58^ 
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Fig. 3.28 Variation of Shear Stress with Distance for 

Different Aspect Ratios of Loaded Area (Uo/Bf=1) 
a) b) U/B,=4. c) U/Bt=5 


at greater 


under area?, with smaller aspect ratios (L^/B^<2) and 
depths is not advantageous. 

The effect of Poisson''s ratio, v , on I is studied in 

Fig. 3.30 which presents the variation of I with the length, 

L /B , of the strip for different values of p and for the strip 

placed at a depth, U^/B^=l, below a square area (L^/B^=l). It is 

evident from the figure that with increase in p , there is an 

6 

increase in the values. For ^^ 2 : 0.2 there is no significant 

improvement in I values for strips with lengths, L /B >2.5 or 3. 

W Cl T T 

This again reiterates the fact mentioned earlier that longer 
strips (L /B.>2.5) under smaller areas do not prove advantageous. 

Jr I 

The increase in I values with increasing values of p is due to 

sc s 

the fact that at a depth, U /B=l, the stresses are predominantly 

o f 

positive. Higher interfacial stresses are possible for soils with 
higher values of Poisson's ratio since lateral displacements are 
higher. As a result, soils with p =0.5 show maximum I values. 

So far, single strips placed centrally below the loaded area 
were studied. Strips placed away from the centre and in groups 
also affect the SRC considerably. In the next two sections, 
effects due to two strips and effects of multiple strips are 
studied. 

3.6.3 Two Strips at Depth U 

o 

The half-width of the strip in this case is 0.0256^,. A 
single strip of half-width, 0.058^,, placed centrally below the 
loaded area, as was studied in the previous section is equivalent 
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Fig. 3.29 Effect of Length of Stop on Isc for 

Different Aspect Ratios of Loaded Area 
(Uo/B,=2, i^s=0.3) 


? 0.006 


Fig. 3.30 Effect of Length of Strip on Isc for 

Different values of Poisson's Ratio of Soi 

(Uo/Br=1. U/B,= 1) 





a distance 


r S , 
y 


to two strips of half -width, 0.0256^., placed at 

equal to 0.0258^., i.e. the half-width of the strip, from the 

X-axis. At this spacing the two strips are contiguous. Settlement 

reduction coefficients obtained for two strips placed touching 

each other at the centre are compared with those obtained for a 

single strip (Table 3.3). The results agree well. Table 3.3 

implies that a single strip of half -width equal to 0.05B and two 

f 

contiguous strips of half-width, 0.025B give the same I values. 

f sc 


Table 3.3 I Values for a Strip Considered as Single and 
sc 

as Two Contiguous Ones 



I for U /B, 

sc o f 

=1, L /B,=l, V =0.3 

If s 


L /B 

r 1 

Single Strip 

B /B =0.05 

r f 

Two Strips 

B /B,=S /B, =0.025 

r f y f 

Difference 

7. 

1.0 

0.002296 

0.002390 

4. 10 

2.0 

0.(X>3357 

0.003474 

3.50 

3.0 

0.003381 

0.003495 

3.46 

f 

O 

0.003417 

0.003536 

3.48 

5.0 

0.(X)3480 

0.003600 

3.45 


It is to be noted that SRC values presented in the figures 

are the settlement reduction coefficients due to one of the strips 

among the pair at a distance, S . To obtain the SRC due to both 

y 

the strips, the SRC values from the figures are to be doubled 
since the strips are symmetrically placed about the x-axis. 

Fig. 3.31 shows the plot of the normalised shear stresses, 
T/q, versus the distance, x/B^, along the strips for different 



values of the distance. 


X-axis and 


strips from the 

for strips of length, L/B~2, placed at a depth, U^/B,=l, beneath 
a square area. For all S^/B, values the stresses originate from 
0, and reach a maxitnum at a distance, x/B^ of around 1. An 

interesting point worth noting is that maximum stresses are 
observed for the distance, S^/B^=0.25, while strips placed closer 
to the x-axis (S^/Bj.=0.1) show lesser values of stresses. The 

maximum shear stress, r/q, for S^/B^=0.25 is of the order of 1.13 
at a distance. x/B^=l. Further, if this figure is compared with 
Fig. 3-14 and the curve of shear stresses for a single centrally 
placed strip is superimposed on Fig. 3.31, it is seen 

that the shear stresses for the single strip case lie slightly 

below the curve for S^/B=l. As S^/B^ increases beyond 0.25 the 
stresses decrease quite rapidly, the stresses for the distance of 
the strip, S /B , from the x-axis equal to 2 being very small. 

jr A 

This trend of the stresses can be explained as follows. For 

strips placed closer to the x-axis (S /B <0.25) the interference 

y f 

of one strip on the other is high due to close spacing. For 

^y/®f.~0»025 (contiguous strips) the interference is maximum. As 

this distance increases from 0.025 to 0.25, the influence of one 

strip on the other decreases. Consequently there is an increase 

in the stresses- Beyond a distance, S , of 0.25B, the lateral 

y f 

displacements due surface loading decrease resulting in lower 

shear stress mobilisation. For large S /B, values (S /B >1.5) the 

y f y f 

lateral displacement is negligibly small and hence stresses arc 
relatively smaller compared to closer spacings. 


The 

variation of 

the 

normalised 

tension, T/qB*, with 

the 

distance. 

x/B , along 

the 

strip for 

different values of 

the 
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distance, of the strip from the x-axis and for strips of 

length, L^/B^ 2, placed beneath a square area (L^/B^=l) at a 

depth, 1 , is depicted in Fig. 3.32. The stresses being 

positive over the full length (Fig. 3.31), the strip is subjected 

to pure tension throughout its length. As a consequence of these 

stresses, tension is maximum for strips placed at a distance, 

Sy/B^=0.25, and minimum for Sy/B^=2.0. The maximum tension, 

T/qB^, at the centre for strips at a distance, S /B=0.25 is 0 076 
' y f 

while the same for strips at a distance, S /B =2 is 0.014. 

y f 

The effect of the distance of the strips from the x-axis 

(S /B ) on I for various depths, U /B,, and for strips of 

length, L^/B^=2, placed beneath a square area is presented in Fig. 

3.33. For closer spacing from the x-axis, i.e. S <0.25B , I 

y f so 

increases monotonical ly upto a depth, U =B, and then decreases 

o f 

with increase in depth upto 2B,- I for a distance, S =0.25B 

f sc y f 

increases from 0.00025 at a depth of 0.25B^ to 0.0023 at depth of 
B and then decreases to 0.0013 at U /B=2. As the distance, S , 

f of y 

increases beyond 0-25B^, there is a change in slope at shallow 
depths and a decrease in the maximum value of I . The depth at 

SC 

which maximum settlement reduction occurs also increases with the 

increase in the distance, S /B . The I values decrease with 

y f sc 

increase in S /B since the stresses decrease with increase in 
y t 

S /B (Fig. 3.31). In reasoning for slope changes at shallow 
depths, apart from the aspects considered in Fig. 3.19, the ratio 
S^/B^ also plays an important role. It is interesting to note 
that larger settlement reduction results from relatively distantly 
spaced strips placed at greater depths (U^/Bj.>l .5) . 



lension T/qB, 


Fig. 3.31 Variation of Shear Stress with Distance for 
Different Distances of Strip from x— axis 
(U/B,=2. U/B,= 1. Uo/B,= 1. i/s=0.3) 



Fig. 3.32 Variation of Tension with Distance for 

Different Distances of Strip from x— axis 








Fig. 3.34 presents the variation of I with the distance, 

sc 


s /B , for different lengths 

V f 

of 

the strip. 

L /B . 

r f 

placed 

at a 

depth, U^/B^=1.0, below a square 

loaded 

area. 

It is 

observed 

that 

the I values increase with 

sc 

S /B, upto 
y f 

S /B^ 
y f 

equal 

to 0.25. 

For 


larqer values of S /B,, I values reduce sharply with increase in 

S /B At distances beyond 1.75B , the settlement reduction is 

^ I 

negligible. For all lengths (L^/B^) a maximum is observed at 

a distance of S /B^=0.25. This is because the shear stresses for 

strips at this distance are maximum (Fig. 3.31). The variation in 

I values with S /B, is a direct outcome of the variation of 
sc y f 

stresses- The maximum I values are 0.0016 for L /B =1, 0.0023 
for L /B =2 and 0.0024 for L/B=5.0. The increase in for 

r f T t 

L/B>2.0 is minimal which indicates that there is no advantage in 

r f 

providing strips longer than 2B^ at a depth of 


In studying the effect of S^/B^ on settlement reduction it 
is imperative that for maximum benefit from strips of length, 
L /B =2, placed at a depth, U /B =l, under square areas, the 

r f j 

optimal choice of the distance from the x-axis, S^/B^ is 0.25. 

3.6.4 Two Strips One Below the Other 

The interference of two strips one below the other and the 
consequent effect on the shear stresses mobilised are studied here 
under. The settlement reduction resulting from these stresses are 
compared with those obtained from the principle of superpositi 
i.e. I values obtained by considering each strip as a s' g 

SC 

strip. The ratio of the I obtained for the present case 

SC 

that obtained from the principle of superposition is ter 







ress variation for two strips 


Fig. 3.35 depicts the shear st 

of length, L/B^=2, at depths, U^/B^=0.5 and U^/B^=1.0 below a 

square loaded area. The dotted lines are the shear stresses for 

single strips at the respective depths. It is evident that the 

stresses for the two strips, one below the other are smaller than 

those for single strips. The difference in the stresses is much 

higher for the strip at a shallower depth (U/B=0.5). The 

o f 

decrease in the shear stresses is due to the interference of one 
strip on the other. As a consequence the settlement reduction 
coefficients of points on the surface are smaller. 

Table 3.4 presents the values for two strips placed one 
below the other at depths, ^nd for various lengths, 

of the strips. Two sets of results are presented. In the 
first set, one strip is placed at a depth, U^/B^=0.75 and in the 
second it is placed at a depth, U^/B^=1.5. The depth, 
both sets is varied from 0.25 to 2. It is seen that the least 

value of I in the first set is when the second strip of length, 

0 

L /B =2 is placed at a depth, U /B =1 and equals 0.769. The 

r t if 

corresponding value for the second set is for a strip of length, 
L^/B^=5 at a depth, U^/B^=2 and equals 0.825. In both the sets it 
is observed that at shallow depths the ratio, I^, is marginally 
greater than 1 which indicates that the I^„ values obtained by 
considering the interference of each strip on the other is greater 
than that obtained from the principle of superposition. This small 
increase is because of the following reason. The shear stresses 
on the second strip at shallow depths (U^/B^^O.25) are 
predominantly negative. These negative stresses cause the points 



Table 3.4 Values for Two Strips Placed Centrally Below the 

Other 



U 

o 

t 
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11 

m 

75 


U /B 
o 

f = 

5 

\ L /B, 

N. r I 

1 

2 

5 

\l /B^ 

N. r 1 


2 

5 

U /B^ 




1 t \ 




0.25 

0.956 

1.1 

1.09 

0.25 

0.99 

1.04 

1.1 

0.50 

0.800 

0.8 

0.815 

0.50 

0.976 

0.96 

1.02 

1.00 

0-818 

0.769 

0.785 

0.75 

0.962 

0.917 

0.94 

1.50 

0,962 

0.917 

0.942 

1.00 

0.926 

0.86 

0.863 

2.00 

O 

t 

0.965 

0.984 

2.00 

0.93 

0.858 

0.825 


on the first strip at depth to move away from the centre. 

Hence, the net horizontal displacements of points along the first 
strip are higher than the corresponding displacements for a single 
strip. As a result the stresses mobilised are higher. 
Consequently, the settlement reduction of points along the surface 
are also higher. This results in value greater than 1. 

3.6.5 Multiple Strips 

Settlement reduction due to two pairs of strips at the same 

depth is studied herein. I _ values due to these two pairs, one 

separated by a distance, 2S^^/Bj., and the other by a distance, 

2Sy^/B^, are computed. values for strip pairs at S^/B^ and 

S /B, from x-axis are computed independently and summed up. The 
y2 f 

ratio (I ) of I for the first case with respect to the one 

I> SC 

obtained due to summation is calculated. The difference of 

two spacings (S -S ) is termed as S . It is checked whether 

y2 yi ^ 





the principle of superposition is valid for multiple strips. 

The variation of shear stresses, r/q with distance, x/B for 

f 

multiple strips of length, LyB=2, placed at a depth, U /B^=l, 
below a square area is presented in Fig. 3.36. The solid lines 
indicate shear stresses for multiple strips, placed at distances, 
S /B =0.025 and S /B =0.5, and the dotted lines are the shear 
stresses for a single pair of strips at the respective distances, 
S^/Bj^, from the x-axis. The shear stresses mobilised for multiple 
strips are smaller than those for a single pair of strips. 
Stresses on one pair influence the displacements of nodes along 
the second pair resulting in smaller net displacements of the 
nodes. As a result the stresses mobilised are smaller. These, in 
turn affect the settlement reduction of points along the surface, 
as can be seen in the following figures. 

Fig. 3-37 depicts the variation of I with S for different 

r* yD 

lengths, L /B,, of the strip placed at a depth, U /B =1, below a 
r f Of 

square area. I values increase from about 0.54 to 0.99 with 
n 

increasing values of S . For strips with length, L /B =1, I 

y0 r f D 

equals 0.9 for S =0.5, whereas for strips with L /B ^2, I equals 

yl> r f D 

0.9 for S =1. It is clear that for multiple strips the principle 
y» 

of superposition is not valid unless S is greater than about 0.7 

yD 


for 

L /B =1 and around 1 

T f 

for L /B >2. 

r f 

When the S is 

yD 

smal 1 , 

all 

the 

four strips have 

equally high 

influences on 

each 

other 


resulting in low shear stresses. 


Fig. 3.38 shows the variation of with 


for different 



Fig. 3.37 Effect of Difference of Spacing on lo 
for Different Lengths of Strip 
(Uo/B,= l. U/B,= 1, ^'s-0.3) 


Uo/B, 






beneath 


depths of placement, of strips of length, L /B =2, 

r f 

a square area. I increases from 0.53 to 0.99 with S 

yD 

depths. With depth, U^/B^ increasing from 0.5 to 1 the decrease 
in the values is noticable while it is marginal with further 
increase in depth. The principle of superposition appears to be 
valid for S >0.7 for U /B <1 and for S >1 for U /B >1. 

y*> o J yl> 0 £ 

From these two figures it is clear that for close spacing of 
strips the principle of superposition is not valid. An appropriate 
correction factor based on the spacing of strips is to be 
considered in order to obtain the settlement reduction 
coefficients due to multiple strips from those of single strips. 

3.7 Conclusions 

A parametric study, using the elastic continuum approach has 
been developed to study the effects of the aspect ratio of the 
loaded area, length and depth of a strip reinforcement on the 
shear stresses mobilised at its interface with soil and on the 
reduction of surface settlements. The effects of the distance of 
the strip from the centre and the Poisson's ratio of the soil on 
the quantities mentioned above have also been examined. The shear 
stresses mobilised at the interface oppose the displacements of 
the soil due to the surface load. It is perceived that these 
stresses are responsible for the reduction in surface settlements- 
The results obtained from the parametric study depict that an 
optimum choice of the parameters has to be made in order to 
extract the maximum benefit from rigid strip reinforcements. 
Based on the results, the following conclusions are enumerated. 



1. A short strip placed at a shallow depth 0.35) is 

ineffective in reducing settlement at the surface. It is 
subjected to compression since the shear stresses mobilised 
are negative, i.e., they arc directed outward. 

2. Maximum settlement reduction of the centre of the loaded area 
due to a single strip is observed for a strip placed at a 
depth in the range of 0.758^. to below an area with aspect 
ratio upto 3, and at about 1.5B^ below longer rectangles. 

3. For a strip placed at a depth equal to B^ below a loaded area 
with aspect ratio upto 3, the settlement reduction coefficient 
increases appreciably with the length of strip upto L =2B . 

r f 

With further increase in the length of the strip the increase 
in SRC is marginal. 

4. As expected, the results show that a longer strip (L /B >3) is 

r f 

effective if placed below a longer rectangular area upto a 
depth equal to 

5. At depths, U^/Bj,2:2, increasing the length of the strip beyond 
3B^, for any aspect ratio of the loaded area does not 
contribute to any increase in the settlement reduction 
coefficients. 

6. As expected, strips placed in incompressible soils (saturated 
undrained condition) show maximum reduction. 

7. For a pair of rigid strips of length, L^/Bj,=2, placed 
symmetrically about the x-axis and at a depth, U^>0.5B^, 
maximum reduction in surface settlements is observed if the 
centre of each strip lies at a distance of 0.25B^ from the 
X-axis. 

8. For two strips placed one below the other the principle of 



superposition is valid if the vertical spacing of the strips 

is greater than 0.5B^. A suitable correction factor will have 
to be used for closer vertical spacings. 

9. For multiple strips at the same depth, the principle of super- 
position is valid if the difference in the distances of each 
pair from the x-axis is greater than 0.7B^. For closer 
spacings (S^/B^<0.7) a suitable correction factor has to be 
used. 



CHAPTER 4- 

SOIL—RICID STRIP REINFORCEMENT INTERACTION “ NORMAL STRESSES 

4.1 Introduction 

In chapter 3, the rigidity of the strip in the axial 
direction was considered. In this chapter the strip is assumed to 
be rigid in the vertical (normal) direction. As seen earlier, the 
applied load on the surface causes the soil below to move downward 
and outward. A deformable strip will deform as shown in Fig. 4.1 
(dotted line) but a strip assumed to be rigid undergoes a rigid 
body displacement in the vertical direction. As a result, normal 
stresses are mobilised along the soil-strip interface. The 
stressess are symmetric about the z-axis due to the symmetry of 
the loading. Since the normal stresses alone are considered it is 
assumed that the rigidity of the strip in the horizontal (axial) 
direction (i-e- the shear interaction) has no effect on the normal 
displacements. 

The mobilised inter facial normal stresses represent the 
difference between the normal stresses on the top and the bottom 
faces of the strip. Because of the rigidity of the strip, the 
normal stresses act upward near the centre and downward at the 
sdges, thus forcing the strip to undergo a uniform displacement. 
The net result of these stresses is to push the soil on the 
surface upward near the centre of the loaded area. Consequently, 
there is a reduction in the settlement of the surface. Thus, the 
reinforcing strip helps in reducing foundation settlements. 



4.2 Single Strip 


4,2.1 Problem Definition 

The analysis of soi 1 -reinforcement strip interaction for 
normal stresses is very similar to the analysis for shear stresses 
discussed in chapter 3, but with a difference in the mechanism. 
In this chapter the normal stresses only are considered. A rigid 


reinforcing strip 

of size, 2L x2B , is 

r r 

placed at 

a depth. 


centrally 

below a 

rectangular area of 

size. 

2L x2B, 

t f 

, transmitting 

a uniform 

load of 

intensity, q, (Fig. 

4.21. 

The 

width of 

the 

strip (0. 

IBj.) is 

relatively small 

and 

its thickness. 

t , 

r 


negligible. 


4.2.2 Analysis 

The elastic continuum approach is resorted to in modelling 
the soil -strip interaction. The vertical displacements of the 
points at the soil -strip interface are calculated using 
Boussinesq''s solution for vertical displacement at a point within 
an elastic half-space due to a point load acting on the surface. 
The Boussinesq''s problem is as shown in Fig. 4.3. The vertical 
displacement of point. A, as shown in the figure, due to load P is 


P(l+P ) 
s 



s 




(4. 1) 


where R=y r* + z*, r and z arc respectively the radial distance 
and depth of point A, with respect to the load P , and are 
respectively the modulus of deformation and Poisson^'s ratio of the 


soil. 
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Fig. 4.1 Strip Deformation 



The vertical displ 
strip, due to thr 
integrating fq. 4. 



acemrnts, p^, of points along 
uniform surface load, 

1 over the loaded area as 


q ( 1 +!-> ) 
e: 

SnE'-R 
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2ti-i> y + 


dA 


the 

is 


reinforcement 
calculated by 


(4.2) 


wtiere dA is an c'lenriental area on the loaded surface. 


The discretization of the loaded area is the same as in 
section 3.2.2. The vertical displacement, , of the centre of 

Z I 

element i (Fig. 4.4), due to the surface load is obtained by 
carrying out the integration of Eq. 4.2 numerically as 




qtl+y ) 
s 

2'nL~R 

s 


2(1 -n^) 



dA 


(4.3) 


where y =■/' x * 4 ^ y*, and R=i/ r ^ . NDn“diiriensional ising all 

o 

length parameters with half the width, B^., of the loaded area, 
Eq. 4.3 is expressed as 



t( 

I q 


(4.4) 


where is a dimensionless influence coefficient that depends on 

i 

the aspect ratio, of the loaded area, the depth, U^/B^, of 

the strip, the Poisson^'s ratio, of the soil and the location 

of node i. The vector of the vertical displacements of all the N 
nodes along the half-length of the strip is obtained by evaluating 




Eq. 4.4 for all the nodes as 



(4.5) 


where the vectors 


{ p' } and I l' j 


are of size N. 


Eq- 4.5 gives the vertical displacement profile of the soil 
without any reinforcement at depth if a rigid strip is placed 
at this depth, normal stresses arc mobilised at the soil-strip 
interface so that the strip undergoes a rigid body displacement. 
To obtain these stresses, the vertical displacements due to these 
stresses are to be computed. It is achieved by using Mindlin*^s 
solution for vertical displacements due to a vertical force acting 
within an elastic half-space. The vertical displacement of a 
point, A(x,y,z) due to a vertical force, Q, (Fig 4.5) acting at a 
depth, c, in a semi -inf ini te elastic continuum is 


Q(l+w ) 


* BrrE 




s m 


(4.6) 


(3-4v ) C8(l-v )^-(3-4v )] (z-c)* 

where G = +■ ? — + + 


R 


R 


R‘ 


(3-4v )C(z+c)^-2cz3 6cz(z+c)* 

f + 

R* r" 


R^= -/^ 


* + y*+ (z-c)* , R= -/x^ + y*+ (z+c)* 


7 are the co-ordinates of 


and X jr y I 
displacement of any point 
on an elemental area, dA, 


point A. The vertical 
along the strip due to normal stress, a, 
is 


dp^ 


xydACl+v ) _ 

6nr ( T^T 

K C 


(4.7) 


To obtain the vertical displacements, pj along the 

Z 

reinforcing strip due to the mobilised normal stresses. Eg. 4.7 is 
integrated over the area of the reinforcing strip. 


The vertical displacement, p . of node i, due to normal 

ZIJ 

stress. O', acting on element j (Fig. 4.6) is 

j 

p (1 +1-' ) r ‘v 

= j §ng ; n^ [ G ] c, dA (4.8) 

dA 

where all the parameters are as mentioned in Eg. 4.6 with c=U^= 2 , 
and X and y being the horizontal distances of node i with respect 
to the eleff»ental area dA of the element. 


Eg. 4,8 is integrated numerically to obtain the influence of 

the normal stress, o' , on the displacements of the N nodes. As in 

J 

the case of shear stresses, the vertical displacement of node i is 
not only affected by the normal stress on element j alone but also 
by its image on element j"". The vertical displacements of node i 
due to stresses on both these elements are written as 


r 
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(4. 10) 
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* _. * . . 
where and p^^^ are displacements 

are respectively the dimensionless 

normal stresses acting on elements 

functions of the parameters, L /B , 

node i and element j. The vertical 

act in the same direction and hence 


of node 

i and, I * 

and 

I ^ 




'' j 

influence 

coefficients 

due 

to 

j and j 

r 

I.* and 

r 

are 



J 


U^/B . V 

O f s 

and the locations 

of 


^4 ^ 

displacements, p and p ^ , 

Z'-J zij 


are additive. 


Combining Egs. 4.9 and 4.10, the vertical displacement at node i 
is obtained as 


*'■1 
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•f 
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(4.11) 


The displacefficnt , , at node i due to stresses on all the N 

elements is 



j = i 


(4. 12) 


The vector of the vertical diplacements of all the nodes along the 
half-length of the strip is 



J *6 i 


Eq. 4.13 is rewritten as 


(4. 13) 



(4. 143 





where the vectors. 



1 

f ) 





}■ and i 

[ O' X, are of size N, 

[■1 

is 

a 


square matrix of size N. 


The solution for the normal stresses are obtained by satisfying 

the compatibility of displacements at all the nodes at the 

soil-strip interface. In section 4.2.1, it was stated that since 

the strip is rigid it undergoes a uniform displacement in the 

vertical direction. Hence, the net vertical displacement of the 

strip is equal to the rigid body displacement, 6 , i.e. 

o 

{ ^ - {'’I } = { 


Combining Egs- 4*5, 4*14 and 4.15^ one gets 






(4. 16) 


Eq. 4.16 gives N equations while there are N+1 unknowns, viz. N 
normal stresses and one rigid body displacement, 6^. The 
additional equation is obtained by satisfying equilibrium at the 
strip-soil interface, i.e. 


N 


^ ^ dA = 0 


(4.17) 


t =5 1 


where ^ F is the summation of all forces acting on one half of 
the strip in the z-direction and dA is area of each element along 



2-axis 


the strip. Since the stresses are symmetrical about the 
the summation is carried out along half the length of the strip 
only. Eqs. 4.16 and 4.17 are rearranged and written as 

[ l" ] I o'/q J + } t4.18) 

^ = 0 193 

Eqs. 4.18 and 4.19 give N+1 equations for N+1 unknowns. They are 
solved by the Guass elimination technique for N values of 
normalised stresses, and the value of the normalised rigid body 
displacement, <5 E /B,q. 

OS f 

Once the normal stresses mobilised at the strip -soil 
interface are computed, the settlement reduction coefficients 
along the surface are calculated on the same lines as in the 
previous chapter. The mechanism for settlement reduction is that, 
since the stresses are acting in the upward direction over the 
middle of the strip they cause the soil above the strip to move 
upward. As a result the surface heaves with a maximum heave at 
the origin, 0. In order to arrive at this heave profile of the 
surface, Mindlin^s solution for vertical displacements due to a 
vertical point load within an elastic half space is used. 
Referring to Fig. 4.5, the vertical displacement of point A, due 
to a vertical load, Q, is given by Eg. 4.6. The vertical 
displacement of point, B (x,0,0), on the surface is arrived at by 
substituting z=0 in Eg. 4.6 which reduces to 



(4.20) 
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where R = x* + 


Eq. 4.20 IS integrated over the area of the reinforcing strip in 
order to evaluate the vertical displacements of points on the 
surface due to normal stresses at a depth, U , as 
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dA (4.21) 
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where R = Vx - * 


U and G = c 

o s 2il+v y 


The vertical displacement of any point k, on the surface due to 
normal stress, a , on element j is given by 
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xkj 
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s e J ^ K 


8(l-i> )^ 4(l-w )U* 

s ^ s o 


.3 




dA (4.22) 


, X, . being the distance along the x-axis 

k j o k j 


where "R = -/ 

between point k, and the centre of element j. Eq. 4.22 is 
rewritten as 


*kj 


B r 
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(4.23) 


Displacement at point k, due to normal stress on the image 
element, j' on the left half of the strip is 



(4.24) 
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uhsre and are dimensionless displacement coefficients 

for the influence of normal stresses on the elements J and J' and 
are functions of the ratios, the location of 

the respective elements- 

Corribining Eqs. 4.23 and 4.24 and summing up for stresses on all 
the elements along the strip, the displacement of point k, is 

written as 
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(4.25) 


j = 1 


The vector of the vertical surface displacements is 




(4.26) 


As mentioned in chapter 3, the vector, 
vector , {^} of size N. The matrix ^ j 


is of size, and 


is of size N^. x N. 


Eq. 4.26 is rewritten as 
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where 

I 

n 

1 is a vector 

of the Settlement Reduction Coefficients, 

def ined 

in 

chapter 

3 as 




(4.28) 
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The vector of the SRC for points on 

the 

surface is 
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(4.29) 

where -f I 1 

1 J 

^ is obtained as 








(4.30) 

where matrix ^ j vector | 

are 

known. 



4.3 Two Strips at Depth, U 

o 

4.3.1 Problem Definition 

Two strips of size, 2L k 2B (Fig. 4.7), are placed symmetric 

r r 

about the x-axis at a depth, U , within the soil mass- Each strip 

o 

is placed at a distance, S , from the x-axis. The problem 

y 

definition is the same as in section 3.3.1 with the only 
difference being that the vertical displacements of the strips are 
considered. 

The normal stresses mobilised along both the strips are 
equal due to the symmetry of placement of the strips. The vertical 
displacements of the nodes along the length of one strip are 
influenced by the stresses acting on the second strip, in addition 
to the influence of the stresses acting on the first strip itself. 



4.3.2 Analysis 


Eq. 4.3 IS used to calculate the vertical displacements 

along the length of the strips due to surface loading with 

f ^ 2 

r=y X + ts^-y> . These displacements being identical for both 
strips, computations are made for one strip alone. The vector of 
the displacements of all N nodes is given by 




(4.31) 


where 


vectors 




are of size N. 


Elements of 



in this case are dependent, in addition to the parameters, U /B , 

of 

L /B,, and v , on the ratio, S /B . 
f f s y f 


The vertical displacements of the nodes due to normal 
stresses acting on the strips are calculated using Eg. 4.6. The 
vertical displacements of node i, due to normal stresses, acting 
on elements j, and its image j' (strip 1), j" and its image j"' 
(strip 2) (Fig. 4.8), are obtained as 


and 



(4.32) 


(4.33) 


(4.34) 


(4.35) 




Rg. 4-8 Discretization of the Strips 



where and are dimensionless coefficients for the 

influence of the noririal stresses on elements 

on the displacements of node i and are functions of the parameters 

•-/Bf' the locations of the node i and 

element j. The vector of the vertical displacements of the N 
nodes due to normal stresses on all the elements is 


The vectors 

and O’ 

^ are 

of size N, 

and 

['1 

is a square 
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matrix of size N 
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Satisfying the compat ib i 1 i ty of displacements at the nodes as in 
Eq. 4.16 one obtains 







(4.37) 


The additional equation of force equilibrium as in Eq. 4.17 is 

N 

^ ^ dA = 0 (4.38) 

i ■ 1 

where the parameters have the same meaning as in Eqs. 4.16 and 
4.17. On rearranging Eqs. 4.37 and 4.38, N+1 equations for N 

unknown stresses, and 1 unknown rigid body displacement, 6^, 
become available as 



(4.39) 






(4.40) 


Solution of these equations yields the N values of the normalised 

normal stresses, a/q, and the normalised rigid body displacement 

5 E /B q. 
os f 


Having calculated the normal stresses, the settlement 
reduction coefficients arc computed based on the formulations 
discussed for a single strip. The vertical displacement of the 
points on the surface due to the normal stresses are calculated 
using Eq. 4.20 with R = -/ + S* + U* . The vertical displacement 
of point k, on the surface due to normal stresses on the elements 


jf 

J'-' 

(Fig. 4.8) is 
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and 
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I ® , and I are displacement 
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coefficients for 

the influence of 

the normal stresses on elements 

j» J'r j"r and 


j ^ respectively and are functions of the parameters, 



location of the nodes. 


Corribining Eqs. 4.41 through 4.44, 
of node k due to normal stresses 


the vertical displacement, 
on all N elements is 
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(4.45) 


The vector of the vertical displacements of points along the 
surface is 





(4.46) 


Eg. 4.46 gives the heave profile of the surface. The SRC are 
obtained on similar lines as in Eqs. 4.27 through 4.30. 


4.4 Two Strips - One Below the Other 


4.4.1 Problem Definition 

As in chapter 3, two strips of size, 2L x2B , placed beneath 

r r 

the centre of the loaded area along the x-axis, one at a depth, 
U , and the other at a depth, U , are considered (Fig. 4.9). The 

O i 

assumption regarding the applicability of Boussinesq^s solution 
remains the same. The normal stresses mobilised along each strip 
are different since they are placed at different depths. In 
calculating the vertical displacements of the nodes along each 
strip due to mobilised normal stresses, the influences of both the 
strips on each other are to be considered. 
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4.4.2 Analysis 

The analysis is similar to the one described for shear 
stresses in chapter 3. The vectors of vertical displacements at N 
nodes along strips 1 and 2 due to surface load, q, as calculated 
from Eq. 4.3 are 


and 


(4.47) 


(4.48) 


where and 


{ } """ { } 


the parameters v and the location of the nodes. 

f f s 


are vectors of size N and depend on 

and 


{^ 1 } 


{'0 are also functions of and respectively. 


The vertical displacements of the nodes along each strip due to 
normal stresses on them are computed using Eq. 4.6. The vertical 
displacements of node i due normal stresses on element j along 
each strip are calculated on the same lines as in section 3.4.2. 
The vector of the displacements at N nodes along strips 1 and 2 
due to the normal stresses along each strip are 
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1 " " 1 are dis- 

placement vectors for strips 1 and 2 due to stresses on strips 1 


and 2 respectively. [ ^ j ' [ ^ j' [ 


21 


and 


[ ' “ ] 

[ I and I 


are 

on 


matrices for the influence of the stresses^ h 

the displacements of nodes along strips 1 and 2 , | } and { 

being the normal stress vectors along strips 1 and 2 respectively. 
All the vectors mentioned in Egs. 4.49 through 4.52 are of size N 
and the matrices are square matrices of size N. 


Combining Eqs. 4.49 through 4.52, the vector of vertical 
displacements of N nodes along strip 1 and N nodes along strip 2 
is 




where 


{ p' } and { } 


are vectors of size 2N, and 


(4.53) 


IS 


square matrix of size 2N. 


consists of vectors 

^ O' ^ consists of the vectors ^ . The matrix ^ j 

comprises of the matrices ^ 


expressed in the same fashion as in Eq. 3.68. 



combining Eqo. 4.47 and 4.48 tha vector of vertical diaplacaments 

of the nodes along strips 1 and 2 due to the surface load, q, is 

given by 




(4.54) 


where 
of vectors 
and 


{ vectors of size 2N. 

^ ^*1 ^ ^xz ^ 1 

{ }• 


comprises 

vectors 


Satisfying the compatibility of displacements at all the 2N nodesp 
one obtains 


{».} ‘ 


(4.55) 


where 6 J- is a vector of size 2N consisting of vectors, 
6 -I 1 I for strip 1 and 6 -I 1 >• for strip 2 each of size, N. 

Oil J °n J 

Substituting Eqs. 4.53 and 4.54 in Eg. 4.55, one obtains 




(4.56) 


Eq. 4.56 gives 2N equations in 2N unknowns for the normalised 
normal stresses, { cy/q and 2 unknowns for the normalised rigid 

body displacements, Hence the 

equilibrium equations to be satisfied are 


two 
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(4.57) 
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(4.58) 


The solution of Eqs. 4.56 through 4.58 gives 2N values of normal 
stresses along strips 1 and 2 and the rigid body displacements of 
strips 1 and 2. 


The SRC due to these stresses are computed on similar lines 
as explained in earlier sections. 

4.5 Multiple Strips 

The problem is similar to the one considered in chapter 3. 

Two pairs of strips are considered, symmetrically placed below the 

loaded area at a depth, U . The first pair is contiguously placed 

o 

(S /B^=B /B^) while each strip of the second is placed at a 

f r f 

distance^ S /B. from the x-axis (Fig. 4.10). 

y2 I 

4.5.1 Analysis 

The analysis^ again is similar to the one for shear stresses 
as explained in chapter 3. Only one strip of each pair is 
considered, due to symmetry. The vector of vertical displacements 
of all the 2N nodes along both the pairs obtained from the 
Boussinesq solution for vertical displacements given by Eq. 4.1 
is 




q 


(4.59) 




wh^re the vectors { p; | and { l' } are of siee 2N. 


The vector of vertical displacements of all the 2N nodes due 
to normal stresses acting on both pairs is obtained on similar 
lines as explained in chapter 3. For this purpose Mindlin^s 
equation for vertical displacements given by Eg. 4.6 is used. In 
computing the vertical displacements due to normal stresses the 
influence of the stresses on all elements including the image 
elements along both the pairs of strips are to be considered. The 
vector of the vertical displacements of the 2N nodes due to normal 
stresses, along the first pair and ^ along the second 

pair of strips is 




(4.60) 


where 
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r I 1 and [ I 1 

{ p } and I 


and I “ I are matrices of the influences of the normal 
stresses, ■! J- and ■{ J-, along both pairs of strips on the 


displacement of the N nodes along the first pair of strips. 

I ** jj and ^ j matrices of the influence of the normal 

stresses along both the pairs of strips on the displacement of the 
N nodes along the second pair of strips. The vector {^} 

comprises of the vectors { ^1 } and { "4- 


Compatibility of displacements at the 2N nodes is satisfied 


as in Eq. 4.55 as 



(4.61> 


{ '’1 ] - { ^ } 

where <5^ ^ is a vector of size 2N consisting of vectors, 

1 ^ the first pair of strips and ^ ^ the second 

pair of strips. Substituting Egs. 4.59 and 4.60 in Eg. 4.61, Eg. 
4.61 is rewritten as 
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Eg. 4.62 gives 2N equations in 2N unknowns for the normalised 


normal stresses 


^ <y/q I 


and 2 unknowns for the normalised rigid 


body displacements, <5 E /B,q and 6 E /B,q. Hence the two egui- 

Ol B f 02 S 1 

librium eguations to be satisfied are 
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= 0 


= 0 


(4.63) 


(4.64) 


The solution of Eqs. 4.62 through 4.64 gives 2N values of 
normal stresses along pairs 1 and 2 and the rigid body 
displacements of pairs 1 and 2. 


The SRC due to these stresses are computed on similar lines 


as explained in earlier sections. 



4.6 Results and Discussion 

The analysis aims to evaluate the normal stresses developed 
at the soil-strip interface, and estimate the settlement reduction 
of the surface thereof. As indicated in the analysis, integrations 

of the basic solutions given by Boussinesg and Mindlin are carried 
out numerically. 

For the purpose of integration, the descretizations of the 

loaded area and of the reinforcing strip are the same as those for 

the case where shear stresses alone were considered (chapter 3). 

That is, the loaded area is divided into subelements of size, 

0. 025B^xO- 025B^ , and the size of elements along the strip is 

O.lBxO.lB for centrally placed strips and O.lBxO.OSB for 
11 f f 

strips placed off the centre. The influence coefficient, l’" , for 

vj 

the vertical displacement of node i, due to normal stress on 
element j, is evaluated by dividing the element j into 20x20 
subare-as for |i-j|:<2. Coarser subdivision of 4x4 was found to be 
adequate for cases |i-j|>2. The half-width of the strips, B , is 

r 

maintained at 0.05B^ for centrally placed strips and 0.025B^ for 
those placed off the centre. 

A parametric study brings out the effects of the depth, 
length and distance of the reinforcement from the x-axis, the 
aspect ratio of the loaded area, and Poisson's ratio of the soil 
on the mobilised normal stresses and the settlement reduction 
coefficients along the surface. The ranges of the parameters 
considered are the same as in chapter 3. They are 



1 to 10 


1) L^/B^ 

2) L /B, 

r t 

3) U /B^ 

o f 

4) S /B, 

y f 

5) V 

s 

6 ) B / B 

r f 


1 to 5 
0.1 to 2 


(aspect ratio of loaded area) 
(length of reinforcement) 
(depth of reinforcement) 


0.025 to 2 (distance of reinforcement from x-axis) 
0.1 to 0.5 (Poisson's Ratio of the soil) 

0.025 for S /B >0 

y ( , . 

(width of reinforcement) 

0-05 for S /B = 0 

y f 


Except for cases where the effect of Poisson's ratio, v , is 

s ^ 

studied or unless mentioned all results are for u = 0.3. 


4.6.1 Normal Stresses - Single Strip 

Initially, the effect of the parameters mentioned, on the 
mobilised normal stresses are studied. The variations of the SRC 
are discussed in the next section. 


Fig. 4.11 presents the variation of the normalised normal 
stress, cr/q, with the distance, x/Bj., along the half-length of the 
strip for various depths, U /B,, below a square area (L/B=l), 

Ox f f 

and strip length of L /B=2. A positive stress is one that acts 

r f 

upwards. It is seen from the figure that at a distance, x/B^, in 
the range 1 to 1.15 the stresses change sign. Further, the 
negative stresses at the edge of the strip increase rapidly, 
tending to infinity. This is due to the edge effect. Also, as 
expected the stresses at the centre of the strip are maximum. 

It is interesting to note that for the depths studied, the 
stresses are maximum for a strip at depth, U^/Bj,=0.25. As the 
depth of the strip increases, the stresses decrease, the magnitude 



At very 


of stresses being small at great depths (U /B =2) 

o f 

shallow depths of placement of strip (U^/E^<0.25) , i.e. very close 
to the loaded area the vertical displacements due to the surface 

load are high. At this depth, the interference effect of the 

normal stresses on the displacements of points along the strip are 

also high. Hence, the stresses required to be mobilised so that 

the strip displaces uniformly is small. With increase in depth of 
placement of strip, U^/B^, upto 0.25, the interference effect 
decreases although the displacements due to the surface load are 
high. The stresses mobilised, therefore are higher. With further 
increase in U^/B^, the displacements decrease and hence stresses 
are small and tend to be uniform. The normalised normal stresses, 
©■/q, at the centre for strips at depth, U /B=0.25 and 2 are 4.25 

O f 

and 1.3 respectively. 

The effect of the length of the strip, L /B,, on the 

r f 

normalised normal stresses, cr/q, for strips placed below a square 
area with L^/B^=l, and at a depth, U^/Bj.=l, is presented in Fig. 
4-12. It is observed that for all lengths of the strip, the 
stresses decrease with the distance along the strip, x/B^, and 
become negative beyond a certain distance. This negative stress 
increases with x/B^ and increases very rapidly at the edge. For 
strips with L /B<3, positive stresses are developed over a 

r 1 

distance of 0.6 times the half length, L . For longer strips 

r 

(L/B>3), stresses are negative over predominant lengths of the 

r f 

strips. The normalised stress at the centre of the strip, 
(x/B=0> for L /B =1 and 5 are 1.3 and 5.5 respectively- In the 

f r f 

case of shorter strips the difference between the vertical 


Normal Stress g/q 



Fig. 4.11 Variation of Normal Stress with Distance for 


Different Depths of Placement of Strip 
(U/B,=2, Lf/B,= 1, t's=0.3) 



Fig. 4.12 Variation of Normal Stress with Distance 
for Different Lengths of Strip 
(Uo/8,= 1. U/B,=l i^s=0.3) 


displacements of the points at th«. ^ . 

Kwjiiiis ax trie centre and the edge of the 

.trip due to surface load is soall. As a result, stresses of 

lesser magnitude are sufficient to force the strip to undergo a 

uniforii, rigid body displacement. On the other hand, the gradient 

betoeen the vertical displacement due to surface load at the 

centre and edge of longer strips is high. Conseguently. higher 

stresses are required to make the strip displace uniformly 

throughout its length. The stresses change sign at a distance 

approximately equal to 0.45 to 0.55 times the half length, L , of 

r 

the strip. 

The aspect ratio of the loaded area, appears to have 

a considerable influence on the mobilised stresses, a/^ (Fig. 
4.131. The figure depicts the variation of a/q, with the distance 
along the strip, x/B^, for various aspect ratios, 

strips of length, L /B =2, placed at a depth, U /B=l. For all 

r 1 Of 

L^/B^ values, the variation of stresses is similar. The stresses 
change sign at a distance, x/B^, of 1.1. Longer rectangles with 
L^/B^=10 generate maximum stresses at the soi 1-reinformcement 
interface, the stress, o', near the centre being 5q. The increase 
in stresses from square (L^/B^=l) to rectangular areas with 
L^/B^=2 is significant. With further increase in the aspect ratio 
of the loaded area (L,/B.) the increase in the interaction 
stresses is small. The reason for this is that as the aspect 
ratio increases, the loading at points farther away from the 
x-axis have an insignificant effect on the displacement of points 


along the strip. 



% The plot of the normalised normal stress, C'/^, with 
shown in Fig* 4-14, for different values of Poisson'^s ratio, 
and for a strip of length, L /B =2, placed below a square area 

r I 

(L^/Bj.=l) at a depth, U^/Bj=l. It is clear from the plot that 

Poisson ratio does not affect the stresses mobilised along the 

strip appreciably. For v =0.5, the stresses are marginally higher 

s 

than the stresses for other values of v . The maximum cr/q, 

s 

observed at the centre is 3.4q for v =0.5. 

s 

From these plots of the normal stress variations it is 
evident that strips placed at a' depth, U^/B^, of 0.25 develop 
maximum stresses and stresses are higher for rectangles with 
larger aspect ratios. The stresses are not sensitive to the value 
of the Poisson's ratio. 

The normal stresses developed at the interface are 
responsible for the settlement reduction of the surface. In the 
next section the effect of various parameters defined earlier on 
the settlement reduction of the surface will be discussed. 

4.6.2 Settlement Reduction - Single Strip 

Before studying the reduction of settlements at the surface 
due to the normal stresses, it would be worthwhile to examine the 
basic solution given by Mindlin for vertical displacements due to 
a verical force at depth. The vertical displacement of point, 
B(x,0,0) (Fig. 4.5), due to a vertical force, Q, at a depth, c, xs 
given by Eg. 4.20. Non-dimensional ising all length parameters 

with c, Eq. 4.20 is written as 
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(4.65) 


where G^=E^/2( 1+v^) and is an influence coefficient dependent 

on the normalised distance, x/c and Poisson^s ratio, . 

s 

The variation of with the normalised distance, x/c, for 
different values of Poisson's ratio, v is presented in Fig. 4.15. 
For all values it is observed that is maximum for x/c=0 and 
it decreases rapidly with increase in the distance, x/c. I values 

Z 

are positive for all values of i> indicating surface heave 
throughout. values asymtotically approach zero at large x/c 

values. This figure shows that the points on the surface just 
above the point force experience the maximum vertical displacement 
in the direction of the force and the points further away from the 
force displace by smaller amounts. It is apparent that u =0.1 

S 

shows the maximum I values, and =0.5 the least. This difference 

z s 

arises since the parameter, G , has been used in Eq. 4.65. If the 

s 

modulus, E , is substituted for 6 in Eq. 4.65, I values will be 

3B S z 

less affected by v values. It is also perceived that beyond a 

s 

distance, x/c, of about 3, the I values are very small as 

Z 

compared to the maximum value at x/c=0. This indicates that the 

zone of influence for the force at a depth, c in producing heave 

at the surface is of the order of 3c. The values of I^ at x/c=0 

for V =0.1 and 0.5, are 0.223 and 0.158 respectively. At a 
s 

distance of x/c=3, the corresponding values of I^ are 0.05 and 
0.03 respectively. 


The effect of length of the strip, L^/B^, placed below a 



square area CL,/B=1) at depth ii /p -i _ n. 

f j -i oepxn, u^/B^-1, on the SRC is presented 

in Fig. 4.16. The SRC values are maximum at the centre of the 

loaded area (x^/B^=0) and gradually decrease with increase in the 

distance, for all lengths of the strip. The edge of the 

loaded area shows minimum settlement reduction. It is also noted 

that SRC values increase with increasing lengths of the strip. 

The SRC at the centre for L /B =1 is 0.007 and for L /B =5 is 

r 1 r f 

0-13- Strips with L^/B^=l show a small negative SRC at the edge 

indicating settlement of the surface. The increase in SRC with 

L^/B^ is due to the fact that the stresses for longer strips are 

higher than for shorter ones as seen from Fig. 4.12. Further from 

Fig. 4.15 it is clear that points on the surface within the 

distance, x/c of 2 show significant heave. Beyond this distance, 

the I values are comparatively small. Hence it follows that high 

positive stresses within a distance, x/Bj.=l, cause high settlement 

reduction of the surface and this is true for longer strips (Fig. 

4.12). The decrease in SRC with increase in x./B. is obvious 

because the stresses are a maximum at the centre of the strip i.e. 

at x/B=0 and decrease with x/B,. Shorter strips (L/B=l) show 

negative SRC at the edge because the stresses beyond a distance, 

x/Bj.=0.5 are negative which cause the edge of the loaded area, 

X /B=l, to settle, 
f f 

Fig. 4.17 shows the variation of SRC with the distance, 

X /B , along the surface for various depths of placement of 
f f 

reinforcement, U^/B^, and for a length, L^/B^=2, placed below a 
rectangular area with aspect ratio, L^/Bj,=2. All the curves start 
off from a maximum value at x^/B^=0 and SRC decreases with 



Rg. 4.15 Influence Coefficients for Vertical Surface Displacements 
due to Vertical Force at Depth (Mindlin Solution) 


Rg. 4.16 Voriation of SRC with Distance along 
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distance from the centre. The rate of change of SRC with 

maximum for shallow depths and it decreases with increase in 

U /B . At large depths (U /B>1.5) the SRC curve is almost 
of Of 

parallel to the axis indicating uniform heave of the 

surface. SRC values decrease with increase in the depth of 
placement of the strip. It is interesting to note that although 
normal stresses, o/q, for a strip at depth, U^/B^=0.25, are 
marginally higher than those for U^/B^=0.1 (Fig. 4.11), the SRC 
values for the latter are higher. This is because the latter is 
closer to the surface and the stresses have a greater influence on 
the settlement reduction of the points on the surface. As the 
depth of placement of the strip increases, the influence of the 
mobilised stresses in reducing surface settlements decrease 
rapidly. 

The variations in SRC for different aspect ratios, of 

the loaded area arc presented in Figs. 4.18 and 4.19. The 
variation of SRC with the distance, x/B^, along the surface for 
strips of length, L /B =2, placed at a depth, is depicted 

in Fig. 4.18. The maximum settlement reduction is for long 
rectangular areas (L/B =10) while square areas exhibit the least 
SRC values. It is noted that the increase in SRC is significant 
for an increase in the aspect ratio from 1 to 2 

further increase in L/B, values upto 10, the increase in SRC is 
not as Significant. The SRC at the centre of the loaded area with 

an aspect ratio, L/B,=10, is about 0.07 while for a square area 

. .. _i i e. a+ V /R =1 the values of SRC are 

it is 0.044. At the edge, i.e., at 

« - ^ i.r I /R Kiaiatively less significant, 

small and the variation with relatively 



SRC 


Uo/Bf 




In Fig. 4.19 the variation of BRC with distance, for 

strips at depth, U^/B^=2, is depicted. The trends in this plot 
are similar to those observed in Fig. 4.18. The SRC values are 
much smaller than the corresponding values of SRC for strips at 
depth, U^/B^=l, as seen in the previous figure. It is seen that 
the increase in SRC for L^/B^ increasing from 2 to 10 is as 
significant as its increase from L^/B^=l to 2. The SRC values at 
the centre, i.e. x^/B^=0, for square, and rectangular areas with 
L^/B^=2 and L^/B^=10, are 0.006, 0.0096 and 0.013 respectively. 

From Figs. 4.18 and 4.19 it is evident that for all depths, 
long rectangular areas (L^/B^=10) show maximum BRC values. With 
the decrease in the aspect ratio of the loaded area the SRC values 
decrease, with square areas exhibiting the least SRC values. 

Thus far the variation of SRC along the surface due to the 

effect of the various parameters was studied. It is observed from 

the figures that the centre of the loaded area experiences maximum 

settlement reduction. It would also be of practical interest to 

consider the settlement reduction at the centre of the loaded 

area. The effect of the various parameters on the SRC at the 

centre, defined as I , is presented in the ‘succeeding graphs. 

sc 

Fig. 4.20 plots I values with the depth ratio of the 

sc 

strip, U /B , for various values of the aspect ratio, L /B , of 
of * ‘ 

the loaded area and for a strip of length, L^/B^=2. I^^ values 
decrease with increase in the depth of reinforcement, U^/B^, 
all L^/B^ values. They increase with the aspect ratio, 


of 









the loaded area, with long rectangular areas (Lj./Bj.= 10) exhibiting 
highest values. The increase in I for areas with LVB, 

SC f f 

increasing from 2 to 10 is not as significant as the increase for 

areas with L^/B^ increasing from 1 to 2. At depths, U /B>1.5, 

o f 

the value of are negligible. The reason for the rapid 

decrease in the values of with the depth of reinforcement is 

that the magnitude of stresses decreases with depth. In addition 

to this, with increasing depths, the stresses mobilised have a 

decreasing influence on the displacements of points along the 

surface. At a depth, U /B=0.1, the I values for aspect ratios, 

and 10 are 0.245 and 0.31 respectively. At a depth of 

placement, U =B,, the difference in I values for L /B =1 and 10 
of sc f f 

is much smaller, the corresponding values of I being 0.04 and 

SC 

0.07. This plot implies that rigid strips placed at shallow 
depths are advantageous in reducing settlements. Also, strips 
placed below longer rectangles give better results. 

The effect of the aspect ratio, l-j./B^, of the loaded area on 

I for different lengths of the strip, L /B,, is depicted in 
sc r t 

Figs. 4.21 and 4.22. Fig. 4.21 shows the plot for the depth, 

U /B =1. I values increase with the length of the strip, and 
of sc 

the aspect ratio of the loaded area. For shorter strips CL^/B^=11 

the I values are very small and the aspect ratio, L /B r has a 
sc * ’ 

minimal effect on these values. For longer strips, e.g. L^/B^=5, 

the aspect ratio, L^/B^, affects the values considerably. The 

1 values for this length of the strip (L /B =51 are 0.12 and 
sc ’ 

0.25 for aspect ratios, L/B-l and 10 respectively. 
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The Plot for a deeper strip, i.e. for U^/B^=2, is presented 
in Fig. 4.22. Here again, the values increase with the length 
of the strip, L^/B^, and the aspect ratio, of the loaded 

area. There is a continuous increase in I values with L /B . 

ec r f 

However, they are small as compared to the values for the strips 
placed at a depth, U /B=1.0. The values of I for strips with 

* SC 

are of the order of 0.04 for L^/Bj=l and 0-105 for 
L^/B^=10. 

From these two figures it is evident that longer strips 

below longer rectangles at shallow depths result in significant 

settlement reduction at the centre of the loaded area. Strips at 

shallow depths are more advantageous as compared to strips at 

g»"eater depths. As seen from Fig. 4.12, the stresses arc higher 

for longer strips. Also, from Fig. 4.15, it is clear that the 

maximum heave of the surface is experienced within a distance, x/c 

of 2 from the point force at a depth c. Hence positive (upward) 

normal stresses for strips at depth, U^/B^=l, within a distance, 

x=2B^ (c=U^=B^) and at depth, U^/B^=2, within a distance, x=4B^ 

(c=U =2B ) arc responsible for settlement reduction of the 
o f 

surface. Normal stresses are positive over larger distances for 
longer strips. Hence longer strips result in higher values. 

The effect of Poisson^s ratio, v^r on for strips at a 

depth, U /B=l, placed below a square area (L /B =1) is shown in 
Of * * 

Fig. 4.23. The I values increase with length, For all 

ssc 

values. The effect of v is felt only in case of strips with 

s 

length, L/B^>2. For L/Bj=5, values for p^= 0. 1 and 0.5 

(saturated undrained soils) are 0.16 and 0.11 respectively. 



From Figs. 4.20 to 4.23, it is evident that the settlement 
reduction coefficient at the centre of the loaded area, I , is 

SC 

affected considerably by the length of the reinforcement, L /B , 

r f 

depth of reinforcement, the aspect ratio of the loaded 

area, L^/B^ and the Poisson^s ratio of the soil, v • Hence, in 

s 

designing the reinforced foundation bed, an optimum choice of all 
these parameters is to be made. 

The effect of the length of the strip, L /B,, on the 

r 1 

normalised rigid body displacement of the reinforcement strip, 

<5 E /B.q, for different depths of placement, U /B,, and for a 
square loaded area (L^/B^=l), is depicted in Fig. 4.24. 
decreases with depth, U /B,. It is also seen that 6 E /B q 
decreases with increase in the length of the strip, L /B . The 

r 1 

rate of change of rigid body displacement of strips, with depth, 

U /B , also decreases with the increase in L /B . For a strip of 
of r t 

length, L^/Bj,=l, the values of the rigid body displacement of the 

strip, <5 E /Bq at depths, U /B=0.1 and 2 are 1.78 and 0.85 
o s f of 

respectively. The corresponding values for longer strips with 

L /B =5 are 0.72 and 0.5 respectively. At shallow depths CU^/Bj.<l) 

r f 

the chsingc in the rigid body displacement of the strip for lengths 
of strip, L /B , increasing from 1 to 2, is much more significant 
than with further increase in A® explained earlier, the 

difference between the vertical displacements due to surface 
loading at the centre and edge of the strip is high for longer 
strips and low for shorter strips. The average value of 
displacements for a shorter strip is therefore higher than that 
for a longer one. The uniform displacement value which a rigid 



E,/Bf.q 
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strip maintains lies in the vicinity of this average value. This 
explains why the rigid body displacement, <5 E /B q, is maximum for 

OS { 

strips with lengths, L /B =1 and minimum for L /B =5. Further, as 
the depths of placement of strip^ increases^ the vertical 

displacements of points are smaller and consequently, the rigid 
body displacement decreases. 

Fig. 4.25 presents the plot of the rigid body displacement, 

of the strip versus the length, of the strip for 

various aspect ratios, L^/B^r of the loaded area. This figure 

reiterates the observation that S E /B^q decrease with length, 

o s f ^ J 

L /B,, made with respect to the data presented in Fig. 4.24. This 

is true for all values. It is also observed from this figure 

that the rigid body displacement increases with L,/B,. «5 E /B,q, 

f f o s f 

for L /B,=l is 3.8 for a long rectangular area with aspect ratio, 

r f 


L /B =10 
f f 

and 1 

.23 for a square 

area 

(L /B=l). 
t f 

The increase 

in 

<5 E /B a 
os f 

with 

L^/B^ is because 

higher 

stresses 

are mobilised 

for 


strips below longer rectangles due to high vertical displacements. 


4.6.3 Two Strips at Depth, 

Fig. 4.26 presents the effect of distance of the strip from 

x-axis, S /B , on the normalised normal stresses, cr/^, for a strip 
y f 

of length, L /B =2, placed below a square area (L /B =1) at a 
r 1 

depth, U /B=l. For all S /B, values the stresses are maximum at 
Of y * 

the centre i.e. at x/B^==0, decrease with the distance, x/B^ and 
increase rapidly on the negative side towards the edge of the 
strip. The stresses change sign at a distance, x/B^ of around 
1.1. Normal stresses are maximum for S^/B^=0.5. At larger 
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*9* 4-25 of Length on Rigid Body Displacement of 

ri I 'In Aspect Ratios of Loaded Area 

(Uo/Bf=1, ^s=0.3) 


Fiq 4 26 Variation of Normal Stress with Distance for 
Different Distances of Strip from x-axis 
(U/B,=2. L,/B,= 1. Uo/B,= 1. Vs=0,3) 





are very small. The stress 


distances the stresses 

<r/q, for tWD contiguous strips, i.e. S^/B, =0.025, at tha centre 

(x/B=0) is 3.3. The stress at the centre for S /B=0.5 and 2 are 

y { 

4.6q and 1.2 q. At closer spacings of the strips, that is for 

distances, S^/B^<0.5 the interference of the normal stresses over 

each other is high and as this distance increases upto S /B=0.5 

y { 

the interference effect decreases. As a result stresses increase 

upto S^/B^=0.5. For S^/B^>0.5 the vertical displacements of the 

points along the strip due to the surface loading are themselves 

small. Consequently the stresses decrease with S /B . 

y t 

The plot of SRC at the centre, i.e I , versus the distance 

SC 

of the strip from x-axis, S /B^, for various lengths, L /B , of 

y f r f 

the strip placed at a depth, U /B=l below a square loaded area is 

of 

shown in Fig. 4.27. The I for all lengths of strips, increases 

SC 

upto a certain value and then decreases with S /B,. The maximum 

y f 

I is observed at a distance, S /B. of around 0.15 to 0.25, 
oc y f 

depending on the length of the strip. Further, the I values 

SC 

increase with L /B, and this is obvious because the stresses 

mobilised are higher for longer strips (Fig. 4.12). The maximum 

values of I of 0.005 and 0.081 for L /B=l and 5 are observed at 

sc r 1 

distances, S /B =0. 15 and 0.25. At S /B =2 the I values for the 
y I y * 

corresponding lengths are practically zero and 0.01 respectively. 

The maximum I is observed at a distance, S /B of around 0.25 
sc y I 

because the stresses developed are maximum at this distance as 
perceived from Fig. 4.26. From this figure it is clear that for 
rigid strips at a depth, equal to B^., maximum settlement 

reduction accrues if the two strips are placed each at a distance 


of 0.25B, from tho x-axis. Further, longer strips result in higher 
settlement reduction. 

Fig. 4.28 depicts the variation of I with depth ratio, 
U /Bf, for various values of the distance, S /B., of the strips of 

y f 

length, L^/B=2r from the x-axis placed below a square area. It 

is seen that beyond a depth of 0.5B^ the I values decrease for 

all values. At shallow depths tU^/B^<0.25) there is a 

change in the slope of the I curve with increase in the values 

of from 0.025 to 0.25. For S /B <0.25, the maximum I is 

y * y f sc 

observed at a depth of 0. IB. whereas for S /B>0.25 it occurs at a 

f y f 

depth of 0.25B^. Further, in the depth range of 0.2B^ to 0.6B^, 
I is maximum for strips placed at a distance S /B =0. 1 and for 

SC y f 

depths, U /B,>0.6, it is maximum for S /B=0.25. This shows that 
o t y i 

the distance, S /B at which maximum settlement reduction occurs 
y f 

is a function of depth of the strips. At a depth, U /B=0. 1, 

o f 

contiguous strips CS /B =B /B =0.025) exhibit maximum I . At 

y f r f sc 

greater depths CU /B,>1.5), I values are very small and the 

o f sc 

effect of S /B, is insignificant, 
y f 

Figs. 4.29 and 4.30 show the effect of the parameters on the 
normalised rigid body displacement, 6 E /B q. The variation of 

O o 1 

6 E /Ba versus the distance, S /B, of the strips from the x-axis 
o s y f 

for various lengths of the strip, L /B , placed at a depth, 

r 1 

U /B =1 below a square area is presented in Fig- 4-29. The rigid 
o f 

body displacement decreases with for all lengths. The rate 

of change of 6 E /B,q is maximum for shorter strips and it 

OS! 1 

decreases with increase in the length of strip. 


The maximum value 












“^served for contigoous strips. It is also evidant 

from tha figura that tha rigid body displacamants dacraasa with 

increase in length, of the strip, as observed in Figs. 4.24 

and 4.25. For contiguous strips the value of 5 E /B q are 1.23 

o s 

and 0.59 for L/Bj=l and 5. At S^/B,=2 the corresponding values 
are 0.6 and 0.38 respectively. The decrease in <5 E /B q with the 

OS f 

distance, is obvious because the vertical displacement of 

the points along the strip decrease with S /B 

y f' 

Fig. 4.30 shows the plot of <5 E /b a with the depth ratio, 

U^/B^, for various distances from x-axis, S^/B^, and for strips of 

length, L^/B^=2, placed below a square area. The rigid body 

displacement decreases with depth for. all S /B values. With 

y f 

increase in S /E , 6 E /B.q decreases. The rate of change of the 
y 1 o s 1 

rigid body displacement is high for smaller distances, S /B , from 

y f 

the x-axis. The decrease in <5 E /B,q with S /B, is due to the 

o s f y » 

fact that the vertical displacements due to surface loading 

decrease with increasing S /B, values. The decrease in rigid body 

y f 

displacement with depth is again due to the decrease in vertical 

displacements with U /B, values. 

of 

4.6.4 Two Strips One Below The Other 

The interference of two strips placed one below the other is 
investigated next. The stresses are calculated and the 
corr espond ing effect on the settlement reduction of the centre of 
the loaded area on the surface studied. 


The comparison for normal stresses for strips one below the 
other at depths, U^/B^=0.75 and with those for single 



influence on each other. As a result, the stresses mobilised are 
smaller as compared to those mobilised if considered indepen- 
dently. Further, this difference in stresses is higher for longer 
strips. The consequence is that the settlement reduction reduces 
as compared to strips considered independently. It is also evident 

that even for a vertical spacing of Au=1.75 the value of I for 

1 > 

1 and 5 are 0.72 and 0.66. From this table it is clear that 
the principle of superposition cannot be applied for two strips 
placed one below the other at least for the depth range considered 
in the analysis. A suitable correction factor will have to be 
used to calculate the settlement reduction coefficients due to two 
strips acting together from those of strips acting independently. 

4.6.5 Multiple Strips 

Two pairs of strips are considered at the same depth and the 
influence of each strip on the other is investigated. One pair is 
contiguously placed (S^^/B^=0.025) while the distance of the 
second pair, S /B from the x-axis is varied. The ratio of the 

yz t 

1 for the multiple strips to that obtained from the summation of 

SC 

the I for a single pair is calculated. The effects of the 

SC 

lengths of the strips, L /B, and depth of placement, U /B on the 

r f o I 

ratio, I (defined subsequently), is studied in this section. 

1 > 

The comparison of normal stresses for multiple strips at 
distances, S /B =0.025 (contiguous strips) and S /B =0.5 with 

' yi f y ’ 

those for a single pair of strips at the respective distances is 
presented in Fig. 4.32. Strips considered are of length, L^/B^=2, 
Placed at a depth, U^/Bj=l, below a square area. 


It is seen that 
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the normal Btreasos for the multiple stripB are emaller than those 

tor a single pair considered at the respective distances IS /B>. 

y t 

This again is due to the influence of each pair of strips on the 
other. The vertical displacements along each pair to be 
counteracted by the normal stresses are reduced and as a result 
the normal stresses mobilised are smaller. 


In this case the ratio, I , is defined as 

D 


^or two pairs of strips at S and B 

j ^ y 1 yZ 

» t _ for strips at S + I for strips at S 

SU yl SC y2 


(4.66) 


Fig. 4.33 presents the variation of the ratio, I , with the 

D 

horizontal spacing, S (=S /B,-S /BJ for various lengths of 

yl> yl f y2 f 

the strips, L /B^ placed at a depth, U /B=l, below a square area. 

r f Of 

increases with increase in the spacing for all lengths of the 
strips- For shorter strips the increase in I with S is rapid 

D yD 

as compared to longer ones. At a spacing, S =0.4, I_=0.9 for 

yD D 

strips of length, L /B =1 while for longer strips (L /B =5) a high 

r f r 1 

value of I is observed at large values of S . This is because, 
j> yi> 

at close spacings the influence of each strip on the other is high 
and this is particularly high for longer strips. The stresses 
mobilised are consequently low, resulting in smaller settlement 
reduction coefficients as compared to a single pair considered in 
an earlier section. 


The variation of with for various depths, 

Placement of the strips of length, L/B=2, below a square area is 
presented in Fig. 4.34. It is interesting to note that at closer 
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Fig, 4.33 Effect of Spacing Difference on 
!n for Different Lengths of strip 
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Epacings, for greater depths is higher while this trend 
reverses with inc-ease in the sparing, The reason tor this 
is that at shallow depths, strips placed at closer spacings have a 
higher influence on each other as compared to distantly spaced 
strips. On the other hand, the effect of distantly spaced strips 
becomes predominant at greater depths, as seen in Fig. 4.28. 

From these two figures it is clear that I is a function of 

D 

the length, depth, U^/B^ of the strips and the spacing 
between them. Based on these parameters a suitable factor is to 
be selected from these figures to obtain settlement reduction due 
to multiple strips from those obtained by summing up the 
coefficients for individual pairs. 

4.7 Conclusions 

An analysis using the elastic continuum theory for a rigid 
strip displacing uniformly in soil is proposed. The normal stress 
interactions are studied and evaluated. The effects of the 
parameters mentioned in section 4.6 viz., the aspect ratio of the 
loaded area, length, depth and distance of the strip from the 
x-axis on the normal stresses mobilised at the strip-soil 
interface and the reduction in surface settlements, resulting 
thereof have been studied. It is seen that, in forcing the strip 
to undergo a uniform vertical deformation normal stresses are 
mobilised at the strip-soil interface which act upward at the 
centre and downward at the strip edges. These stresses, in turn 
push the soil up, bringing about a reduction in surface 
settlements. The results of the parametric study indicate a 
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number of facts which are listed below. 

1. A rigid strip placed close to the surface (U /B<0.1) gives 

of 

rise to maximum reduction of the settlement of the surface. 

2. For any depth of placement, for shorter strips (L/B<31, the 

r f 

increase in settlement reduction is appreciable with increase 
in the aspect ratio of the loaded area upto 2 while it is 
marginal for further increase in the aspect ratio. 

3. For any depth of placement, longer strips show higher 
sett lenient reduction of the surface since the mobilised normal 
stresses are positive over larger lengths. 

4. Similar to the results of chapter 3, for a pair of rigid 
strips placed symmetrically about the x-axis and at a depth 
U >0.5B, maximum reduction in surface settlement is observed 

O 1 

if the centre of each strip lies at a distance of 0.256^, from 
the x-axis. 

5. A suitable factor is to be used for obtaining the settlement 
reduction due to two strips placed one below the other and for 
multiple strips from those obtained from the principle of 
superposit ion. 

6. For a rigid strip, the settlement reduction coefficients due 
to normal stress interactions are much higher than those due 

to shear stress interactions. 



CHAPTER 5 

STRIP-SOIL INTERACTION : EFFECTS OF EXTENSIBILITY, FLEXIBILITY, 
COUPLED SHEAR AND NORMAL STRESS INTERACTIONS AND SLIP 

5.1 Introduction 

In the preceding chapters the strip was considered to be 
rigid such that no axial elongation was permitted (Chapter 3) or 
it was forced to undergo a uniform vertical displacement (Chapter 
4). The shear and normal stresses mobilised thereof, considered 
independently, were calculated and the consequent reduction in 
surface settlements computed. In this chapter, firstly the 
extensibility and the flexural stiffness of the strip are taken 
into account separately. Unlike the case of an inextensible one 
(Chapter 3), the horizontal displacements of the points along the 
extensible strip are equated to the elongations of these points 
due to tensions in the strip which are an outcome of the shear 
stresses mobilised at the interface. A flexible strip is assumed 
to behave like a beam and hence the flexural deformations of 
points along the strip are equated to the soil displacements. The 
shear and normal stresses mobilised thereof, were computed as also 
the resulting settlement reductions of points along the surface- 
The combined effect of the horizontal and vertical rigidity for a 
rigid strip is then considered wherein the inextensibility and the 
uniform vertical translation of the strip are incorporated 
simultaneously. An allcjwance is made for possible slip along the 
strip~soil interface by limiting the shear stresses mobilised to 
the frictional resistance available at the interface due to 
mobilsed normal stresses and the overburden stress at that depth. 



5.2 Extensible Strip - Single Strip 


5.2.1 Problem Definition 

The definition of the problem is similar to the one in 
chapter 3 with the difference that the strip (of size, 2L x2B 

r r 

placed centrally beneath a loaded area of size, 2L x2B 

f f 

transmitting a uniform load, q, at a depth, U ) is extensible 

o 

(Fig. 5.1). The width of this strip, 2B , as explained earlier is 
relatively small and its thickness, t^, negligible. The strip 
being extensible, the net soil lateral displacements resulting 
from the applied surface load and the mobilised shear stresses are 
not zero but equal to the axial elongations of the strip. 
Satisfying the compatibility of displacements for points along the 
strip, the shear stresses iMsbilised at the strip-soil interface 
are computed. 


The tension in the strip and settlement reduction 
coefficients of points along the surface are calculated as before. 
Again, the shear stresses developed are assumed to be constant 
over the width of the strip and to vary along its length. 

5.2.2 Analysis 

The discretization of the loaded area and the strip for 
numerical integrations are the same as in chapter 3. The 
Boussinesq'^s and Mindlin'^s equations are used in computing the 
lateral displacements due to the applied load on the surface and 
the mobilised shear stresses respectively- 

Thc vector of horizontal displacements of the N nodes along 
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thp half-length of the etrip d„e to the surface load, g as 
obtained by integrating Boussinesg's eguation over the loaded 


area, is 


where 


{ vectors of size N. 


t5.1) 


The vector of horizontal displacements of the N nodes along 
the half "length of the strip due to shear stresses nobilised along 
the whole length of the strip as obtained from integration of 
Kindling's equation for each element is 




(5.21 


where 




is an influence coefficient matrix for the influence 


of the shear stresses on the horizontal displacements of the N 
nodes. It is a square matrix of size N, and ^ {^} are 

vectors of size N, being the vector of the mobilised shear 


stresses. 


being the vector of the mobilised shear 


The net lateral displacement of all the nodes, as mentioned 
earlier is equal to the elongations of the strip. Thus, the 
compatibility equation is 




(5.3) 



nodes due to the 


where | j is a vector of displacements of the 

elongation of each element of the strip and is of size N. The 
elongations are due to the tension in the strip resulting from 
mobilised shear stresses. The force equilibrium of the element i 
in Fig. 5.2 is vnritten as 


k k k k 

or 

AT = -T dA = -T dl 2B 

k k k r 


(5.4) 


(5.5) 


THe nEgat ivE sign indicatss that ths tfinsion in thE strip 

decreases with distance along the x-axis. The tension, T , to the 

k 

left of elenient k, would then be 


N 


N 


i *s k i * k 


dl 2B 


(5.6) 


The elentental elongation, hp , of the element is 

£K 




Ek 


T +(T +AT ) dl 
k k k 

A~E 2“ 

r r 


(5.7) 


V *' r r 


(5.8) 


, T, ^ . dl 2B 


(5.9) 


i =k-t- 1 


A is the cross-sectional area of each element (assumed constant) 

r 

and E is the modulus of deformation of 


the strip. 


The 


Definition Sketch — Extensible Strip 




Fig. 5.2 q) Stresses on the Strip 
b) Tension in an Element 


1 /d 


cross -sectional area, A^, of the strip is 


A «= t 2B 

r r r 


(5.101 


t^ being the thickness of the strip. Substituting Eg. 5.10 in Eg. 
5.9 and rewriting 




Ek 




(5.11) 


The total elongation of node k is egual to the sum of all the 
elemental elongations from nodes 1 to k, and is written as 


Ek 


= P, 


Ek-1 




Ek 


(5. 12) 


Combining Egs. 5.11 and 5.12, the vector of the elongations of all 
the N nodes is written as 


where ^ J is a sguare elongation coefficient matrix of size, N, 
and is given by 

0.5 1.0 1.0 1.0 1.0 

0.5 1.5 2.0 2.0 2.0 

0.5 1.5 2.5 3.0 3.0 

0.5 1.5 2.5 3.5 N-0.5 



Non-dimensiDnalising the length parameters with B,, Eg. 5.13 U 

rewritten as 



Tt~7B7TT 

r f r 




Substituting for 




in Eq. 


5.3 


(5.14) 


or 




['■){'} 


rt~"7B ) E 

r f r 


1 '.](') 


(5.15) 




where C = (dl/B^)* and is the dimensionless 

s“f 

elongation factor which takes into account the axial stiffness of 
the strip and the modulus of deformation of the soil. 


Eq. 5.16 gives N equations for the normalised shear stresses 
mobilised along the N elements of the strip, which are solved by 
the Gauss elimination technique. 

The tension in the reinforcement and the SRC are calculated 
as explained in chapter 3 for a single rigid strip. 

5.3 Extensible Strip - Two Strips at Depth 

5.3.1 Problem Definition 

The problem definition remains the same as in chapter 3. 



The horizontal displaeement vector for the dieplacemente of all 

the N nodes of one strip due to surface load as expressed br Eq. 
3.40 in chapter 3 is 




(5. 17 ) 


The vector of horizontal displacements due to shear strt 
as expressed by Eq. 3.46 is 




C5. 18) 


where all the vectors and matrix 


as in Eqs. 3.40 and 3.46. 


[■1 


have the same definitions 


The elongation vector 




being a function of the 


stresses and the material properties of the strip, is the same as 
in Eq. 5.14. Satisfying the ccmipatibility of displacements at all 
the nodes one obtains 




(5. 19) 


<5.20) 


the normalised shear stresses 


, { T/q ). 


are obtained by solving 


the N equations. The difference in Eq. 5.16 and Eq. 5.20 is that 
the matriv Ft*'! and vector ( T in Eq. 5.20 are functions of 


z[l'] 


and vector 


{■'}■ 


S^/B^ also, while they were functions of 



1/6 


only in Eq. 5-16. 

The tension and SRC are - 

**re calculated as explained in chapter 3 

for the same problem. 

5.4 Flexible Strips - Single Strip 

5.4.1 Problem Definition 

The strip considered is assumed to deform in the vertical 
direction, i-e., it is assumed to behave like a beam. As in 
chapter 4, a strip of size, 2L^x2B^, is placed centrally at a 
depth, U^, below a rectangular area on the surface transmitting a 
load of intensity, q (Fig. 5-3). 

The surface load causes the strip to deform in the vertical 
direction as shown in Fig. 5.4. Due to this deformation, normal 
stresses are induced in the soil. The total vertical displaceoMent 
of the soil is equal to the displacement due to the loading plus 
the displacement due to normal stresses generated. Satisfying the 
compatibility of the net displacements of the soil adjacent to the 
strip with those of the strip assumed to behave like a beam, the 
normal stresses are arrived at. The SRC along the surface due to 
these normal stresses are then calculated on the same lines as in 
chapter 4. 

5.4.2 Analysis 

The elastic continuum approach is used to obtain the 
vertical displacements. The loaded area is discretized in the 
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same manner as in chapter 4. For tho 

convenience of applying 

boundary conditions for the beam (strip), the strip is discretized 
as shown in Fig. 5.5. Due to symmetry, only half the length of 
the strip is considered. The strip is divided in such a way that 
two nodes are available at the two ends as shown in the figure. 
The other nodes are placed such that all the nodes are equidistant 
from each other and each node forms the centre of an element. 
Hence, excepting the first and the last elements, all other 
elements arc of the same length, dl. The first and the last ones 
are half -elements. The number of ntxies, N = L /dl + 1. 

r 


Using Boussinesg^s equation, the vector of vertical 
displacements of all the N nodes due to the surface load as 
expressed in Eq. 4.5 of chapter 4, is 




(5.21) 


The vertical displacements of all the N nodes due to the 
normal stresses mobilised at the interface is cAtained using 
Mindlin''s solution, and the vector of the displacements as defined 
in Eq. 4.14 of chapter 4, is 




(5.22) 


where ^ ^ ^ the vector of normal stresses. The meanings of 
vectors matrix j remain unaltered. In the 
case of an inflexible strip the beam displacement vector was equal 
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to the rigid body displacement vector (Eg. 4.15) 

The total eoil dieplaceaent vector, | p' ^ ,he difference of 
Eqs. 5-21 and 5-22 expressed as 


or 







['■](') 


<5. 23) 


<5.24) 


For calculating the strip displacementp the beam equation is used. 
The loading on the beam is the normal stresses mcrtiilised at the 
soil -strip interface. The beam equation is 


I = M 

"d 


<5.25) 


where w is the beam deflection in the z direction, M is the moment 

and El is the flexural rigidity of the beam. The moment, M, is 

r r 

produced due to the stresses acting over the length, of the beam 
away from the point. The moment of the normal stresses about node 
i <Fig. 5. A) is 

dA dl 

M. = &. — — + 

2 4 

where x is the distance between node i and the centre of the area 

over which the normal stress, is acting and dA=2B^dl is the 

area of the element. 


H 


I 


& dA X 


<5.261 
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The ver'ical displacement of an, node i, from Eq. 5.25, expressed 
in thE finite difference form is 


, w. -2w. +w. . 

: I f— l::! h — ^ „ 


(5.27) 


Substituting Eq. 5.26 in 5.27 one obtains 


w. -2w. +w. 


E I 

r r 


^ w. .«-w . ^w. .. ax 

I J ' - - ] 

N-1 

J cxj^dl 2B^ j^(k-i)dlj + o-^ — 2B^ 


4 

dl dl 


k = i + i 


(5.28) 


Eq. 5.28 is written in matrix form as 


E I 

r r 


dl 


[■■]{■}■['.](•) 




(5.29) 


where I I, I is the finite difference coefficient matrix and 
is moment coefficient matrix. 




The boundary conditions for Eq. 5.25 are as follows. 

At X = O, slope, ^ = 0 
and 

at X = L , moment, ^ = 0 

dx 

The displacement of node 1 expressed in the finite difference form 
is 


(5.30) 


(5.31) 



IBl 


w -2w '♦•w 
2 1 


■ W •V 


(5.32) 


where w^ (Fig. 5.7) is the displacement of the node 2* at a 
distance dl to the left of node 1. From the boundary condition at 
node 1 (Eq. 5.30), one obtains 


w = w 
2 2 


(5.33) 


Substituting Eq. 5.33 in Eq. 5.32 the equation for displacement of 


node 1 is 


-2w +2w 


I r ^ = n 

rfl* i ‘ 


(5.34) 


The second boundary condition will be applied later in terms of 


the normal stresses mobilised. Hence the matrices 
are of size N-lxN while vectors ^ ^ 


[ ^ ] ““ [ k ] 

are of size N 


where 




-2 2 0 0 
1-210 
0 1-21 


0 1-21 


Non-dimensional ising the length parameters with the half- width of 
the loaded area, Eg. 5.29 is written as 




(5.35) 
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Fig. 5.5 Discretization of the Strip 
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Fig. 5.6 Normal Stresses 



Fig. 5.7 Node 2'' for Calculation of Moment at Node 1 



of the 


E I 

where K = — ^ ^ 

(g-> 

Eq. 5.35 gives the vertical displaceeent vector | w/B, 

points along the beam. Satisfying the coepatibility of 
displaceoents at the nodes by equating the net soil displaceoents 
and the strip displacenents the following equation is obtained 




C5.36) 


At this juncture, a dimensionless parameter, K = =r^- — is defined 

e 

as the flexibility ratio which takes into account the flexural 
rigidity of the strip. Eg. 5.36 is rewritten as 




(5.37) 


where 


0 '/^ ^ and 


is of size N-lxN and 


arc vectors of size N. 


The second boundary condition defined by Eg. 5.31 is now applied 
in terms of the stresses, calculating the moment at 

node N, stresses mobilised along the vrihole length of the strip, 
i.e., 2L are to be considered. For every element i, there is an 

t 

image element i"' over which the stresses are the same. The moment 



as 


at the 


node 


is written 


"h = ".T 28,[(‘N-l)dl-51 ] * [2L^-,(N-ndl-Sl,]] 

H- 4 

J o-.dl 2B^ J^(N-i)dl + [^2L_^-(N-i)dljj + 


^ ^ K-'4 ]] 


(5.38) 


Boundary condition 2 (Eq- 5.31) requires H to be zero at x=L 

: 

Hence equating Eq. 5.36 to zero one obtains 


H~ 1 


2 


O'. 

V 



(5.39) 


Eq. 5.39 is expressed in dimensionless form as 



i *a 


<7 

0.5 ( — - 
q 



(5.40) 


Eq. 5.40 implies that the force equilibrium for the strip is also 
satisfied. Substituting Eq. 5.40 in Eq. 5.37, N equations for the 
nornalised normal stresses, o'/q J are obtained. Once the 

stresses are obtained, the displacements are calculated using 
either Eq. 5.22 or Eq. 5.35. The settlement reduction coefficients 
(SRC) of the points along the surface are calculated in the same 
manner as explained for rigid inflexible strips in chapter 4. 

5.5 Coupled Effect (Horizontal and Vertical Rigidity) 

Rigid Single Strip 

Thus far, the rigidity of the strip in the horizontal and 



directions were considered separately. The horizontal 
and vertical displaceiwnls of the nodes alons the strip, 
independently due to the shear and normal stresses apart fro. 
those due to the surface load were considered. In this section 
the influence of the shear stresses in causlns vertical 
displacements and normal stresses in causina horizontal 
displacements are also coupled. The reinforcina strip now acts as 
an inextcnsible and inflexible one. 

5.5.1 Prc*lciB Definition 

The definition of the problem is the same as that for a 
single strip, i-e., a strip of size, 2L x2B is placed centrally 

r T 

below a rectangular area transmitting an uniform load of intensity 
q, at a depth U^. The strip being rigid does not experience any 
elongation while it undergoes a uniform vertical translation due 
to the soil displacements. 


5.5.2 Analysis 

The discretizations of the loaded area and the reinforcing 
strip are the same as explained in chapter 3 or 4. The vectors of 
the horizontal and vertical displacements of the N nodes along the 
half-length of the strip due to the surface load, computed from 
Boussincsq'^s equations are 


and 






(5.41) 


(5.42) 


the mobilised 


The horizontal displacements of the N nodes due to 
shear stresses | j. and the normal stresses | a J- 

Mindlin's equations are 


computed 


from 


and 


{ ^ [ 'Ir ] { - } 

s: 


(5.43) 


(5.44) 


Similarly, the vertical displacements of the N nodes due to the 
mobilised shear stresses and normal stresses obtained from 
Hindi in^s equations are 

s 

and 

{ C } ‘ -F ['-]{" } 

It is to be remembered that, as explained in section 4.2.1 
the nor ami stresses considered are the difference of the stresses 
acting on the top and bottom faces of the strip. 


Combining Egs. 5.43 and 5.44, the vector of horizontal 
displacements of the N nodes due to the shear and normal stresses 
are 





(5.47) 



where | | is a vector o{> size N, | p J is a stress vector of 

size 2N and j is a matrix of size N x 2N. Further 

{p} = Ll^j 

f 1 = f i’' *> I*" 1 

j_ X J L X'T x& j 

Combining Eqs. 5.45 and 5.46, the vector of vertical displacements 
of the N nodes due to the stresses is 

B 

where { < } * { } ^ { } 

and 

Satisfying the compatibility of displacements for the nodes as 
Ci) the net horizontal displacements equal zero, and (ii) the net 
vertical displacements are equal to the rigid body displacement, 

6 , one obtains 
o 


(5.52) 


(5.53) 


(5.48) 


(5.49) 


(5.50) 


(5.54) 


and 


{ ■ "o{ ‘ ■ 


(5.55) 



Hence 


and 




(5.56) 


(5.57) 


Eq. 5.56 gives N equations for 2N stresses (N values of r/q j- 

and a/^ ] each) and Eq. 5.57 gives N equations in 2N stresses 

and a rigid body displacement. Combining Eq. 5.56 and 5.57 and 
satisfying the force equilibrium relation as in section 4.2.2, one 
obtains 2N+1 equations in N values of shear stresses, r/q j-, and 

normal stresses, o-Zq each, and a rigid body displacement, 6^, 


as 




(5.58) 


where 




'"I 

I 


r t"" • T*' 1 

XT ■ x<y 



X 


[‘ J 


!*• 


r I 



L Z 


1- *T xty J 


(5.59) 


and 


{ I' f = [ulii J 


(5.60) 


These equations are solved for the stresses on and the uniform 
displacement of the strip. Once the stresses arc obtained the 
settlement reduction coefficients of points along the surface are 
calculated as explained in sections 3.2.2 for shear stresses and 


4.2.2 for normal stresses. 



5*5-3 Slip at the Ifitcrfacc 


It IS essential that the shear stresses moblised at the 
strip -soil interface do not exceed the frictional resistance 
available at the interface due to the normal stresses and 
overburden stresses at that depth so that slip is prevented. This 
condition may be expressed as 

f O' ^ T 

^ ^ 5 . 61 ) 


where fj is the coefficient of friction at the strip-soil 
interface, — ^ — "TT” normalised overburden stress, and o/q 

and T/q are the normalised normal and shear stresses mobilised at 
the interface. 


If this condition is not satisfied, the shear stresses are 
limited to the frictional resistance as 


T 



C5.62) 


Eq. 5.61 is checked for stresses at all the elements. 
Wierever the condition is violated, Eq. 5.62 is applied. The 
shear stresses on elements where slip occurs are now known and 
these values are substituted in Eq. 5.58 which is solved for the 
stresses on the remaining elements. 



5-6 Result* and Discussion 


The effects of the extensibility and flexural rigidity of 
the reinforcing strip on the stresses mobilised and settlement 
reduction of the points on the surface, are studied herein. As 
discussed in the section on analysis, the net displacement of the 
soil at the soil -strip interface is equated to the elongation or 
the flexural deformation of the strip, as the case may be. The 
discretizations of the loaded area, and of the reinforcing strip 
for shear interaction arc the same as in section 3.6. The strip 
is discretized for flexure differently as discussed in section 
5.4.2. 


5.6.1 Extensible Strip 

As in the case of a rigid strip (Chapter 3), the stresses 
mobilised at the interface are evaluated at the centre of each 
element. Apart from the parameters studied in section 3.6, the 
elongation ratio, K , defined as K =E t /E B, is also considered, 
where E and t arc the elasticity modulus and the thickness of 

r r 

the strip respectively- K =0 indicates absence of reinforcement 

T 

(unreinforced soil) implying zero stress mobilisationr while K^=oo 
implies presence of a rigid strip. The range of the ratiOr K^r in 
the parametric study is 5 to 5000. The upper limit of is 

restricted to 5000 because the results obtained for this value 
agree well with those for a rigid strip considered in chapter 3. 
Further, for this value of K^, the elongations are negligibly 
small. The Poisson's ratio, is maintained at 0.3. 

Fig- 5.8 presents the variation of the normalised shear 
stress, T/q, with the distance, x/B^, along the half-length of the 



strip for different values and tor a strip of length, L /B =2, 
placed below a square loaded area with an aspect ratio, L/b'= 1 at 
a depth, For all values of the elongation ratio, K^, the 
shear stresses are a eaxinum at a distance, x/B=l, i.e. below the 
edge of the loaded area. For a highly extensible strip with K =5 

T ' 

the stresses mobilised are the least, with the edge of the strip 
showing negative stresses. With an increase in the K value the 
stresses increase and are positive throughout the length. For 
strips with K^=5000, the shear stresses are the same as those for 
a rigid strip considered in chapter 3. The maximum shear stress, 
x/g, is 0-B in this case, and is observed at a distance, x/B =1. 

t 

For strips with K^=5, a small negative stress is observed near the 
centre. This negative stress is possibly mobilised due to two 
conflicting criteria. Due to symmetry of the problem, the 
displacements at and near the centre have to be zero or extremely 
small. However, since the tensile force at the centre is maximum 
and the strip is highly extensible, the displacement of node 1 
could be significant. Thus the negative stress is a consequence of 
the need to satisfy these two mechanisms. The net lateral 
displacestents for an extensible strip are high owing to the high 
elongations experienced by the strip. Hence, the displacements 
which the shear stresses have to counteract are less resulting in 
low shear stress values. As increases, the elongations in the 
strip decrease, consequently the displacements to be resisted by 
the interacting or mobilised shear stresses increase 

significantly. 


The effect of the depth of placement, of a strip with 

K^=50, and length, L^/Bj=2, placed below a square area with. 



L,/B,= 1 . on the normal ined sheer Mresnes, r/,, ie depicted in 

Fig. 5.9. This figure shows a trend similar to the one seen in 

Fig. 3.M for a rigid strip. A maeminue, stress of T=0.672g is 

observed for strips placed at depth U /n -i x 

neptn, u /B -1 at a distance, x/B =1 

‘ f 

along the length of the strip. In <-nn*.ra*r+ ^ - -j x - 

=*'-rip. in contrast, rigid strips show a 

maximum stress of 0.8q at the same depth (Fig. 3.14). Because of 
the eetensibilitr of the strip, the net displacements are 
relatively large resulting in smaller shear stress mobilisation at 

the soil -strip interface. 


The variation of the tensile forces, T/gB^, induced in the 
strip of length, L /B ==2, placed under a square area with L/B=l, 

r I f f 

at a depth, U^/B^=l, for various values is presented in Fig. 
5.10. The tension for the strip with, K^=5 is the least owing to 
low shear stress mcAiilisation. Strips with K =5000 exhibit a 

T 

tension of O. lOSqB* at the centre which compares very well with 
that for a rigid strip (Fig. 3. IB). For highly extensible strips 
with K,^=5, the tension at the centre is 0.038qB^. Near the edge, 
these strips experience a small compression due to the 
mobilisation of negative stresses there at (Fig. 5.8) 

Fig. 5.11 depicts the effect of the elongation ratio, on 
the elongation, of the strips of length, L^/B^=2, placed 
below square area (L^/Bj=l) at a depth, As expected, strips 
with an elongation ratio, show maximum extension. The total 
elongation, <5E^/Bj.q, of these strips is of the order of 0.067. 
The elongation is maximum for these strips at a distance, x/B^=1.5 
along the length of the strip. Beyond this distance, it decreases 
marginally as the strip is subjected to a small amount of 
compression in this zone. Elongations of the strips decrease with 
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Fig. 5,8 Variation of Shear Stress with Distance 
for Different Elongation Ratios 
(U/Bf=2. Uo/Br=1, »^s=0.3) 


Fig. 5.9 Va notion of Shear Stress with Distonpe for 
Different Depths of Plocennent of Strip 
(U/B,=2, U/B,= 1, Kt=50, ^s=0.3) 
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FiQ. 5.10 Variation of Tension with Distance 
^ for Different Elongation Ratios 

(U/Bf=2. Lr/B,= 1. Uo/Bf=1. J^s=0.3) 
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Fig. 5.11 


Variation of Elongation 


as can be anticipated. 


increasing values of K 

The variation of SRC with distance, along the surface 

for various values for strips of length, L^/B^=2, placed below 

a square area at a depth, is presented in Fig. 5.12. The 

trend for all the curves in this figure arc similar to that for 

rigid strips placed at the same depth (Fig. 3.203. The SRC at the 

centre (x /B =03 and at the edge of the loaded area (x/B=13 for 
11 f f 

a highly extensible strip with K =5 are 0.0014 and 0.0001 

T 

respectively while those for a strip with K =50 are 0.00285 and 

T 

0.0007 respectively. As K increases, the SRC values tend to 

T 

those (0.00335 and 0.00093 of a rigid incxtensible strip. The 
increase in SRC values with is a consequence of the increase in 
the shear stresses with K . 

T 

The effect of the elongation ratio, K , on the settlement 

T 

reduction coefficient at the centre, I , for various depths of 

SC 

placement, U^/B^, and for strips with length, L^/B^=2, placed 
below a square area is presented in Fig. 5.13. I is a maximum 

SC 

in the depth range, U /B. of 0.75 to 1 depending on the values of 

o f 

K . For low values of K (K =53 a maximum I value of 0.0015 is 

T T T SC 

observed at a depth U /B=0.75. At very shallow depths, 

o f 

U /B=0.25, the I values are practically zero. As K increases, 

o f sc ^ 

the I values at all depths tend to those of a rigid strip- The 
sc 

maximum I value for strips with K =5 are around 0.45 times that 
sc T 

for a rigid strip. 


Fig. 5.14 shows the plot of the settlement reduction 
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gc' centre of the loaded area versus the 

length, L /B.r of the strip for various aspect ratios, L /B of 

r 1 f f 

the loaded area and for a depth ratio, U /B =1 and elongation 

o f 

ratio, K^=50. It is evident from the figure that the length of 

the strip, L^/B^, and aspect ratio, of the loaded area 

affect the I values considerably. The trends in this plot arc 
»c 

similar to the ones observed for rigid strips. The I values, in 

8C 

this case have a smaller magnitude as compared to rigid strips. 
This is because, the stresses mobilised at the interface are 
smaller. The explanations for the trends observed, given for a 
rigid strip hold good in this case also. 


The effect of the elongation ratio, K^, on the for 

strips located at various distances from x-axis, 

strips of length, L^/B=2r placed below a square area at a depth 

U /B=l is studied in Fig. 5.15. Here it is to be remembered that 

o f 

since two strips are placed off-centre, the half-width of the 
strips considered are 0.0258^. and the values presented are for 

one strip only. The trends in the curves are very similar to 
those observed for rigid strips tFig. 3.34). The distance, S^/B^, 
where the maximum occurs is a function of K^. For strips with 

small values, this distance is 0.15 whereas for rigid strips it 
is 0.25. At large values (S^/Bj.>1.5) the values 

negligibly small. The maximum value of for strips with K^-5 

is 0.00075 which is around 0.33 times that for a rigid strip. 

The effect of the length of the strip, 

for an extensible strip to that of a rigid one, 
sc 


of the I 





^«c<*/^*c<r>' various values and for strips placed below 
square areas at a depth, U^/B^=l, is depicted in Fig. 5.16. The 
ratio, increases with and tends to a value of 1 
indicating that the strips behave as rigid ones. At low K values 

TT 

(K^=5), the values of the ratio are 0.56 and 0.4 for strips of 
length, L^/B^ equal to 1 and 2 respectively. Based on this plot, 
^sc for extensible strips can be calculated from those for 
rigid strips by choosing the appropriate correction factor. 


Based on the plots discussed so far it is clear that the 
variation of tension in the strip or the settlement reduction 
coefficient for an extensible strip follow the same trends as 
those for a rigid inextensible one. It is only the magnitude of 
these values that differ depending on the values of the elongation 
ratio, K , and the length, L . From Figs. 5.8 and 5.9 it is 

* r 

evident that the stresses mobilised are small for extensible 
strips resulting in low values of SRC and tensions in the strip. 
For design purposes, the SRC values for extensible strips may be 
obtained by multiplying the corresponding values for a rigid 
inextensible strip by an appropriate correction factor based on 
the value of the elongation ratio, K , and the length of the 

T 

reinforcement. This factor is low for very low values of K . For 

T 

strips with the ratio, 500>K^>100, a correction factor of 0.9 may 
be used for any length of the strip. For strips with K ^500 the 

T 

correction factor can be taken as 1.0. 


5.6.2 Flexible Strip 

The strip, as mentioned in section 5.4 is modelled to behave 
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as a beam. The second order differential equation given in Eq. 
5.25 is solved using the finite -difference technique. The stresses 
and SRC obtained with dl=0.1B^ were compared with those obtained 
with a finer discretization of 0.056^.. There wasn't any 
significant improvement in the values as a consequence of this 
finer element size. The nodes were selected in such a manner that 
the boundary conditions for the beam may be satisfied without any 
difficulty. It is to be noted at this Juncture that due to edge 
effect, the stresses of a point at the edge of the strip tend to 
infinity. Hence, if the node is considered to be exactly at the 
edge of the strip, i.e, at x=L , then the stresses resulting from 

r 

the solution of the simultaneous equations given by Eq. 5.37 are 
highly sensitive to the value of the stress at the edge of the 
strip which is very high. As a result, the variation of stresses 
along the strip is erratic. In order to eliminate or minimise 
this inconsistency, the last node i.e, the node is considered 
to be not at the edge but at a distance of dl/4 from the edge 
(Fig. 5.17). The stress at this point is definitely a finite 
quantity and hence the solution obtained for normal stresses is 
stable. For the beam equations the boundary conditions for x=L 

T 

are applied at node N which is just at a distance of dl/4, i.e, 
0.025B^ (O. lBj,/4) from the edge. 
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Fig. 5.17 Discretization for Beam 
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The effects of the various parameters mentioned in chapter 

4, are studied on the normal stresses mobilised and the SRC at the 

surface. The main parameter considered in addition to these, is 

the flexibility ratio of the strip, K , defined as E I /£ 

* r r s f 

indicates no (or infinitely flexible! reinforcement while 
*'ep>'esents a rigid strip. The range of K considered in this 
analysis is 10 * to 10*. The upper limit of K is limited to 10* 

B 

since for this value the strip behaves as a rigid one. All 

parameters are studied for a Poisson '■s ratio, v =0.3 

s 

Fig. 5.18 presents the variation of the normalised normal 
stress, cx/q, with the distance, x/B^., along the strip for various 
values of the flexibility ratio, K , and for strips of length, 

B 

L /B =2 placed below a square area at a depth U /B =1. The 

r f Of 

stresses for strips with low K values (K =10~*> are practically 

B B 

zero. As K increases, the stresses increase and tend close to 

B 

those observed for rigid strips. For intermediate values (K =1), 

B 

the stresses, as expected are maximum at the centre (x/B^=0) and 
decrease with distance, x/B^. The stresses change sign close to a 
point at a distance, x/B^=l, i.e near the edge of the loaded area 
and tend to very large negative values at the edge of the strip. 
The peaking of stress is due to the edge effect. For strips with 
K <10*, the stresses are less compared to those for a rigid strip 

B 

and consequently the edge stresses are also smaller. For strips 
with low K values the flexural deflections are high. These 

B 

deflections tend to the vertical displacements of the soil due to 
the surface loading alone (as K decreases). The difference 

B 

between the soil displacements and the strip deflections being 
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small, the normal stresses mobilised in order to counteract this 
^ in displacements is small. As increases, the 

flexural deflections decrease and for K^>10*, the deflections are 
uniform as for a rigid strip. The difference between the soil 
displacements due to applied surface loading alone and the actual 
strip deflections are high, and hence, the stresses mobilised also 
are high. 

The effect of the depth of strip, U /B , on the normalised 

o f 

normal stresses, o/q, for strips with K =1 and length L /B =2, 

B r f 

placed below a square area CL^/B^=11, is depicted in Fig. 5.19. 
The variation of the stresses with distance, x/B^, along the strip 
is similar to those for a rigid strip (Fig. 4.11). Maximum 
stresses are observed for strips at shallow depths. At a depth 
U^/Bj,=0.25 the stress, a/q, at n/B^—0 is 3.1 while for a rigid 
strip it is of the order of 4.25. With increasing depth of the 
strip the stresses decrease and consequently the peak negative 
stresses at the strip edge also decrease. 

Fig. 5.20 presents the variation of the strip deflections, 
wE /B q, with the distance, x/B , along the strip for various K 

S f f B 

values, and for strips of length, L /B =2, placed at a depth, 

r f 

U /B=l, below a square area. It is seen that for strips with low 
o f 

K values (K =10~^) the deflections tend to the displacements of 

B B 

the soil due to surface loading. At the centre, the deflection. 


wE /B q, is 1.34 while at the edge it is 

u f 

0.6. 

For 

strips 

with 

K =1, the corresponding displacements 

B 

are 

1.12 

and 

0.82 

respectively. With increasing to 10*, 

the 

deflections 

equal 
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Fig. 5.18 Variation of Normal Stress^ with 

Distance for Different Flexibility Rotios 
(U/Bf=2, U/Bf=1, Uo/Bf=1. i/s=0.3) 
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those for a rigid strip and are equal to the rigid body 
displacement, <5 E /Ba, of 0.99. 

OS I 

The effect of flexibility ratio, K on the settlement 

B 

reduction coefficient, SRC, on the surface for strips of length, 
L placed below a square area at a depth, U /B=l, is shown 

in Fig. 5.21. The SRC is fnaximum at the centre and decreases with 
the distance, x /B , along the surface for all K values. For 

if B 

strips with very low values of the flexibility ratio (K =10“*), 

B 

the SRC is negligibly small. This is obvious because the stresses 
are very small (Fig. 5.18) for the highly flexible strips. As K 

B 

increases, the flexural rigidity of the strips become prominent, 
the stresses mobilised increase and consequently the SRC values 
increase. For strips with K =10*, the SRC values coincide with 

B 

those for a rigid strip (Fig. 4.16) with the centre and edge 
exhibiting SRC values of 0.048 and 0.013 respectively. The 
corresponding values of SRC for strips with an intermediate value 
of the flexibility ratio (K =1) are 0.29 and 0.008. The SRC at the 

B 

centre for strips with Kg=l is about 0.6 times that of a rigid 
inflexible one. 

Fig. 5.22 presents the effect of on the SRC at the centre 

of the loaded area, I , for different depths, 1) /B,, of the strip 

sc o f 

of length, L^/B^=2 placed below a square area. The trends of the 

curves are again similar to those of a rigid strip (Fig. 4.20) 

with maximum I being observed at shallow depths. The value 

decreases rapidly with depth of placement of strip, since the 
stresses mobilised decrease and the influence of the strips in 


Fig. 5.20 Variation of Strip Deflection with Distance 
for Different Flexibility Ratios 
(U/Bf=2, U/B,= 1. Uo/Br=1, i's=0.3) 


Fiq. 5.21 Variation of SRC with Distance along 
Surface for Different Flexibility Ratios 
(U/Bf=2. U/Bf=1, Uo/Bf=1, i/s=0.3) 
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causing settlement reduction reduces. Highly flexible strips show 
negligibly small values. Strips with K =10* show an I value 

* SC 

of 0.19 at a depth, U^/B^=0.25, which is equal to the value for a 

rigid strip. As the depth increases, the effect of K on I 

B sc 

decreases. At great depths, (U /B.il.5), 1 values are very 

small. The value for strips with K^=l is 0.135 at a depth, 

25, which is about 0.7 times the value for a rigid strip. 

The effect of the aspect ratio, of the loaded area on 

values, for different lengths of the strip, L /B , and for 
strips with K =1, and placed at a depths U /B=l is depicted in 

B Of 

Fig. 5.23. It is clear from the figure that I values increase 

sc 

with the length of the strip, L /B., and the aspect ratio, 

r 1 f f 

of the loaded area. As observed in Fig. 4.21 longer rectangles 
show higher I values for all lengths of the strip. The increase 

SO 

in the values from L^/Bj.= l to 2 is significant while the 

increase is less for further increase in the aspect ratio from 2 

to 10. For rigid strips (Fig. 4.21) it is seen that I values 

sc 

increase monotonically with the length of the strips, L /B , while 

r f 

for strips with K =1 the increase in I is gradual beyond a 

B SC 

length, L /B, of 2. The reason for this may be attributed to the 

r f 

low normal stress mobilisation for flexible strips. The I values 

SC 

are less compared to those for rigid strips. For the length of 
strips, L /B =5, the I values for aspect ratio, L/B=l and 10 

r f SC 1 I 

are 0.035 and 0.053 respectively, while the corresponding values 
for a rigid strip are 0.125 and 0.245 respectively. Thus for 
strips with K =1 the I values for L /B =1 and 5 arc about 0.28 

B fiC r f 

times and 0.22 times the values for a rigid strip. 
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The variation of with the distance of the strip, 
from the x-axis for various values and for strips with length, 
L /B =2 placed below a square area at a depth, U /B =1, is 

r f Of 

presented in Fig. 5.24. It is to be remembered here that, for 
strips placed of f -centrally the half-width of the strips, B^ is 
0.025B^. Further, the settlement reduction is calculated for one 
strip only. For the estimation of the effect of both the strips 
on the settlement reduction, these values are to be doubled. In 
the figure, the trends are similar to those observed for rigid 
strips CFig. 4.27). Maximum I values occur for strips at a 

SC 

distance, S /B=0.25, from the x-axis. As S /B increases beyond 
y f y * 

0.25, I drops rapidly. The reason for this trend is as discussed 

SC 

for rigid strips in Fig. 4.27. Strips with very small values 

show negligible I^^ values. I^^ for strips with intermediate 
values IK^=l>r placed at a distance Sy/Bj=0.25 is 0.021 while for 
strips with K =10* it is 0.0295 which coincides with that of a 
rigid strip. 


The effect of the length, L /B,, of the strips on the ratio 

r f 

of I for a flexible strip to that of a rigid one, ^sca/^sc(r> 

SC 

for various flexibility ratios, is presented in Fig. 5.25. It is 
evident that the length of the strip affects this ratio 

considerably. For shorter strips with L^/Bj.=l, the ratio, 

I /I >0.9 for K >0.5 while for strips with L /B =5 it is 

SC(f> SC<r> B 

equal to 0.9 for K =100. For K =1, strips of length, 
rigid while the ratio is around 0.2 for strips of length L^/Bj=5. 
The reason for this is that for shorter strips the difference in 
vertical displacement at the centre and edge of the strip due to 
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Rg. 5.24 Effect of Distance of Strip from x-axis 
po Isc for Different Flexibility Ratios 
(Lr/Bf=2, Uo//Bf=1, Lf/B/=1, i's=0.3) 


Ra. 5.25 Effcet of Flexibility Ratio on Ratio of 'sc fo'' 

^ Flexible and Rigid Strips, for Different Strip Lengths 

Uo/Bt=1. Vs=0.3) 
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surface loadins is small and hence the flexural rigidity required 
to deform the strip uniformly is small. On the other hand, this 
difference in displacements for longer strips is high and hence a 
higher flexural rigidity is required to cause the strips to 
displace uniformly. This figure can be used as a design chart. 
Knowing the flexibility ratio of a strip, the settlement reduction 
at the centre can be evaluated by multiplying those values for a 
rigid strip by a suitable correction factor obtained from this 
figure. 

From the illustrations seen thus far, it is clear that for 
very flexible strips the stresses mobilised are very small. The 
normal stresses tend to those for a rigid strip with increase in 
the flexibility ratio. Consequently, as K increases the 

displacements at the soil-strip interface also change from those 
of soil displacements due to the applied loading alone to a 
uniform rigid body translation. The variations in the calculated 
quantities, viz., stresses and settlement reduction coefficients 
for the flexibility ratio ranging from 10 * Cfor a highly flexible 
case) to 10* ta rigid strip) are similar to those considered in 
chapter 4. It is also perceived that the relative rigidity of the 
strip does not only depend on the flexibility (E^I^) but also on 
the length of the strip, L . K =1 is sufficient to make a strip 
of length, L /B=i behave as a rigid one while a value of K^=500 
is necessary to make a strip of length, L^/B^=5 behave as a 
sufficiently rigid one. Fig. 5.25 can be conveniently used for 
strips below square areas at a depth, U^/B^=l, to predict the 
values for flexible strips from those of a rigid one. 



5,6t.3 Coupled Shear and Normal Stress Interactions 

The coupled effect of the shear and normal stresses for a 
rigid strip arc considered. This implies that the effect of the 
shear stresses on vertical displacements and vice-versa are also 
accounted for. The shear and normal stresses mobilised at the 
interface arc evaluated. 

The variation of the normalised shear stress, T/q, with 
distance, x/B^ along the length of the strip of length, L^/B^=2, 
placed at various depths, U^/B^, below a square area tLj./B^=l) is 
presented in Fig. 5.2it. It is evident that the shear stresses 
follow exactly the same trend as for the case where the shear 
interactions alone were considered (Fig. 3.14). The magnitude of 
the stresses are marginally smaller as compared to Fig. 3.14 but 
for all practical purposes it may be stated that there is no 
difference in the magnitudes of the shear stresses. 

Fig. 5.27 depicts the variation of the normalised normal 
stresses, o/q, along the length of the strip described in the 
previous paragraph. Here again, it can be seen that the stresses 
follow the same trend as for the case where the normal stress 
interactions alone were considered (Fig. 4.11). The difference in 
the magnitudes of the stresses, too is marginal. From Figs. 5.26 
and 5.27 it is clear that the effect of shear stresses in causing 
vertical displacements and normal stresses in producing horizontal 
displacements is negligible. It, hence follows that mechanisms 
considered individually result in accurate values of shear and 


normal stresses. 


Normal btress a/q 


Fig. 5.26 Variation of Sheer Stress with Distonce for 
Different Depths of Strip (Coupled Effect) 
(U/Br=2. U/B,= 1. ^5=6.3) 



Fiq 5 27- Variation of Normal Stress with Distance for 
Different Depths of Strip (Coupled Effect) 
(U/B,=2. U/B,= 1. ^5=0.3) 
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The effect of both the stresses in causing reduction of 
settlements of points along the surface is depicted in Fig. 5.28. 
This figure presents the variation of SRC with distance, 
foT strips of various lengths, L /B , placed below a square area 

r I 

at a depth, U^/B^=l. The curves follow the same trend as those 
for a rigid strip undergoing uniform vertical displacement 
(Fig. 4. 16>. The magnitudes of the SRC values are also the same 
as those observed in Fig. 4.16. In fact, the magnitudes are equal 
to the summation of the SRC values for a rigid strip where shear 
(Fig. 3.20) and normal (Fig. 4.16) stresses are taken into account 
separately. It so turns out that the SRC values due to shear 
stresses are much smaller than those due to normal stresses and 
hence the SRC values for the coupled case are very close to those 
due to normal stresses alone. 

It can hence be concluded that settlement reduction of 
points on the surface for a rigid strip result predominantly due 
to the normal stress mobilisation. For design purposes it is 
sufficient if SRC values due to normal stresses alone arc con- 
sidered. Further, these values will be on the conservative side. 

5.6.4 Slip at the interface 

Fig- 5.29 shows the variation of the normalised shear 
stress, T/q, with the distance, x/B^, for different coefficient 
of friction values, (Jr for an overburden stress ratio, j'Bj/q=0.5. 
The results arc for strips of length, L^/Bj=2, placed at a depth, 

U /B =1, below a square loaded area. The dotted line depicts the 

o f 

shear stress variation for the coupled case without slip. It is 
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Fig. 5.28 Variation of SRC with Distance along Surface 
■for Different Strip Lengths (Coupled Effect) 
(Uo/Bf=1, Lf/Bf=1, ^5=0.3) 
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Rg. 5.29 Variation of Shear Stress with Distance 

for Different Coefficient of Friction Values 
(U/B,=2. U/Bf=1. Uo/B,=:1. yB,/q=0.5. Vs=0.3) 




seen that a portion of the strip beneath the edge of the loaded 
area , i.c., at x/B^= 1, predominantly experiences slip. Shear 
these points reduce to low values because the normal 
stresses in this region are very small and the condition that 
shear stresses cannot exceed n times normal stresses (given by Eg. 
5.61) is to be satisfied. For a low coefficient of friction value 
(A^0*25), the strength available is small and slip occurs over a 
large distance (0. 4<x/B^^l .5) . Further, the stresses on elements 
which do not fail are higher than those observed for the no slip 
case. The stresses on the elements along which slip occurs are 
forced to the limiting values resulting in a redistribution of 
stresses on the adjacent elements. As a result the stresses on 
the adjacent elements increase. As expected, the shear stresses 
are higher for higher coefficient of friction values. 

Fig. 5.30 shows the normalised shear stress variation for 
different overburden ratios, j^B./q, for strips of length, L /B,=2, 

f r I 

placed at a depth, U /B,=l, below a square surface loaded area. 

o f 

The results are for a coefficient of friction, ju, of 0.5. The 
elements within the region, 0.65x/B^51.5, depending on the 
overburden stress ratio, experience slip. The region shifts to 
the left with an increase in the value of j^B^/q. The stresses 
along elements which do not undergo slip are once again higher 
than those for the no slip case. 

The variation of the normalised shear stress, T/q, with 
distance for various depths of placement of the strip is presented 
in Fig. 5.31. The results presented are for a coefficient of 
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Fig. 5.30 


Variation of Shear Stress with Distance 

for Different Overburden Ratios 

(U/Bf=2. Uo/Br=1. ^=0.5, Vs=0.3) 


Fig. 5.31 


Variation of Shear Stress with Distance 
for Different Depths of Strip Placement 
(U/Bf=2, U/Bf=1, /^=0.5, yB(/q=03, i^s=0.3) 
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friction, m= 0.5, and overburden ratio of 0.5. At a depth, 
0^/B^=0.25, no element along the length of the strip experiences 
slip. For a strip at depth, U /B=0.25, the value of shear stress 

O f 

beneath the edge of the loaded area is very small unlike the 

values of shear stresses for strips at depths greater than 0.25. 

Also, the normal stresses for strips closer to the surface are 

higher. Hence, none of the elements of the strip at U /B=0.25 

of 

experience slip. With increasing depths of the strips the normal 
stresses tend to become somewhat uniform and the points where the 
shear resistance is least move closer to the centre of the strip. 
Therefore, slip occurs along more number of elements and on 
elements relatively closer to the centre. 

Fig. 5.31 shows the variation of SRC with distance, 
along the surface for various values of the coefficient of 
friction, /j, for an overburden stress ratio, 2'Bj,/q=0.5. The 
results are for strips of length, L /B =2, placed below a square 

r f 

area at a depth, U =B . It is seen that there is a marginal 

o f 

decrease in SRC values with decreasing values of the friction 
coefficient, ^i. SRC at the centre of the loaded area (x^/Bj,=0) 
for a^= 0.25 is 0.046 which is 0.97 times that for the no slip case 
to. 0473). This shows that the slip at the interface does not 
affect the settlement reduction of points along the surface. It 
seen earlier that the settlement reduction due to shear 
stresses is much less than that due to normal stresses. Hence, 
slip failure along few elements below the edge of the loaded area 
and the consequent redistribution in stresses does not affect the 

It is thus evident that the failure of the 


settlement reduction. 



reinforcement in not mobilising the elastic shear stresses, does 
not lead to any significant decrease in the reduction of surface 
settlements. 

The shape of the r/q versus distance curves reflect the 
influence of slip along the reinforcement. Uhere slip occurs, the 
shear stress is limited to {Jt times the total normal stress. The 
difference between the elastic shear stress and the limiting value 
is redistributed to the adjacent elements. The redistribution 
effect does not appear in the shear stresses on elements farther 
away from the elements along which slip has occured. From these 
results it is clear that the elements on which the normal stresses 
are small, fail in shear. The reduction in shear stresses along 
these elements results in an increase in the shear stresses along 
the adjacent ones. The number and the distance from the centre, 
of the elements which fail by slip depend on the coefficient of 
friction, /n, the overburden stress ratio, and the depth of 

placement, of the strip. It is also noted that at greater 

depths, relatively a larger number of elements fail by slip, 
because the normal stresses are uniformly distributed and thus 
smaller in magnitude. 

5.7 Comparison with available literature 

Brown and Poulos (1981) have presented a finite element 
analysis of foundations on reinforced soil . The problem consists 
of a strip footing on the surface of a soil layer reinforced with 
strip reinforcements. Four layers of reinforcements have been 

flexural stiffness and transmitting axial 


considered with zero 



forces only. The parameters of their analysis expressed in terms 
of the parameters selected in this present analysis are as 
follows: aspect ratio of the loaded area, L^/Bj=10 (strip 

footing), length of the reinforcement, L /B =5, depth of the first 
layer, U /B=0.33, For the soil, E'/c''=1000 and v =0.3. The 

VI s 

authors have conducted the analysis for a linear density ratio 
(LDR) of 0.5. This corresponds to strips of width ratio, B /B = 

r f 

0,025 at a spacing of 0. 105B^ throughout the length of footing. 


For the parameters mentioned above, the settlement reduction 
coefficients at the centre of the loaded area (I ) are calculated 

8C 

for each strip considered independently. For the configuration, 
I , computed as the sum of the values obtained individually is 

SC 

0.4272. Choosing a correction factor for interference effect 
(Fig. 3.23) as 0.5, the net I =0.2136. The reduction in settle- 

SC 

ment at the centre of the footing is calculated from Eg. 3.35 as 




q 


This is valid for loading levels within the elastic range. At a 
loading level of q=16c'' the reduction is 


_ 0.2136xl6c'x2(l+0.3) 
_g_ roocE’ 


0.0089 


From Fig. 5.30 the corresponding settlement for the unreinforced 
case is about 0.0756^.. This implies a reduction in settlement of 
11.951. At a loading level of q=20c'', the reduction is 1351. These 
reductions are for a single layer of reinforcement. 



a depth 


For the second layer of reinforcements at 

^sc *=o"*P‘^ted as the sum of the values due to 
single strips is 0.166. Considering both layers and choosing a 
'factor of 0.4- for horizontal and vertical interference 
the net is 0.237. The settlement reduction for the stress 

levels of 16c' and 20c' are 137. and 14.57. respectively. The 
reductions obtained are for rigid strips as the elongation 
stiffnesses arc not known. The actual settlement reductions for 

elastic strips will be smaller than those obtained. The results 

» 

compare reasonably well with those reported (Table 5.1a). 

Shankar iah and Narahari (1988) conducted model tests on a 
sand bed reinforced with G. I. strips. The test details are shown 
in Table 5.1b. Reinforcements were placed in both directions, 
occupying an area of 43 cm. square. The parameters of the tests 
in non-dimensional form are: aspect ratio of the footing, Lj/Bj,= l, 
length of the strips, L^/B^=5, depth, U^/B^=0.47, width of each 
strip B /B =0.058. Spacing between each strip = 0.82B,. 

r f f 

The settlement reduction coefficient at the centre of the 
footing, I , due to shear stresses and normal stresses mobilised 

SC 

along the strips are calculated for single strips independently. 

I due to shear stresses for strips in a single direction is 
sc 

calculated as 0.0069. Considering both directions, the I is 

SC 

0.0138. Since the spacing between the strips is of the order of 
0.82B^, a correction factor due to interference of 0.9 is taken. 
Hence, the net I due to shear stresses is 0.0124. I due to 
normal stresses for strips in both directions is calculated as 



0.905. Choosing an interference correction of 0.8, the I is 

sc 

0.724. Based on the flexural ratio of the strips, a correction 

factor of 0.3 is taken which gives a net I value of 0.217 due to 

sc 

normal stresses. The total I due to horizontal and vertical 

SC 

rigidity is 0.2294. The reduction in settlements due to 

reinforcements is calculated from Eq. 3.35. The elasticity modulus 

E of the soil is calculated as the slope of the load settlement 
s 

curve given for the unreinforced case within elastic range as 2000 

t/m^. The computed reduction in settlement is 0.51 mm at a 

loading of 40 t/m* and 0.26 mm at a loading of 20 t/m*. These 

values correspond to an improvement of lOX and 9.6X respectively 

whereas the observed improvements at these loading levels are 

(Table 5.1b) around 20%. The reasons for the difference in the 

predicted and observed values arc possibly due to the following 

reason; (i) assumed value of E ; (ii) for strips with low flexural 

s 

rigidity the meiTibrane effect is also important at large vertical 
def or mat ions; and tiii) the depth of the sand bed is finite and 
hence the solution based on a semi -inf inite layer will be 
approximate . 

Fragaszy and Lawton (1984) conducted laboratory model tests 
to study the load -settlement behaviour of reinforced sand. The 
test set up details are given in Table 5.1c. In the dimensionless 
form, the parameters of the model tests are; footing aspect ratio, 
L^/B^=2, length of strips, L^/Bj=2.675, width of strips 
B /B=0.11, interval between strips = 0.45B . The depths are 

r f 

0.22B^, 0.45B^ and 0.67B^. The value due to shear stresses 

for the strips in three layers considered as single strips is 
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Table 5.1 Reduction in Surface Settlements; Comparison with 
Selected Results 


(a) Brown, B.B. and Poulos, H.B. , (1981) 


No. of 

layers 

q=16c 


q=20c 

Predicted 

Actual 

Predicted 

Actual 

1 

11.97. 

SI 37. 

137. 

SI 47. 

2 

137. 

SI 47. 

14.57. 

S167. 

(b) 

Shankar i ah. 

B. , and 

Narahar i. 

R., (1988) 

Footing Reinforcement 

Depth 

of E of Settlement Reduction 

Dimensions 

Details 

first 

Soil 

q=20t/m* q=40t/m* 


L B L B Spacing layer ^ Prcd. Act. Pred. Act. 

cm cm cm cm cm cm t/m 

8.5 8.5 43 0.5 3.5 2 2000 9.67. 207. 107. 207. 


tc) Fragaszy, R.J. , and Lawton, E. , (1984) 

Footing Reinforcement No. of Depth of 

Dimensions Details strips in Spacing layers 


LB LB 

cm cm cm cm 

Thickness 

cm 

a layer 
cm 

cm 

cm cm 

u 

a 

cm 

45.6 22.8 61 2.54 

0.0254 

9 

5.08 2.54 5.08 

7.62 


E of Soil 


Settlement 

Reduct ion 


KPa 

q=20KPa 

D =617. 

T* 

q=20KPa D =707. 

r 


D =617. D =707. 

r r 

Predicted Actual 

Predicted Actual 


1260 1850 

317. 

357. 

327. 417. 

i_ 


L = Length, B = Width, Pred. = Predicted, Act. = Actual 



0.0303. Multiplying this value by a factor of 0.6 to account for 


the extensibility, the net I is 0.0182. The I value due to 

sc sc 

normal stresses for all the strips is 4.1226. Choosing a 

correction factor of 0.2 for flexural deformations, the net I is 

sc 

0.825. Summing up the two I values and using an interference 

SC 

correction factor of 0.45 the total I is 0.379. At a loading of 

sc 

20 KPa and for 61% relative density the predicted reduction in 
settlements is 0.178 cm (31%) whereas the observed reduction is 
around 0.2 cm (^5%) (Fig. 5.33). Under the same loading and for 
a relative density of 70% the predicted reduction in settlements 
is 0.125 cm (S32%) as against the observed value of 0.16 cm 
(^1%). Thus there appears to be a reasonable agreement between 
the results of Fragaszy and Lawton and the reduction predicted 
from the present analysis. At larger relative densities the 
discrepancy increases. This is probably due to increased 
confinement action at greater relative densities which are not 
considered in the analysis. 

5.8 Conclusions 

The analysis developed in chapters 3 and 4 respectively for 
the interaction of embedded rigid (inextensible and inflexible) 
strips are extended for extensible and flexible ones. In the 
former case, i.e. extensible strip the strip is treated as a 
uniaxial member while in the latter case it is idealised as a 
beam. The effects of the elongation ratio, K^, for extensible 
strips and flexibility ratio, K^, for flexible strips on the 
stresses mobilised at the strip-soil interface and the consegueni 
reduction in surface settlements have been studied. The combinec 


settlement (% OF FOOTING WIDTH) 


Fla 5 32 Variation of SRC with Distance along Surface 
F.g, 5.32 Coefficient of Motion J^alues _ 

(L^/Bf=2, U/Bf=l. Uo/Bf=1, rBf/q-0'5* i^s-0.^) 


BEARING PRESSURE (kPa) 


BEARING PRESSURE (kPo) 



Frnnaszv Qod Lowton 'S (1984) 

Fig. 5.33 Comparison with^ Fragcs y 


Experimental Results 
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effect of considering simultaneously the rigidity of the strip in 
the horizontal and vertical directions has also been accounted for 
and provision made for possible slip at the interface. The 
settlement reduction of the points along the surface arc computed 
considering the combined effect and for cases where slip is 
accounted for. Conclusions drawn from the results obtained from 
the parametric study are listed below. 

1. The variations in the SRC for an extensible strip follow the 
same trend as those for a rigid one (chapter 31 and range from 
negligible values to those of a rigid one depending on the 
value of the elongation ratio, K^. 

2. The depth at which the SRC at the centre of the loaded area is 
a fflaximum, ranges from 0.75B^ to for ranging from very 
low values to that of a rigid strip. 

3. For a pair of strips placed symmetr ically about the x-axis, 
the distance of the centre of each strip from the x-axis for 
which the settlement reduction is maximum lies in the range of 
0. IBj^ to 0.25B^ for the ratio K^, varying from very low values 
to that of a rigid strip. This is applicable for both 
extensible and flexible strips. 

4. A suitable correction factor is obtained for the SRC at the 
centre of the loaded area for an extensible strip from the 
corresponding value for a rigid one depending on its length 
and the value of its ratios K (extensible) or K (flexible). 

T B 

5. For a flexible strip the variations in SRC follow the same 
trends as those for a rigid one (chapter 4). 

6. For the case where the coupled effect of rigidity of the strip 
in the longitudinal and normal directions is considered, the 



shear stresses reduce marginally whereas the normal stresses 
remain unaltered as compared to the case where the rigidity of 
the strip in both the directions are considered independently. 
Consequently the resulting values of SRC are the sums of the 
values observed for shear and normal stresses considered 
separately. 

7. The region along the strip where normal stresses mobilised are 
smallr fail due to slip. The settlement reduction coefficients 
are not significantly affected due to slip as the stresses get 
redistributed over the adjacent elements. 

8. The model is able to predict the settlement reductions due to 
strip reinforcements quite satisfactorily. 



CHAPTER 6 

SOIL-REINFORCEMENT SHEET INTERACTION : AXI-SYMMETRIC CASE 


6.1 Introduction 

The previous chapters dealt with the analysis of interaction 
between soil and long, narrow strip reinforcements beneath 

rectangular loaded areas. Another form of reinforcement commonly 
used in practice is reinforcement in the form of sheets. Sheet 
reinforcement is often used to design embankments over soft soils 
and to reinforce highway subgrades. In the latter case, the 
contact area between the wheel and the subgrade can be 
approximated by a circle. Thus, a simple approach to evaluate the 
influence of sheet reinforcement is to consider an axi -symmetric 
problem in which the applied surface load and the reinforcement 
are circular areas. Thus, in this chapter , the influence of a 
single circular sheet placed at some depth, on settlement 
reduction at the surface is evaluated. The sheet firstly is 
treated as rigid, transversely ( inextensibleJ and normally 

(inflexible). Subsequently, analyses of extensible and multiple 
layers of reinforcement beneath a uniformly loaded circular area, 
are presented. 

The BDUssinesq''s equation for radial displacements due to 
point load on the surface is integrated over the circular loaded 
area to obtain the radial displacements for points along the 
sheet. Hindi in"s equation is used to calculate the radial 
displacements at various points of the sheet due to the radial 





shear stresses mobilised. Satisfying the compatibility of radial 
displacements, and solving the resulting equations, yield the 
shear stresses mobilised at the soil -reinforcement interface. For 
an inextensible sheet, the net displacements are zero while for an 
extensible one they equal the radial elongations of the points 
along the sheet. In order to compute these elongations, the 
governing differential equation for a circular membrane subjected 
to tension is used. Following this analysis, rigidity of the 
sheet in the vertical direction, i.e. a rigid disc is considered. 
Normal stresses mobilised at the soil -sheet interface, as in the 
case for strips, are computed. The disc being rigid, experiences 
a uniform vertical displacement. Consequently the stresses act 
upward over the central portion and downward over the latter 
portion. The settlement reduction of points along the surface due 
to the stresses mobilised are then calculated. 

6.2 Horizontal Rigidity (Sheet Reinforcement) 

6-2.1 Problem Definition 

In this section, the rigidity of the sheet i.e. its 
inextensibility alone is considered. The problem consists of a 
circular sheet of reinforcement of radius, R, placed symmetrically 
below a circular loaded area of radius "a" transmitting a uniform 
load of intensity, q (Fig. 6.1). The surface load causes the soil 
to move downward and radially outward. The outward radial 
displacements of the points at the soil-sheet interface are 
resisted by the shear stresses mobilised. These stresses, in turn 

upward, thereby causing the surface to heave. 


push the soil above 


resulting in settlement reduction of points along the surface. 

fc.2-2 Rigid (Incxtcnsible) Sheet 

The radial displacement of a point within the elastic 

continuum with co-ordinates, (r,z) due to a point force on the 

surface is given by Eg. 3.1. In order to obtain the displacements 

of points along the soil -sheet interface at depth, U , Eg. 3.1 is 

o 

integrated numerically over the loaded area. For this purpose the 
loaded area is divided into n^ elements of length, Ar, along its 
radius and n^ elements of A© radians along the tangential 
direction. Since the loading is symmetric, integration is 
performed over half the loaded area only, i.e., over 180 degrees. 
The sheet is divided into N elements, each of length dr, along its 
radius and the shear stresses mobilised are assumed to be constant 
over each element. 


The global radial diplacement of the centre of any element, i 
on the sheet, due to the surface load is given by 
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I 


'■ f f * 


q ( 1 +1^ 3 
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2rTE R 
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rU 


(l-2v )r 
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(R+U ) 
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Cosa. . 
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( 6 . 1 ) 


where AA . is the area of element, j over the loaded area, ot . is 
J 

the angle between radial line joining element j and node i, and 
the axis along which the displacements are desired. The meanings 
of all other parameters are the same as in Eg. 3.1. 
Nondimensional ising all length parameters with the radius of the 
loaded area 'a' and calculating the radial displacements at all N 



nodes along the sheet r Eg* 6.1 is expressed in vector forift as 


where a coefficient vector for the influence of the 

surface load on the displacements of the nodes along the radius of 
the sheet. and (I'l are vectors of size N. 


The radial displacements of the N nodes due to the shear 
stresses mobilised along the radius, R, are computed using 
Hindi in •'s equation for horizontal displacements due to a 
horizontal force within the elastic continuum. This equation is 
integrated numerically over the radial and circumferential 
directions to obtain the displacements of all nodes due to shear 
stresses along the sheet. For this purpose, the area of the sheet 
is further divided along the circumferential directions into N' 
elements of d© radians each. The area of each element is then 
equal to r .d6dr where r. is the distance to the centre of each 

j j 

element from the centre of the sheet and dr is the length of each 
element. Fig. 6.2a shows the shear stresses over the sheet. Due 
to symmetry of the problem, the shear stresses are independent of 
©, the tangential direction, i.e. for 0*5 6 5 360 , but vary with 
the radial distance, r. 


It is to be noted here that the Mindlin's solution for a 
horizontal force beneath the surface is in terms of the cartesian 
co-ordinates. In order to calculate the radial displacements of 
the N nodes along the radius of the sheet, the components of the 




forces over each element are to be taken. The force on element j 
(Fig. 6.2b), is given by 


F. = T. r. dO dr 

The components in the x and y directions are 

P Cos© (6.4) 

and 

F . = F. Sin© (6.5) 

yj j 

where © is the angle made by the radial line joining the centre of 
the sheet and node J, with the x-axis. 

The radial displacements of the nodes are now calculated for 
the forces defined by Eqs. 6.4 and 6.5. The nodes are considered 
along the x-axis and the radial displacements of these points with 
respect to the global x-axis may as well be termed as horizontal 
displacements. The displacements are calculated at the points 
along the x-axis as shown in Fig- 6.2b. 


The horizontal displacements of the nodes due to the force. 


F , acting on any element j, is calculated from the Hindi in''s 

xj 

equation (Eq. 3.10). The displacement of node i due to F . is 
written as 



F . ( 1 ) x 

X J s 

SttE (1-v ) 
s s 



( 6 . 6 ) 



where x is the distance along the x-axis between nodes i and j 
while all other parameters are as defined in Eq. 3.10. 


Substituting Eqs. 6.3 and 6.4 in Eq. 6.6 one obtains 



T . ( l+v> )x 
j s 

SrrE 1 1 -v> ) 
s s 



r . Cosd d© dr 

J 


(6.7) 


In order to obtain the horizontal displacement of node i due 

to force, F , Hindi in^s equation for displacements along x-axis 
yj 

due to a force, Q, parallel to y-axis is used, and given as 


Q(l+w ) 
s 

P = X y 

** SttE (1-w ) 
s s 



where 

H 



t3-4v ) 
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R® R,(R,+z+c)® 

2 2 2 


( 6 . 8 ) 


R = ■/ x^ + y^+ (z-c)* 
1 

R = Vx^ + y^+ (z+c)^ 
2 


and X, y and z are the co-ordinates of the point where the 
displacement is to be calculated due to the horizontal force. 
Substituting for the force, in Eq. 6.8, the horizontal 

displacement of node i is 



F (l+v ) 
y J s 

8nE (1-v ) 
s s 


y 



(6.9) 




Substituting for F^. in terms of the stress, t.. 


one obtains 


Ty T , (1+p )xy 

'^xij “ InE (1-v ) 1^1'', (6.10) 

B B \ J J 


Eq. hmh, in dimensionless form is 
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KIJ 




T. d© dr 
J 


( 6 . 11 ) 


r 

where is a dimensionless coefficient for the influence of the 

force component, F . along x-axis due to the stress, t , on the 

j 

displacement of node i and •'a^ is the radius of the loaded area on 
the surface- Eq. 6.10, expressed in the dimensionless form is 


p I.^ T. d& dr 

XV J E VJ J 


( 6 . 12 ) 


r 

y 

where is a dimensionless coefficient for the influence of the 

force component, F along y-axis due to the stress, t., on the 

yj J 

displacement of node i. The net horizontal displacement of node i 
due to the stress on element j is obtained by summing up Eqs. 
6.11 and 6.12 as 


r a 

^ r r .. 

1 T. d© dr 

(6. 13) 

XVJ E^ 

L M '-j , 

1 J 



The displacement 

of node i 

due to stresses on 

all 

the 

elements 

along the radius 

and along 

the circumference 

of 

the 

sheet is 


obtained by integrating Eq. 6.13 over the area of the sheet as 



T. d© dr 

i 


(6.14) 


R 2 tt 


1 J t '‘i " ) 


0 o 


Due to symmetry the integration is carried out only over half the 
area of the sheet, as 


n 


1 t 




0 0 


1/ 1 T. d9 dr 
'•j J j 


(6. 15) 


Integrating numerically. Eg. 6.15 is expressed as 

r r 

pr-2^ ^ I [ 'u " ■** 

® j=l k=4 

where N^=N^/2, is the number of elements along half the circum- 
ference of the sheet. The horizontal displacements of all the N 
nodes in vector form is 

where ^ is a dimensionless coefficient matrix for the 

influence of the stresses on the displacements of the nodes and is 
a function of the depth ratio, U^/a , radius of the sheet, R/a and 
Poisson^'s ratio v . 

For a rigid sheet, the net radial displacements equal zero, i.e. 


{/} 


(6.18) 



Substituting for | | and | P*^ J one obta 


IDS 


[ ] { r/, } = { l' } 


(6. 19) 


Eq. 6.21 gives N equations, the solution of which gives the shear 
stresses mobilised along the reinforcement . 

The settlement reduction of points on the surface are 
evaluated using Mindlin''s equation for vertical displacement of a 
point within the elastic half-space due to a horizontal force 
within the continuum given by Eq. 3.26. In order to apply the 
equation, the force component of the stresses have to be resolved 
parallel to the axes corresponding to the cartesian co-ordinates 
CFig. 6.2b). The points on the surface where the settlement 
reductions are calculated are points with co-ordinates (x,0,0), 
i.e., they lie along along the global x-axis on the surface. The 
settlement reductions of these points alone are calculated. 


Substituting z—0 in Eq. 3.26 one obtains 


QCl+v )x 




(l-2v ) 
s 

RtR+c) 


where 


R = /x^ + y^ + 


y + c 


( 6 . 20 ) 


Eq. 6.20 gives the vertical displacement of the point on the 
surface due to a horizontal force, Q, at a depth, c. The vertical 
displacement of node k along the surface due to the force, 
acting on element j (Fig. 6.2b) at a depth is given by 


238 


K 
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jckj 


x= 


F . ( 1 ^ X, . 

X J S k j 



(1-2p ) . 

-] 

R<R+U ) ' 


( 6 . 21 ) 


2 2 

X, .+ y, .+ U , X, . and y. . are the distances between the 

*c j k j O kj kj 

centre of element j and node k along the x and y axis respective- 
ly. Substituting Eq. 6-3 and Eg. 6-4 in Eg- 6-21 one obtains 



T . 
J 


(1+1^ )X, . 
S k J 

SfrE 

s 



(1-2p ) 
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R(R+U ) 
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r. Cos© d© dr 
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( 6 . 22 ) 


Eq. 6.22 is expressed in dimensionless form as 



^ i;;. T. d© dr 

kj J 


(6.23) 


where G = E /2(l+v ). 
s s s 


The vertical displacement of node k due to the force acting on 
element j at a depth is given by 
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il+v 
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(6.24) 


where R=-/ x^ .+y* .+ U*. The terms, Xj^ . and Eq. 6.21 are 

interchanged in Eg. 6.24 since the force, F^. is parallel to the 
y-axis. Substituting Egs. 6.3 and 6.5 in Eg. 6.24., one obtains 



5;^^ 



I r. Sin© d© dr 

R(R+U^) ^ ^ 

o 


(6.25) 


Eg. 6.25 is expressed in dimensionless form as 



t6.26) 


*kj b kj j 


The net vertical displacement of node k due to shear stress t , is 

j 

obtained by summing up Eqs. 6-23 and 6.26 as 




a 

“F" 
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d© dr 


(6.271 


The vertical displacement of node k, due to shear stresses on all 
the elements on the sheet, is obtained by integrating Eq. 6.27 
numerically over the area of the reinforcement as 


M N 


e 






d© dr 


(6.281 


j=i t=i 


The integration is performed over half the area only because of 
symmetry. The vertical displacements of all the nodes along the 
surface in vector form is 




(6.291 


The settlement reduction is desired at points along the radius of 
the loaded area, which is divided into points. Hence ^ 

a vector of size and j is a matrix of size (N^xNl. 

Eq. 6.29 is rewritten as 




q 


(6.301 
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I j is a Settlement Reduction Coefficient vector of size N 


which is obtained as 


{■•}■[■;]{') 


(6.31) 


The vector of SRC along the surface is 



(6.32) 


From Eq. 6.30, the reduction of settlements of points along the 

surface due to a circular reinforcement sheet at depth, U , is 

o 

obtained. 


6.2.3 Extensible Sheet 

The net radial displacements of the nodes for the rigid 
( inextensible) sheet, as seen in section 6.2.2 were equal to zero. 
For an extensible sheet the net radial displacements due to the 
surface load and the shear stresses are equal to the elongations 
of the points on the circular sheet in the radial direction. In 
order to compute the elongation of the sheet, the same is 
considered to be a thin circular membrane. The thickness of the 
sheet is negligibly small in comparison to its radius. 

Fig. 6.3 shows the forces on an element of the circular 
membrane in polar co-ordinates. Satisfying the force equilibrium 
of an element, two equations of equilibrium in the radial and 
tangential directions are arrived at, as 



r 

w- 


(6.33) 


1 &T 


and 


1 9a 


W 
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+ F = 0 
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_ + __ + = 0 


(6.34) 


where F^ and F^ are the body forces in the radial and tangential 
directions respectively. 


In the case of circular reinforcement, i.e. in an axi- 
symmetric case, a^ and vanish and a^ is independent of &, due 
to symmetry. The term F can be substituted by the force per unit 

NT 

volume on the element due to the mobilised shear stress, r. The 
body force, F in the tangential direction is zero. Eg. 6.34 is 

T 

identically satisfied. Eq. 6.33 reduces to 



m 


where t is the thickness of the sheet. 


The strain components in the radial and tangential directions in 
terms of the displacement or the elongation, u, are 


s 

r 


du 

dr" 


and 



(6.36) 


Stress components expressed in terms of the strain components are 
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E E ^ 
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(6.37) 


where v is the Poisson^s ratio of the reinforcing sheet, and E^, 

m 


its the elastic modulus. Substituting Eqs. 6.36 and 6.37 in 
Eq. 6.35 the governing differential equation is obtained as 

^d^u du r*(v>^-l)T 

r — _ + r — - u = p-? (6.38) 

dr dr m 


Eq. 6.38 is expressed in finite difference form for the N nodes 
along the radius of the sheet. The equation for node i written in 
finite difference form (Fig. 6.4) is 
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(6.39) 


Eq. 6.39 written in matrix form is 

m 

where ^ finite difference coefficient matrix for the 

extensible sheet. 


The boundary conditions are 

At r = Or elongation, u = 0 (6.41) 

At r = R, the radial stress is zero, i.e-jcr = 0 (6.42) 

r 

r=0 and r=R correspond to the nodes, O and N+1 (Fig. 6.4). These 
nodes are at a distance of dr/2 from their adjacent nodes while 
all the other N nodes are equidistant ly spaced at a distance, dr. 
Hence, second degree polynomials in u are fitted through the first 







and last three nodes such that they satisfy the boundary 
conditions. The constants of the polynomials are obtained in 
terms of elongations of the three nodes. The first and last rows 
of the coefficient matrix, ^ j thus obtained. 

The net soil displacements from section 6.2.2 are 




(6.43) 


The sheet displacements as expressed by Eq. 6.40 are 




a(v -1) 

Tn 

_ 


(6.44) 


The compat ibi 1 i t y of displacements of all nodes is satisfied by 
equating the net soil and sheet displacements* Substituting Eq- 
6.43 in Eq- 6.44, one obtains 




(6.45) 


where 




E t 

is an identity matrix and K^= relative 

elongation ratio of the reinforcement sheet. Eq. 6.45 gives N 
equations for the shear stresses, r/q ^ mobilised at the soil- 

sheet interface which are solved for the same. 


The SRC of points along the surface are calculated as 


described in section 6.2.2. 



6.3 Vertical Rigidity (Rigid Inflexible Disci 

The problem definition is the same as in section 6.2.1 with 
the difference that the rigidity of the sheet in the vertical 
direction alone is considered. The disc is assumed to undergo 
uniform vertical translation due to the loading on the surface and 
as a result normal stresses are mobilised at soil-sheet interface. 
The discretisation of the loaded area and the reinforcing sheet 
are the same as explained in section 6.2.2. 


The vertical displacements of the nodes along the sheet are 
calculated using Boussinesq's equation defined by Eq. 4.1. The 
displacement of node i at a radial distance, r. on the sheet due 

X 

to the surface load is obtained by integrating Eq. 4.1 over the 
loaded area as 
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where the terms n , n^, and AA are as defined in section 6.2.2 

f © j 

and other parameters are as defined in section 4.2.2. Eq. 6.46 
expressed in vector form is 

a coefficient vector for the influence of the 
surface load on the vertical displacements of the nodes along the 


radius. 



The displacements of the N nodes due to the normal stresses 
mobilised at the interface are computed using Eq. 4-6. The 
diaplaccmcnt of node i due to the stresses over all the elements 
is calculated by integrating Eq. 4.6 along the radial and 
tangential directions of the sheet as 
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Performing the integration numerically and over half the area 
because of symmetry, Eq. 6-48 is written as 
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where is the number of elements along half the circumference of 
the sheet. The vertical displacements of all the N nodes in 
vector form are 


where [ j is a dimensionless coefficient matrix for the 

influence of the stresses on the vertical displacements of the 

nodes and is a function of the depth, U /a, radius of the sheet, 

o 

R/a, and Poisson^s ratio of the soil, v . 

For a rigid sheet, the net vertical displacements equal the 
rigid body translation, 6^ of the sheet, which is written as 



Substituting Egs. 6.47 and 6-50 in Eq. 6.51, one obtains 
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(6.52) 


Eq- 6-52 gives N equations, while there are N unknowns in o' J 

and a rigid body translation- The additional equation is obtained 
by satisfying the force equilibrium along the radius of the sheet 
which is 
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Eqs 6.52 and 6.53 are rewritten as 
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(6.55) 


(6.55) 


Eqs. 6.54 and 6.55 are solved for the normal stresses, o/q and the 
rigid body translation, <5^. 


These normal stresses push the soil upward, thus bringing 
about a reduction in settlements of points along the surface. In 
order to compute these reductions in settlements Eg. 4.20 is 
integrated over the area of the circular disc. The vertical 



displacement of node k, on the surface due to ths mobilised normal 
stresses on the sheet is given by 
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The vertical displacements of all the nodes on the surface in 
the vector form is 
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(6.57) 


(N xN). 
f 


is a vector of size and 
Eq. 6.56 is rewritten as 



is a matrix of 


size 




(6.58) 


I is the SRC vector of size N^- From Eq- 6.58, the reduction 

of settlements of points along the surface due to a circular 

r einf orccment disc at depth U experiencing vertical displacement 

o 


alone, is obtained. 



6.4 Results and Discussion 


6.4.1 Rigid Sheet and Disc 

It is intended to evaluate the shear and normal stresses 
mobilised at the soi 1 -reinforcement interface along the radius of 
the circular reinforcement. The Boussinesq^s equation is 
integrated numerically over the cicular loaded area to obtain the 
radial displacements in the case of shear stresses and vertical 
displacements in the case of normal stresses. The area is 
discretized into n^ subelements of length, 0.02a, each along the 
radius and n^ subelements of 0.4 degrees each along the 
tangential direction. Further discretization did not improve the 
solution. Due to symmetry, the integration is performed over 
half the loaded area. The displacements are obtained at the 
centre of the N elements along the radius of the reinforcement. 

In order to compute the influence coefficients for the 
displacements of points along the radius of the reinforcement due 
to the stresses mobilised at the interface, the reinforcement is 

divided in the tangential direction into elements of 2 degrees 

o 

each (A0=O.O349 ). The radius is divided in such a way that the 

area of each element remains the same. In doing so, the length 

of each element j, CAr .) decreases with increase in the distance 

1 

from the centre. If the number of elements along the radius are 
large, the centre of the last element gets located very close to 
the edge. For eg., if a sheet of radius, R/a=l is divided into 
10 elements, the centre of the lO*"^ element lies at a distance of 
0.973a from the centre of the sheet. For a rigid reinforcement 



the edge stresses increase rapidly (as is evident in the strip 
case) and hence care is taken to position the centre of the last 
element sufficiently away from the edge to prevent the stress 
there from being infinitely large. It was seen that for R/a<4, 5 
elements over a distance of r/a=l were sufficient and for R/a=5r 
6 elements over the same distance were required. Hence, the 
number of elements equal 5xR/a for sheets upto radius, R/a=4 and 
30 for sheets with R/a=5. The discretization is carried out such 
that no further improvement in the influence coefficient due to 
stress on element j, on the displacement of any node i, is 
observed. The stress over each element is assumed to be constant. 
For the aforementioned discretization the length of the first 
element is comparatively large and the centroid of the element 
lies at a considerable distance from the centre. It hence becomes 
imperative to extrapolate the stresses over a distance from the 
centre to the centroid of the first element. The variation in 
normal stresses over this region may not be significant (as seen 
for the strip case) and the assumption that stress over the 
element is constant stands valid. On the other hand, shear 
stresses increase from zero at the centre to the calculated value 
at the centroid of the element. The assumption, hence does not 
stand to reason atleast for the first element. To solve for the 
stresses within this region, element 1 is discretized further 
into subelcments and the displacements of these subelements due 
to the surface load, due to the stresses over these subelements 
and those computed for the other elements, are calculated. The 
stresses are evaluated by solving the equations resulting from 
the displacement compatibility. 



A par amp t r i c study enumerates the effects of the depth, 

U /a, and the radius, R/a, of the reinforcement on the stresses 

0 

mobilised and the settlement reduction resulting thereof. The 
first few figures illustrate the effects of the parameters on the 
shear stresses mobilised at the interface and the resulting 
settlement reductions for a rigid sheet. All results are 
obtained for a Poisson's ratio of the soil, v^=0.3. 

Fig. 6.5 illustrates the numerical instability with 
elements of equal lengths. The figure plots the variation of the 
normalised shear stress, r/q, with distance, r/a, for different 
radii of the sheet placed at a depth, U^/a. The instability is 
reflected on the stresses mobilised along the elements close to 
the centre of the sheet. Elements of equal lengths imply that 
the area of elements increase with distance along the radius of 
the sheet- As a result, the influences of shear stresses on 
elements farther away from the centre, on the displacements of 
nodes close to the centre are very high. Moreover, these 
influences are much higher than the influence of the stresses on 
the first few elements, on the displacements at these nodes. The 
stresses mobilised along the first few elements, hence show an 
erratic behaviour. With increase in the distance from the 
centre, this instability reduces. 

Fig. 6.6 shows the variation of the normalised shear 
stress, T/q, with distance, r/a, along the radius of the sheet 
for various depths of placement, U^/a, and for a sheet of radius, 
R/a=2. The trends are similar to those of the strip reinforcement 
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Rg. 6.5 Variation of Shear Stress with Distance 

for Various Radii of Sheet (Uo/a = 1, i/s=0.3) 


Fig. 6.6 Variation of Shear Stress with Distance 

for Various Depths of Placement of Sheet 
(R/a = 2. i/s=0.3) 




case- Shear stresses increase from zero at the centre (r/a=0) 
and reach a maximum value at a distance of around 0.8a to ' a ' 
(i.e.r below the edge of the loaded areal. At shallow depthSr 
U^/a=0-25, the stresses are positive only over distances in the 
range, r/a=0.<lE> to 1.2. A maximum stress of 0.22q is observed at 
a distance of r/a=l for the sheet placed at a depth of 0.5a. As 
the depth increases the stresses decrease and become uniform over 
the radius of the sheet- The maximum stresses occur ing below the 
edge of the loaded area for the sheet at shallow to intermediate 
depths (0.25<U^/a^ 1) is due to the fact that maximum radial 
displacements occur at this distance. 

The variation of the normalised shear stress, r/q, with the 

distance along the radius, r/a, for various radii, R/a, of the 

sheet and placed at a depth, U /a=l, is depicted in Fig. 6.7. It 

o 

is evident that for sheets of radius equal to that of the loaded 
area the stress increases monotonically and at the edge it 
reaches a high value (T/q=0.58). As the radii of the sheets 
increase, the positive shear stresses remain unaltered with a 
maximum stress of around 0.2q occur ing below the edge of the 
loaded area (r/a=l). For greater radii of the sheets tR/a>2) the 
stresses change sign beyond a distance of around 2.5a to 2.7a. 
This is because the radial displacements beyond this distance are 
negative and as such the shear stresses are directed outward. 
The stresses at the edge for sheets of radii, R/a>2, drop down 
sharply to about -O.OSq for R/a=3 and about -0.06q for R/a=5. 


The effect of the shear stresses in reducing surface 



5 et tlcuicnts is studied in the subsequent figures. Fig. 6.8 shows 
the variation of SRC with distance along the radius of the loaded 
area on the surface for various radii of the reinforcement placed 
at a depth of la. The maximum SRC is observed at the centre of 
the loaded area which decreases with the distancer r^/a. The 
increase in SRC with increase in the radius, R/a, of the sheet 
from 1 to 2 is significant. There is a marginal decrease in the 
SRC with increase in the radius, R upto 4a. This decrease is due 
to the negative stresses mobilised at the interface beyond a 
distance of around 2- 6a along the radius of the sheet. The SRC 
at the centre of the loaded area for sheets with R/a=l is 0.033 
while for sheets with R/a=2 it is about 0.05. 

The variation of SRC with distance for various depths of 

placement of the sheet of radius, R/a=2, is presented in Fig. 

6.9. Similar to the results for the strip case, sheets placed at 

very shallow depths (U /a=0.25) show negative SRC, i.e. the 

o 

points on the surface undergo additional settlement. This result 
indicates that sheets of intermediate radii tR/a£2) placed at 
shallow depths are detrimental to the foundation. The maximum SRC 
is observed for sheets placed at a depth of 0.75a to 1.0a. As the 
depth increases, the SRC decrease and become uniform over the 
distance along the surface. The maximum SRC at the centre is of 
the order of 0.05 for sheets placed at a depth of 0-75a to 1.0a. 

Fig. 6.10 shows the variation of the SRC at the centre, 

I , with the depth, U /a, for various radii of the reinforce- 
»c o 

ment . For all radii the I is a maximum at a depth of about 

sc 
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Fig. 6.7 Variation of Shear Stress with Distance 

for Various Radii of Sheet (Uo/a= 1 , 1 / 5 = 0 . 3) 


Rg. 6.8 Variation of SRC with Distance along Surface 
far Various Radii of Sheet (Uo/a=1, 1 / 5 = 0 . 3) 






0.75a. The a shallow dfjpth, U^/a=0.25, is negative. It 

increases rapidly with increase in depth and reaches a maximum 

value at about U =0.75a. Beyond depths of 0-75a, I reduces 

^ IS c 

gradually with U . The rapid increase in I from shallow to 

o sc 

intermediate depths (0.25<U /a:<0.75) can be attributed to the 

o 

change of stresses from negative to positive in this depth range 

and the mobilisation of maximum positive stresses at a depth, 

U SO. 75a- A maximum I of about 0.057 is observed at a depth of 
o sc 

0.75a for a sheet of radius 5a. For sheets with radius equal to 
that of the loaded area, i.e., R/a=l, the maximum I is 0.038 

SC 

and is observed at a depth, U =0.75. At greater depths (U /a=2) 

o o 

the increase in I with increase in radius of the reinforcement 

sc 

is not appreciable for sheets with radius, R>2a. 

On having looked into the variation of shear stresses along 
the radius of the sheet and its effect on the settlement 
reduction of points on the surface, the rigidity of the 
reinforcement in the vertical direction is considered in the 
subsequent figures. The effect of the parameters mentioned 
earlier on the normal stresses mobilised at the interface of the 
reinforcement disc and the soil are studied. The reduction of 
settlements of points along the surface as a consequence of the 
mobilised normal stresses arc also evaluated. 

Fig. 6.11 presents the variation of the normalised normal 
stress, cx/q, with the distance, r/a, along the radius of the disc 
for various depths of placement of the disc of radius, R/a=2. 
The variation is similar to the trend observed for strip form of 


reinforcement 


he normal stress decreases from a maximum at the 


centre » changes sign and approaches a large negative value at the 
edge of the disc. Discs placed at shallow depths show maximum 
normal stresses due to high vertical displacements due to the 
loading. As the depth increases, stresses decrease and tend to 
become uniform. This is due to decrease in vertical displacements 
of points at depth resulting from the surface load. 

The effect of the radius of the reinforcement disc, R/a, on 

the normalised normal stresses, o'/g, for the disc placed at a 

depth, U /a=l, is presented in Fig- 6.12. For a disc of radius 
o 

equal to that of the loaded area the normal stress near the 
centre of the disc is of the order of 0.52q. The stress 
decreases and becomes negative at a distance of 0.8a. Due to the 
edge effect, the stress near the edge is around -0.92q. As the 
radius of the disc increases, the positive stresses do not show 
an appreciable change. For any radius greater than 2a, the 
stresses become negative beyond a distance, r/a=1.8. The maximum 
negative stress at the edge decreases with increase in the radius 
of the disc. The maximum stress near the centre of the disc is 
of the order of 0.68q for R/a5:l. 

Fig, 6.13 shows the variation of SRC with the distance 
along the surface for various radii, R/a, of the disc placed at a 

depth, U /a=l. The SRC at the centre of the loaded area is a 

o 

maximum for all radii of the disc. The SRC increases with 
increase in the radius of the disc. The positive stresses 
increase with increase in the radius from R/a=i to 2. As a 
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Fig. 6-11 Variation of Normal Stress with Distonce 
for Various Depths of Placement of Disc 
(R/a=2, Vs=0.3) 
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Fig. 6.12 Variation of Normal Stress with Distance 
for Various Radii of Disc (Uo/a=1, 



result the increase in SRC is significant with the radius, R/a, 
increasing from 1 to 2. With further increase in the radius the 
rate of increase in SRC decreases. This is because the increase 
in stresses with increase in the radius is marginal. 

The effect of the depth of placement of the disc, U /a, on 

o 

the SRC of points on the surface for a disc with a radius, R/a=2 

is presented in Fig. 6.14. The SRC at the centre of the loaded 

area, again is a maximum for all depths. The stresses being high 

for the disc placed at a shallow depth CU^/a=0.25), the SRC at 

that depth is high. The variation in SRC with the distance along 

the surface, r^/a, is high for shallow depths. As the depth 

increases the SRC tend to become uniform. This is because the 

stresses at greater depths tend to become uniform. The SRC at 

the centre is about 0.385 for a disc placed at a depth, 

U /a=0.25, while that for a disc at a greater depth (U /a=2), is 
o o 

0.02. The decrease in SRC with depth is due to the decrease in 
stresses with depth. 

Fig. 6.15 depicts the variation of the settlement reduction 
coefficient at the centre, I , of the loaded area with depth, 

SC 

U /a, for various radii of the disc. As expected, the I 
o sc 

decreases with depth, the disc with larger radii showing higher 

1 values. A maximum I of 0.58 is witnessed for a disc of 

sc sc 

radius, R/a=5 at a depth, U /a=0.25. The rate of increase in I 

o sc 

decreases with an increase in R/a. This, again is due to the 
fact that the increase in stresses with increase in R/a from 1 to 

2 is significant while the rate of increase in stresses beyond 





this radius decreases with R/a. 

In Fig. 6-16 the c-ffect of the depth of placement of the 

reinforcement disc, U /a, on the uniform translation, 6 E /ag, of 

o os 

the disc for various radii of the disc is presented. The uniform 

vertical translation decreases with increase in depth and radius 

of the disc. When the radius of the disc is small, the 

difference between the vertical displacements at the centre and 

the edge of the disc due to loading is small while for a larger 

radius this difference is large. In order to maintain a uniform 

displacement the disc with a smaller radius undergoes a larger 

uniform translation as compared to the one with a larger radius. 

As the radius increases, the rate of decrease in the uniform 

translation, 6 E /aq, with depth also decreases. At a depth, 
o s 

U^/a=0.25, the uniform vertical displacement of the disc with 

radius, R/a=l, is 1.38 while for a larger radius, R/a=5, it is 

0.28. At a depth, U /a=2, the corresponding values are 0.8 and 

o 

0.27 respectively. 

From these figures it is clear that the settlement 
reductions due to shear stresses are much smaller than those 
resulting from normal stress mobilisation. It is advantageous to 
place a circular reinforcement rigid in the vertical direction 
at shallow depths while a reinforcement highly flexible in the 
vertical direction but rigid axially is best utilised if placed 
at a depth, U =0.75a to la, where ' a ' is the radius of the loaded 


area. 
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6.4.2 Extensible Sheet 

As seen in Fig. 6.5, equal lengths of elements result in 
erratic values of shear stresses along the first three elements. 
This was the reason why elements of equal areas were selected for 
rigid ( inextensible) sheets and inflexible discs. For an 
extensible sheet, the elongations are expressed as a differential 
equation which are solved using the finite difference technique. 
The formulation for the finite difference calculations are based 
on equal lengths of the elements. Hence, shear stresses are 
computed using elements of equal lengths. It is known that the 
shear stresses should vary smoothly starting from zero at the 
centre of the sheet. To smoothen the stresses along the first 
three elements, a quadratic curve for shear stresses is fitted 
between the centre of the sheet, and elements 4 and 5. 

The elongation ratio, K^, is varied from 10 ‘ to 10*. For 
K >10* there was no difference in the values of shear stresses 

T 

and the settlement reduction coefficients as compared to those 
obtained for K =10*. Hence, the upper limit of was restricted 

T 

to 10*. 


Fig. 6.17 depicts the variation of the normalised shear 
stress, T/q, along the radius of the sheet for different values 
of the elongation ratio, K^, and for a sheet of radius, R/a=2, 
placed at depth U^/a=l. It is evident that as the value of 
incr=as^s, thE stresses also increase and tend to those of a 
rigid sheet for K^=10‘. Stresses, for a highly extensible sheet 
(K =l0-‘) are very small owing to large elongations of the sheet. 
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For sheets with K <1 the ^treceoir j 

^ tne stresses at the edges are negative. The 

strBESES are maximum at a distance equal to the radius of the 

loaded area, and a maximum stress, r/q, of around 0.2 is observed 

for sheets with K^=lo‘ which equals that of a rigid sheet 
(T/q=0.2). 


The variation of the normalised shear stresses, T/q, with 

distance, r/a, for different depths of placement, U /a, and for a 

o 

sheet of radius, R/a=2 and K^=l is shown in Fig. 6.18. The 

variations in shear stresses follow the same trend as those for a 

rigid sheet. The maximum shear stress is observed below the edge 

of the loaded area for sheets at a depth, U^/a=0.5 which equals 

0.165. For sheets at shallow depths (U /a=0.25) stresses over a 

o 

distance, r/a, of 0.75 to 1.25 are positive with a maximum 
positive stress, r/q, of 0.15 occuring below the edge of the 
loaded area. 

Fig. 6.19 shows the effect of depth, U /a, on the 

o 

settlement reduction coefficient, I , at the centre of the 

sc 

loaded area for different values of the elongation ratio, K , and 

T 

for sheets of radius, R/a=2. It is seen that for all values a 

maximum I occurs in the depth range of U /a=0.65 to 0.8. I 

sc o sc 

values increase with increase in K values owing to the increase 

T 

in shear stresses with K . Isc values for sheets with K =10 

T * 

tend to those of a rigid sheet (.Fig. 6.9). At shallow depths 

(U /a=0.25), I values are negative since the stresses are 

o BC 

predominantly negative. 



tonce 

:heet 
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The variation of I with radius, R/a, of the sheet for 

SC 

different K values and for the sheet at a depth, U /a=l is 
presented in Fig. 6.20. For K ^1, I values gradually increase 

T SC 

with R/a increasing upto 2 and with further increase in the 
radius of the sheet the I values are almost constant. For K >1 

SC T 

there is a small decrease in the I values with R/a increasing 

sc 

from 2 to 3. With further increase in R/a, I „ increases 
gradually. This trend is prominent for sufficiently rigid to 
inextensiblc sheets and can be explained in conjunction with the 
results presented in Figs. 6.7, 6.8 and 6.10. In addition to the 
comments made with reference to Fig. 6.8, it may be noted from 
6.10 that I values attain their maximum values at different 

SC 

depths for sheets of different radii. For a sheet with R/a=2, 
the maximum I is for a depth, U /a of 0.85 while for other 

SC ^ 

combinations, I is maximum at U /a=0.75. The slight reduction 

Sc 

in I with increasing values of R/a beyond R/a=2 is due to the 

SC 

increasing stretch of negative stresses tFig. 6.7) over the sheet 
for r/a>2.0. 


The varation of the ratio of the for an extensible 

sheet I > to that for a rigid linextensible) one, ^BC<r>’ 

' SC<E> 

with the elongation ratio, K^, for various radii, R/a, of the 
sheet placed at a depth, U^/a=l is depicted in Fig. 6-21. The 

ratio, I /I , increases with K, for all radii of the 

r d L A setr) * 

. . Crtv w =10* Sheets with R/a=5 show 

sheet and approaches 1 for 10 . 

, ~ vat in I /I , as compared 

marginally higher values of the ratio, sc<r) 

to sheets with R/a=l, while sheets with R/a=2 show the least 
values of the ratio. The reason for this behaviour is due to the 




Fiq. 6.19 Effect of Depth of Plocement of Sheet on Isc 
for Various clongction Ratios (R/a—2, j/,— 0.3} 



Rq. 6.20 Effect of Radius- of Sheet on Ijc :fo'' Various 
Elongation Ratios (Uo/a— 1, i/s— 0-3) 
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decrease in values with increase in the radius, R/a, of the 
rigid sheet f.om 2 to 3 (Fig. 6.20) and follows from the above 
discussion. For K >10 the value of the ratio, I /i for 
sheets of all radii is greater than 0.9 which means that the 
sheets become almost rigid. This figure can be used to calculate 
the settlement reduction coefficients due to an extensible sheet 
placed at a depth, U^/a=l, from those of a rigid one. 


The variation of the ratio, I /I , with K for 

SC(E> SC<r> T 

different depths of placement, U^/a, of the sheet of radius, 
R/a=2 is presented in Fig. 6.22. With increase in the value of 
the elongation ratio, K , the ratio, I /I , increases 

T SC<E> SC<r> 

and approaches 1 for large values of K . It is seen that for the 

T 

sheet placed at a depth, U^/a=0.5, the values of the ratio, 

maximum while they decrease with 

HdE/ mG\Tf 

increasing values of U /a and are the least for the sheet at a 

o 

depth, U /a=2. Results depicted in this figure are the same data 
o 

shown in Fig. 6.19 but with I values normalised with respect to 

SC 

the I values for an inextensible sheet (K =1000). The ratio 
sc T 

I /I reflects the relative variations of the coefficient 

SC<E> SC<r> 

I with the parameters, U /a and K . It may be noted that the 

SC o T 

rate of change of I with depth, for rigid sheets is high as 

SC 

compared to extensible ones. The increase in I for a rigid 

SC 

sheet from U /a=0.5 to 1.0 is sharp as compared to that for a 
o 

sheet with K =1 in the same depth range. As a result the ratio, 

T 

I /I , is high for U /a=0.5 and decreases with depth. 

SC<E> SC<r> O 

From this figure the settlement reduction coefficients for 
extensible sheets of radius, R/a=2, placed at different depths 
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can be calculated from those values of a rigid sheet. 

6.5 Conclusions 

Based on the elastic continuum approach a model has been 
proposed to predict the reduction in settlements of points along 
the surface due to the presence of a circular reinforcement at 
depth. The problem considered is axi -symmetr ic in nature. The 
effects of the radius of the reinforcement, and depth of its 
placement on the shear and normal stresses mobilised along the 
r e inforcement -soi 1 interface and the subsequent reduction in 
surface settlements are studied. The horizontal and vertical 
rigidity of the reinforcement are considered separately. 
Subsequently, the extensibility of the sheet is also accounted 
for. The results of the parametric study are listed below; 

1. Rigid (inextensible) reinforcement sheets of large size show 
maximum settlement reduction in the depth range 0.6 to 0.8 
times the radius of the loaded area. 

2. In the case of rigid (inflexible) reinforcement discs, maximum 
benefit is accrued when the discs of large radii are placed 
close to the ground surface. 

3. Extensible sheets exhibit trends similar to those of an 
inextensible one. Based on the radius, depth of placement and 
elongation ratio of extensible sheets the settlement reduction 
coefficients can be obtained from those of a rigid 

(inextensible) one. 



CHAPTER 7 

SUMMARY AND CONCLUSIONS 


Using the elastic continuum approach an analytical model to 
study the interactions between reinforcements and soil, and to 
predict the reduction in surface settlements due to reinforcements 
below a loaded area, at depth is presented. Two forms of 
reinforcements are considered, viz., strip and sheet or disc 
reinforcements. The mechanisms considered are the shear and 
normal stress interactions. Initially, each of the two mechanisms 
are considered independently. In the case of strip reinforce- 
ments, rigid strips of sufficiently small widths are considered so 
that the stresses vary along the length of the strip alone. The 
compatibility of horizontal and vertical displacements at the 
interface of the strips, considered separately gives rise to 
equations which on solving yield the shear and normal stresses 
mobilised at the interface of the strips. The reductions in 
surface settlements, as a consequence of these stresses are 
evaluated. The analyses consider a single strip, a pair of strips 
at the same depth, two strips one below the other and multiple 
strips at the same depth. In the case of strips one below the 
other and multiple strips the validity of the principle of 
superposition is examined. 

Following the analysis for rigid strips, the same is 
extended for extensible and flexible ones. In these cases the 
elongations and flexural deformations of the strips are considered 
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separately. The equations resulting from the displacement 
compatibility and force equilibrium relationships are solved to 
shear and normal stresses mobilised at the interface. 
The reductions in surface settlements due to these stresses are 

computed. The coupled effect of the horizontal and vertical 

interactions of rigid strips is examined subsequently. In this 
analysis the influence of horizontal stresses on vertical 

displacements and vice-versa are considered. Slip at the interface 
for rigid strips resulting due to the mobilisation of shear 
stresses greater than the available strength is accounted for. 

The analyses for strips is extended for sheet reinforcements 
beneath circular loaded areas taxi-symmetr ic case). Rigid tin- 
extensible) and extensible sheets are considered. The shear 
stresses mobilised at the interface and consequently, the 
reduction in surface settlements are evaluated. Rigid discs 

experiencing uniform vertical displacements arc also considered. 

Based on the results of the analyses for all these cases the 
following general or broad conclusions are drawn: 

1. The normal stresses resulting from vertical stress inter- 
actions for the strips as well as the sheets (axi -symmetric 
case) are much higher than the shear stresses resulting from 
the horizontal stress interactions. As a result, the 

Settlement Reduction Coefficients (SRC) due to normal stress 
interactions are much higher than those due to shear stress 


2 . 


interactions. 

Rigid (inflexible) strips of lengths (2B^ to SB^.) placed at 
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shallow depths reduce settlements significantly. 

3. The optimum location for placing a pair of strips is to place 

them at a distance, S /B=0.25 from the centre of the loaded 

y f 

area. 

4. Slip along the rigid strip occurs at points where the normal 
stresses mobilised are low. Relatively more number of points 
fail on strips at greater depths due to slip, since normal 
stresses are uniformly distributed over it. 

5. Highly flexible but inextensibile strips or sheets if placed 
in the depth range of 0.5 to 1 times the half-width of the 
loaded area provide maximum settlement reduction. The optimum 
length of the strips is 2 to 2.5 times the width of the loaded 
area. 

6. The same conclusions apply for circular reinforcements also. 

The above study can be extended by considering 

(i) the elasto-plastic response of the soil; 

(ii) the membrane effect at relatively large displacements; 

(iii) grid type reinforcement; 

(iv) the elasto-plastic type deformation of the interface and of 
the strips or sheets; 

(V) the soil encapsulated by reinforcement as a relatively stiff 
or rigid inclusion, etc. 
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